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ABSTRACT 


An  analytical  computer  model  has  been  developed  for  the 
analysis  of  reinforced  or  prestressed  concrete  multi-celled  box  girders 
acted  upon  by  loading  combinations  comprised  of  torque,  bending  moment, 
and  shear .  Any  cross-sectional  geometry  defined  by  linear  segments  can 
be  accommodated,  and  generality  of  loading  or  boundary  conditions  is 
ensured  by  the  inherent  characteristics  of  finite  element  modelling. 

The  analytical  model  features  the  incorporation  of  non-linear  material 
behaviour,  aggregate  interlock,  dowel  action,  and  the  warping  restraint 
and  cross-sectional  distortional  stiffness  of  diaphragms. 

To  assist  in  the  evaluation  of  the  analytical  model  performance, 
seven  double-celled  prestressed  concrete  box  girders  were  cast  and 
tested,  five  beams  having  a  rectangular  cross-section  and  the  remaining 
two  a  trapezoidal  cross-section.  All  seven  beams  were  subjected  to 
various  torque,  bending  moment,  and  shear  load  combinations. 

Performance  of  the  computer  model  was  assessed  through  compar¬ 
ison  with  experimental  and  current  theoretical  results.  The  satisfactory 
outcome  of  the  assessment  verified  the  value  of  the  analytical  model  as 
a  flexible,  sophisticated  method  of  analysis. 
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CHAPTER  I 


INTRODUCTION 


1 . 1  Need  for  Research 

Cellular  construction  in  reinforced  concrete  is  commonplace  in 
current  civil  engineering  practice  as  this  method  of  construction  is 
both  functional  and  economically  competitive.  In  addition,  the  flex¬ 
ibility  permitted  by  the  cast-in-situ  technique  enables  box  girders  to 
assume  any  desired  alignment,  as  illustrated  by  the  varied  configura¬ 
tions  of  the  large  number  of  skew  and  curved  highway  bridges.  Where 
large  structures  are  exposed  to  public  scrutiny,  the  aesthetic  appearance 
of  box  girder  design  is  an  especially  valuable  asset. 

As  the  application  of  concrete  box  girder  design  becomes  more 
diverse,  greater  demands  are  made  of  the  civil  engineer  in  the  design 
of  structures  of  increasing  complexity.  Thus,  the  incentive  for  research 
in  the  obscure  fields  of  structural  and  material  behaviour  has  gained 
momentum.  An  example  of  two  such  fields  of  research  that  have  attracted 
international  interest  in  the  past  decade  is  the  torsion  and  shear 
strength  of  reinforced  and  prestressed  concrete  members. 

The  complexity  of  box  girder  behaviour  that  has  frustrated 
designers  of  the  past  has  gradually  been  diffused  since  the  advent  of 
computer  technology.  The  evolution  of  computer-orientated  analytical 
design  methods  has  developed  rapidly  to  the  current  level  of  sophistica¬ 
tion  where  the  capabilities  of  structural  analysis  are  bounded  by  few 
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restrictions.  The  elastic  analysis  of  concrete  box  girder  structures 
is  now  well  refined,  and  research  effort  is  presently  being  focused  on 
the  non-linear,  post— cracking  region  of  behaviour. 

To  this  point  in  time,  non-linear  analytical  methods  of  analysis 
have  been  developed  to  trace  the  complete  load  deformation  path  from  the 
elastic  uncracked  state  through  to  failure.  However,  the  constituent 
behaviour  of  concrete  in  its  cracked  form  has  only  come  under  close 
scrunity  in  recent  years,  and  a  reasonably  comprehensive  understanding 
of  its  behaviour  is  only  beginning  to  emerge.  Thus,  the  logical 
progression  is  the  development  of  an  improved  analytical  computer  model 
that  incorporates  and  reflects  the  more  accurate  representation  of 
material  behaviour. 

1. 2  Thesis  Objective  and  Scope 

The  principal  objective  of  this  thesis  is  the  development  of  an 
analytical  computer  model  that  can  analyze  a  prestressed  concrete  box 
girder  of  arbitrary  cross-section  for  any  loading  combination  of  bending 
moment,  torque,  and  shear.  In  addition  to  estimating  cracking  and 
ultimate  loads,  the  complete  stress-deformation  description  of  the  box 
girder  is  provided  at  any  load  level  in  the  pre-cracked  or  post-cracked 
condition. 


To  assist  in  the  assessment  of  the  analytical  model  performance, 
seven  prestressed  double-cell  concrete  box  girders  were  cast  and  tested 
under  varying  load  combinations  of  torque,  bending  moment,  and  shear. 
Subsequently,  each  beam  was  analyzed  by  the  computer  model  for  the  same 
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respective  loading  conditions.  The  accuracy  of  the  analytical  modelling 
was  evaluated  upon  comparison  of  the  computer  model  results  with  the 
corresponding  experimental  test  results  and  available  theoretical 


predictions. 


CHAPTER  2 


RESEARCH  EVOLUTION 


2.1  Introduction 

This  section  is  not  an  exhaustive  state  of  the  art 
presentation,  but  simply  an  overview  of  research  developments 
demonstrating  the  evolutionary  progression  that  has  occurred  in 
experimentation  and  analysis  in  fields  of  study  pertaining  to  this 
thesis  topic.  The  five  associated  fields  of  research  encompassed  are 
the  three  loading  types  of  bending  moment,  torque,  and  shear  force, 
together  with  the  nature  of  prestressed  concrete  and  the  characteristics 
of  box  girder  behaviour. 

2 . 2  Review  of  Developments  to  Current  State  of  Art 

2.2.1  Experimental  Approach 

In  the  infancy  of  examination  of  reinforced  concrete  behaviour, 
the  complexity  involved  in  formulating  rational  analytical  solutions 
prompted  an  empirical  approach.  Detailed  experimental  programs  were 
undertaken  to  develop  an  understanding  of  behavioural  characteristics, 
and  the  recorded  test  data  was  used  to  develop  and  test  methods  of 
analysis  and  design. 

The  elastic  performance  and  ultimate  strength  of  reinforced 
concrete  in  bending  have  been  investigated  upon  numerous  occasions  in 
the  testing  of  scale  models  of  geometrically  complex  structures.  In 
certain  instances,  actual  structures  have  been  loaded  to  determine 
service  load  response.  An  example  of  a  bending  test  to  destruction  of 
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a  complex  scale  model  is  the  experimental  study  of  Bouwkamp,  Scordelis, 
and  Wasti1,  in  which  a  replica  of  a  typical  two-lane  reinforced 
concrete  bridge  was  cast  and  tested  to  failure. 

An  excellent  example  of  a  field  of  research  where  extensive 

experimental  testing  has  been  undertaken  because  of  the  difficulty  in 

developing  a  rational  analytical  model,  is  the  study  of  the  shear 

strength  of  reinforced  concrete  members.  Widespread  interest  in  this 

topic  has  been  sustained  as  the  "shear  failure"  mechanism  causes  a 

reduction  in  the  strength  of  structural  elements  below  flexural  capacity 

and  a  diminution  of  element  ductility.  To  this  point  in  time,  the 

phenomenon  is  still  not  clearly  understood,  and  consequently  code 

provisions  have  been  structured  to  prevent  such  a  failure.  In  the 

past  ten  years,  investigators  have  turned  their  attention  away  from  the 

shear  failure  mechanism  as  an  entity  in  itself,  and  have  directed  their 

research  effort  toward  closer  examination  of  the  contributing  shear 

components  present  across  a  concrete  crack;  aggregate  interlock  and 

dowel  action.  The  aggregate  interlock  effect  *  *  has  been  researched 

4, 5, 6, 7 

extensively  as  has  dowel  action  ,  and  a  better  understanding  of 

the  physical  mechanics  of  shear  transfer  has  now  emerged.  However,  the 
scope  of  experimentation  in  examination  of  the  latter  two  phenomena  has 
been  restricted  in  several  respects,  and  the  numerous  publications  of 
experimental  results  and  interpretation  are  not  in  complete  agreement. 

A  relatively  recent  state  of  the  art  paper  on  the  shear  strength  of 
reinforced  concrete  members  is  that  of  the  joint  ASCE-ACI  Task  Committee 
4268,  wherein  all  aspects  of  shear  transfer  are  reviewed.  Formulae  in 
the  latter  paper,  as  well  as  those  in  the  shear  provisions  of  most 
building  codes  for  reinforced  concrete,  have  a  purely  empirical  basis. 
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thus  demonstrating  the  valuable  contribution  of  experimental  research 
in  complex  structural  fields. 

Closely  associated  with  the  study  of  shear  is  that  of  torsion 
of  reinforced  concrete  members.  Only  in  the  last  ten  years  has  this 
avenue  of  research  attracted  considerable  attention.  Indeed,  the 
dramatic  explosion  in  technological  advance  in  this  area  is  evident 
upon  comparison  of  the  torsion  provisions  of  the  1963  ACI  Code  to  those 
of  the  1971  ACI  Code  and  CSA  Standard  A23.3.  In  a  similar  procedure  to 
that  adopted  for  shear,  the  development  of  torsion  code  provisions  of 
the  1971  ACI  Code  was  partially  dependent  upon  experimental  results9 »L 0 » 1 1 . 
The  corresponding  torsion  clauses  in  the  more  recently  published  CSA 
Standard  A23.3  are  quite  dissimilar,  a  different  theoretical  model  having 
been  adopted  as  a  basis  for  torsion  derivations.  The  respective 
theories  for  the  two  codes  will  be  addressed  in  the  following  section. 

From  an  experimental  aspect,  test  results  that  were  used  to  substantiate 
the  two  theories  were  in  reasonable  correspondence  but  were  interpreted 
differently. 

Of  paramount  importance  to  the  designer  is  the  nature  of 
interaction  of  all  three  loading  types  bending,  torsion,  and  shear. 

As  in  the  case  of  shear  and  torsion,  comprehensive  experimental 
programs11  were  undertaken  to  establish  interaction  curves  for  torsion 
and  bending, and  torsion  and  shear.  The  only  interaction  formulae  for 
torsion  and  bending  that  were  developed  on  a  purely  theoretical  basis 
are  those  of  Lampert12,  and  the  ensuing  discussion1 3  of  the  proponents 
of  the  empirical  and  theoretical  approach  highlights  the  current  state 
of  the  art  of  this  topic. 
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Experimental  programs1 4 * 1 5  * 1 6  devoted  to  prestressed  concrete 
research  have  been  conducted  in  the  same  manner  as  that  adopted  for 
reinforced  concrete.  Essentially,  the  performance  of  the  two  reinforced 
concrete  types  is  quite  similar  in  that  beams  of  concentrically 
prestressed  and  symmetrically  reinforced  concrete  are  analogous  in 
behaviour,  as  are  eccentrically  prestressed  and  unsymmetrically  rein¬ 
forced  concrete  beams.  Code  provisions  for  torsion  and  shear  of  pre¬ 
stressed  concrete  members  are  empirically  derived,  and  bear  a  close 
resemblance  to  the  corresponding  reinforced  concrete  clauses. 

Behaviour  of  box  girder  members  has  come  under  closer  scrutiny 
in  recent  years  as  a  large  percentage  of  concrete  bridges  are  of  the 
multi-cell  box  girder  type.  The  box  section  has  been  favoured  by  many 
designers  because  of  its  aesthetic  appearance,  efficiency  of  cross- 
section,  and  high  torsional  rigidity.  Since  the  experimental  research 
of  solid  and  single-celled  reinforced  concrete  members  is  applicable  in 
most  facets  of  behaviour  to  box  girder  structures,  research  of  multi-celled, 
hollow  members  has  primarily  been  directed  toward  those  behavioural 
aspects17  peculiar  to  concrete  box  girder  bridges.  Such  aspects  include 
cross-section  distorsion,  warping,  shear  lag18,  and  diaphragm  action. 
Invariably,  these  aspects  have  not  been  isolated  and  studied  individually, 
but  have  simply  been  incorporated  collectively  in  the  testing  of  scale 
models  of  concrete  bridges1’19. 

2.2.2  Development  of  Theoretical  Analyses 

Up  to  the  onset  of  cracking,  stresses  and  deformations  in  a 
determinate  prestressed  concrete  box  girder  of  simple  cross-section  are 
readily  evaluated  for  any  loading  combination.  The  assumptions  that 
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the  uncracked  concrete  behaves  as  a  homogeneous,  isotropic,  elastic 
material,  and  the  contribution  of  the  reinforcement  to  the  girder 
stiffness  is  small,  do  not  constitute  a  severe  approximation  of  the 
actual  structural  behaviour.  However,  if  the  cross-section  is  more 
complex  (multi-celled) ,  longitudinal  warping  restraint  is  present 
(boundary  conditions  or  diaphragm  action) ,  or  the  structure  is  indeter¬ 
minate,  the  simplified  approach  using  classical  elastic  theory  is  not 
accurate,  or  may  indeed  be  non-applicable.  For  such  a  situation,  a  more 
sophisticated  solution  procedure  is  required,  often  in  the  form  of  a 
computer-orientated  method  of  analysis. 

Upon  the  onset  of  cracking,  the  reinforcement  assumes  a  vital 
role  in  carrying  tensile  stresses  and  achieving  stress  redistribution. 

Of  the  three  loading  types,  bending,  torsion,  and  shear,  only  bending 
action  can  presently  be  analyzed  theoretically  to  yield  complete  defor¬ 
mation  behaviour.  At  the  current  state  of  knowledge,  shear  still  defies 
rational  theoretical  analysis,  whereas  several  theoretical  postulates 
have  been  presented  to  predict  torsional  behaviour. 

The  first  recognized  theoretical  model  to  predict  the  torsional 

r\  A 

strength  of  reinforced  concrete  was  that  of  Rausch  ,  whose  model 
consisted  of  a  network  of  bars  to  represent  the  action  of  reinforced 
concrete:  compression  concrete  bars  and  tension  bars  for  reinforcement. 
All  reinforcement  was  assumed  to  yield  simultaneously  at  failure.  In 
time,  this  model  was  modified  by  Cowan21  and  Anderson22  in  assuming  that 
the  concrete  carried  a  torque  at  failure  equal  to  the  cracking  strength 
of  an  unreinforced  beam.  Subsequently,  this  approach  was  questioned  by 
Hsu23,  who  maintained  that  real  behaviour  lay  between  the  two  extremes 
of  that  proposed  by  Cowan  and  Anderson  and  1958  German  Code  approach 


that  assumed  zero  cracked  concrete  torsional  strength.  The  basis  of 
Hsu's  research  was  the  development  of  a  modified  skew  bending  model. 

The  first  skew  bending  model  was  proposed  by  Lessig24,  and  is 
illustrated  in  Fig.  2.1.  The  depicted  failure  surface  is  comprised  of 
an  inclined  compression  zone  and  a  warped  failure  plane  delineated  by 
a  diagonally  inclined  crack  connecting  compression  zone  ends,  the 
inclination  of  the  compression  zone  being  dependent  on  cross-sectional 
geometry  and  ratio  of  longitudinal  to  transverse  steel  areas.  In 
proposing  his  modified  skew  bending  model,  Hsu  drew  attention  to  the 
fact  that  the  Lessig  theory,  being  an  upper  bound  solution,  considerably 
overestimated  ultimate  strength.  Hsu’s  principal  modifications  were 
adjusting  the  shorter  sides  of  the  failure  surface  to  intersect  the 
flange-web  interface  at  90°,  and  neglecting  the  torsional  contribution 
of  the  shorter  stirrup  legs.  Criticism50  of  Hsu's  skew  bending  theory 
stemmed  from  three  reservations: 

1.  Several  constants  were  difficult  to  evaluate. 

2.  Theory  had  to  be  adjusted  for  square  beams. 

3.  Dowel  action  of  longitudinal  bars  was  indeterminable. 

In  recognizing  that  the  skew  bending  theory  was  an  upper 
bound  solution  that  more  accurately  predicted  ultimate  bending  strength 
rather  than  ultimate  torsional  strength,  several  researchers  directed 
their  efforts  toward  refining  the  space  truss  theory  postulated  by  Rausch. 
To  this  point  in  time,  the  most  widely  accepted  space  truss  theory  is 
that  developed  by  Lampert  ,  the  theory  presented  in  a  form  directly 
applicable  to  design  by  Collins  and  Lampert  .  The  space  truss  model 
adopted  by  Lampert  is  shown  in  Fig.  2.2,  consisting  of  intermediate  shear 
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walls  and  longitudinal  reinforcement  considered  to  be  concentrated  into 
stringers  at  the  hoop  reinforcement  corners.  In  the  shear  walls,  the 
stirrups  act  as  posts  and  the  concrete  between  the  inclined  cracks 
provides  the  compression  diagonals.  The  angle  of  the  diagonals  with 
respect  to  the  beam  axis  is  assumed  to  be  constant  for  each  side.  In 
the  walls  that  govern  failure,  the  angle  is  such  that  both  the  longitu¬ 
dinal  and  stirrup  reinforcement  reach  their  respective  yield  points 
simultaneously.  For  this  reason,  the  model  is  a  space  truss  of  variable 
diagonal  inclination.  The  most  attractive  feature  of  Lampert's  space 
truss  theory  is  the  simplicity  of  the  design  equations  that  now  form 
the  basis  of  the  torsion  sections  in  the  Canadian  Code,  CSA  Standard 
A23.3.  However,  the  theory  is  not  completely  general  as 

1.  cross-section  must  be  underreinforced, 

2.  model  cannot  accommodate  shear, 

3.  only  St.  Venant  torsion  can  be  modelled  (no  support  or 

load  warping  restraint) 

Criticism  of  the  Collins  and  Lampert  theory  is  summarized  comprehensively 
in  the  discussion13  of  the  authors'  paper26  wherein  they  stated  that 
the  problem  of  torsion  and  bending  was  basically  solved. 

Thus  far,  only  the  strength  predictions  of  theoretical 
analyses  have  been  reviewed.  The  aspects  of  post-cracking  stiffness 
and  inelastic  member  deformation  have  not  been  addressed.  The  importance 
of  these  two  latter  aspects  is  paramount  since  the  calculation  of 
equilibrium  torsional  moments  in  a  statically  indeterminate  structure 
requires  not  only  statics  but  also  compatibility  conditions.  As  a 
point25  of  clarification,  an  equilibrium  torque  is  one  required  to 
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maintain  equilibrium  in  a  structure,  in  contrast  to  a  compatibility 
torque  that  maintains  structural  compatibility.  Both  Hsu27  and  Collins 

f\  c 

and  Lampert  have  proposed  post-cracking  stiffness  formulae  that  are 
derived  on  the  basis  of  their  respective  theories.  Two  important 
qualifications  of  Hsu’s  theoretical  approach  are  that  the  analysis 
hinges  on  the  empirical  derivation  of  an  equivalent  wall  thickness  for 
solid  and  "thick  walled"  cross-sections,  and  applicability  of  the  theory 
to  non-rectangular  sections  is  not  verified. 

The  theories  of  both  Hsu  and  Lampert  are  restricted  to  the 
pure  torsion  condition. 

2.2.3  Analytical  Computer  Solutions 

In  the  study  of  complex  modern  structures  such  as  box  girder 
bridges,  dams,  and  prestressed  concrete  nuclear  containment  vessels,  the 
experimental  approach  adopted  in  the  past  is  becoming  an  increasingly 
expensive  method  of  investigation.  Thus,  the  empirical  approach  is 
gradually  being  superseded  by  refined  analytical  computer-orientated 
methods.  Not  only  does  the  analytical  method  offer  a  considerable  saving 
in  terms  of  time,  but  the  progressive  improvement  in  computer  technology 
permits  the  development  of  computer  programs  of  increasing  complexity. 

In  the  context  of  the  analysis  of  box  girder  bridges,  the 
numerous  analytical  methods  and  models  that  have  been  developed  fall  into 
two  categories;  those  that  model  behaviour  prior  to  cracking,  and  those 
that  model  the  complete  structural  behaviour  to  failure. 

Of  those  methods  that  have  been  developed  for  the  uncracked 
section,  approximate  methods  of  analysis  based  on  simplified  structural 
behaviour,  such  as  the  equivalent  beam  grillage  or  anisotropic  slab  methods. 
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have  been  used  to  model  structural  systems.  Of  a  more  precise  nature, 
elaborate  methods  based  on  folded  plate  theory  have  been  developed  by 
several  researchers,  the  most  prominent  sustained  research  program 
having  been  conducted  by  Scordelis  et  alia28’29  at  the  University  of 
California,  Berkeley.  In  the  adaptation  of  the  folded  plate  theory  to 
the  analysis  of  box  girder  bridges,  each  component  plate  of  the  box 
girder  is  considered  as  an  assemblage  of  individual  elements,  the  bending 
of  each  plate  normal  to  its  plane  being  analyzed  by  plate  flexure  theory, 
and  the  in-plane  bending  analyzed  by  plane  stress  theory.  These  classical 
theories  necessitate  the  representation  of  the  applied  loading  by  a 
Fourier  Series,  with  the  result  that  computational  effort  is  considerable 
though  less  than  that  of  a  finite  element  solution.  In  an  effort  to 
reduce  the  number  of  equations  and  thus  reduce  computing  time  and 
programming  effort,  the  ’finite  strip  method’  has  been  proposed  by 
Cheung30.  In  this  method,  the  behaviour  of  each  plate  is  approximated 
by  an  assemblage  of  longitudinal  finite  strips  for  which  selected 
displacement  patterns  are  assumed  to  represent  the  behaviour  of  the  strip 
in  the  total  structure.  An  additional  advantage  of  this  simplified 
method  is  that  more  complex  material  properties  such  as  concrete  anisotropy 
can  be  readily  introduced.  However,  both  the  folded  plate  and  finite 
strip  methods  are  very  much  more  restricted  in  their  range  of  application 
than  the  finite  element  method  as  the  two  methods  can  only  be  applied  to 
box  girders  of  constant  cross-sectional  geometry.  The  greater  flexibility 
is  achieved  at  the  expense  of  computational  effort.  In  finite  element 
studies,  it  should  be  recognized  that  the  accuracy  achieved  is  dependent 
upon  assumptions  of  material  properties  and  fineness  of  the  structural 
mesh  subdivisions.  Generally,  results  closely  satisfy  compatibility, 
but  not  necessarily  equilibrium  in  the  continuum  unless  a  sufficiently 
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fine  mesh  is  used.  An  array  of  computer  programs  has  been  developed 
at  the  University  of  California,  Berkeley,  to  analyze  box  girder 
bridges  of  arbitrary  plan  and  general  cross-section. 

Beyond  cracking,  material  behaviour  and  structural  interaction 
are  complex.  Consequently,  analysis  is  almost  exclusively  performed  by 
the  finite  element  method.  Accurate  analytical  determination  of  the 
stress-deformation  condition  of  a  reinforced  concrete  structure  at  a 
certain  stage  of  cracking  is  complicated  by 

1.  Structural  response  is  governed  by  the  interaction  of  two 

component  materials  -  steel  and  concrete. 

2.  Stress-strain  relationships  for  concrete  and  reinforcement 

are  non-linear,  with  the  concrete  exhibiting  anisotropy. 

3.  Shear  transfer  phenomena  of  aggregate  interlock  and  dowel 

action  together  with  bond  slip  must  be  considered. 

4.  Failure  criteria  for  concrete  under  biaxial  stress  conditions, 

dowel  action,  and  bond  slip  have  to  be  established. 

5.  Because  of  non-linear  material  properties,  equilibrium  is  not 

easily  maintained,  and  considerable  iteration  must  be  performed 

within  load  increments. 

6.  Transverse  rigidity  and  warping  resistance  of  actual  diaphragms, 

as  well  as  the  intrinsic  transverse  rigidity  of  the  box  section, 

must  be  included. 

7.  Dramatic  load  reversal  can  occur  upon  application  of  the  first  load 

increment  after  prestress  transfer. 

The  initial  development  of  the  application  of  the  finite 
element  method  to  the  analysis  of  reinforced  concrete  beam  behaviour 
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was  undertaken  at  the  University  of  California.  In  the  earliest 
publication  of  the  application  of  the  finite  element  method,  Ngo  and 
Scordelis  analyzed  simple  beams  in  which  the  concrete  and  steel  rein¬ 
forcement  were  represented  by  two-dimensional  triangular  finite  elements, 
and  special  bond  linkage  elements  were  used  to  connect  reinforcement 
to  concrete.  Linear  elastic  analyses  were  performed  on  the  beams  with 
predefined  crack  patterns  to  determine  principal  stresses  in  the  concrete, 
and  stress  levels  in  the  reinforcement  and  bond  linkages.  Since  the 
initial  publication  of  Ngo  and  Scordelis,  the  flexibility  and  capability 
of  the  finite  element  method  in  this  particular  application  has  advanced 
dramatically  to  the  point  where  most  of  the  seven  stipulated  compli¬ 
cations  have  now  been  successfully  incorporated.  Such  an  analytical 
computer  model  is  that  developed  by  Trikha  and  Edwards34. 

2.3  Definition  of  Thesis  Approach 

2.3.1  Analytical  Model  Aspect 

In  recognition  of  the  complexity  encountered  in  the  analysis 
of  prestressed  concrete  box  girders  in  the  uncracked  and  cracked  states, 
a  finite  element  analytical  approach  has  been  adopted  in  determination 
of  both  strength  and  deformation  characteristics.  In  addition  to 
evaluating  strength-deformation  values  at  both  the  cracking  and  ultimate 
loads,  the  analytical  model  yields  a  comprehensive  description  of  the 
structural  behaviour  for  the  complete  load  range.  Generality  of  loading 
is  assured  in  that  the  external  loads  can  reflect  any  ratio  of  bending 
moment  to  torque  to  shear,  and  the  nature  of  the  loading  pattern  can 
be  altered  at  any  point  in  the  loading  sequence. 
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Analytical  flexibility  of  the  proposed  model  in  a  structural 
mechanics  context  is  exemplified  by  its  ability  to  simulate: 

1.  Interaction  between  the  two  material  components  -  steel 

and  concrete. 

2.  Any  reinforcement  pattern  of  conventional  steel  and  prestress 

strand . 

3.  Non-linear  material  behaviour. 

4.  Presence  of  diaphragms  and  the  cross-sectional  shear  rigidity 

of  the  box  girder  section. 

5.  Shear  transfer  across  concrete  cracks. 

In  essence,  the  strength  of  the  computer  model  approach  is 
its  comparative  freedom  from  analytical  and  structural  constraint. 

2.3.2  Experimental  Aspect 

To  evaluate  the  analytical  model  performance,  a  limited 
experimental  program  was  undertaken  in  which  the  test  specimens  were 
subjected  to  varying  combinations  of  bending  moment,  torque,  and  shear. 
To  be  more  representative  of  concrete  box  girders  employed  in  practice, 
all  test  specimen  cross-sections  were  multi-celled,  and  of  either 
rectangular  or  trapezoidal  shape.  The  higher  degree  of  complexity 
of  test  specimen  cross-sectional  geometry  also  provided  a  more  rigorous 
basis  of  comparison  for  the  computer  model. 
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FIG.  2.1  SKEW  BENDING  MODEL 


FIG.  2.2  SPACE  TRUSS  MODEL 


CHAPTER  3 


ANALYTICAL  COMPUTER  MODEL 


3.1  Introduction 

As  the  complexities  of  structural  behaviour  have  come  under 
closer  scrutiny,  more  exacting  analytical  methods  have  been  developed 
in  an  effort  to  overcome  past  analytical  failings  or  restrictions. 
Research  in  structural  engineering  has  often  turned  to  the  analytical 
model  as  the  most  appropriate  method  of  investigating  complex  structural 
behaviour*  Such  is  the  approach  in  this  thesis  to  the  study  of  the 
action  of  prestressed  concrete  box  girders  under  the  combined  loading 
of  torque,  bending  moment,  and  shear. 

Of  the  analytical  methods  that  are  most  commonly  used  in  the 
analysis  of  concrete  box  girders,  only  four  are  capable  of  representing 
the  most  important  characteristics  of  box  girder  behaviour,  those 
characteristics  being  longitudinal  bending,  St.  Venant  torsion,  trans¬ 
verse  distortion,  longitudinal  warping,  and  shear  lag.  Folded  plate 
theory,  finite  strip  theory,  finite  element  theory,  and  shell  theory  all 
yield  satisfactory  results  in  an  elastic  analysis,  but  only  the  finite 
element  method  is  capable  of  simulating  the  inelastic  behaviour  of 
concrete  in  its  uncracked  and  cracked  states.  Primarily  for  this  reason 
a  finite  element  approach  has  been  adopted.  A  simplified  finite  element 
mesh  similar  to  those  employed  in  this  analysis  is  exhibited  in  Fig.  3.1 

In  the  development  of  the  computer  model,  the  following 
capabilities  have  been  incorporated: 
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1.  Representation  of  any  general,  thin-walled  cross-section 

comprised  of  linear  segments. 

2.  Representation  of  any  in-plane  loading  system  (prestressing 

included) . 

3.  Complete  description  of  beam  behaviour  in  the  elastic  and 

inelastic  regions  up  to  ultimate  failure. 

4.  Accurate  portrayal  of  material  response  to  the  action  of 

torsion,  bending,  and  shear  in  both  uncracked  and  cracked 

concrete  regimes. 

5.  Ability  to  model  both  diaphragm  action  and  the  transverse 

rigidity  of  a  box  girder  without  diaphragms. 

The  coding  of  the  analytical  computer  model  whose  main  program 
and  subroutines  are  listed  in  Appendices  B  and  C  respectively,  is 
written  in  the  Fortran  IV  language  compatible  with  IBM  System/360  and 
System/370.  All  computer  analyses  were  achieved  through  the  use  of  the 
Amdahl  470  computer  at  the  University  of  Alberta  Computing  Centre,  and 
the  several  system  subroutines  utilized  by  the  program  are  public  library 
routines. 

3 . 2  Finite  Element  Types  Employed 

3.2.1  Concrete  Wall  Element 

The  principal  design  characteristic  of  concrete  box  girders  is 
the  structural  efficiency  achieved  by  the  absence  of  the  central  concrete 
core,  thus  significantly  reducing  the  member's  dead  weight  with  a  small 
loss  of  strength.  As  a  consequence  of  the  void,  the  thickness  of 
concrete  walls  that  define  the  girder's  cross-sectional  geometry  is 


19. 


usually  moderately  small  compared  with  the  member's  depth  and  width. 

When  the  wall  thickness  is  comparatively  small  such  that  it  can  be 
defined  as  being  "thin",  the  box  girder  cross-section  can  be  envisaged 
as  an  assemblage  of  flat  plates,  or  in  the  analytical  model  as  a  mesh  of 
appropriate  plane  stress  finite  elements. 


Since  it  is  most  desirable  for  the  sake  of  economy  and  efficiency 


that  the  concrete  box  girder  walls  be  modelled  by  plane  stress  elements, 
the  classification  as  to  whether  a  wall  thickness  is  "thick"  or  "thin" 
is  a  prime  concern.  Thus,  it  is  most  appropriate  that  this  issue  be 
examined  in  detail.  In  the  following  derivation,  a  single  cell,  uniformly 
thick  box  girder  will  be  considered,  as  exhibited  in  Fig.  3.2(a). 


For  a  thin-walled  single  cell,  the  uniform  shear  flow  is  given 


by 


T 


(3.1) 


q  2A' 


in  which  T  =  applied  torque,  and  A'  =  area  enclosed  by  corner  longitu¬ 
dinal  bars. 


Consequently,  the  uniform  wall  shear  stress  is  defined  by 


T 


(3.2) 


V 


2A't 


in  which  t  =  wall  thickness 


K 


s 


(3.3) 
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in  which  K  ^  =  torsion  constant  for  the  equivalent  open  cell,  and  K  ^  = 
torsion  constant  proportional  to  the  additional  torque  arising  from  closing 
an  open  cell. 


Referring  to  Eq.  3.3,  the  St.  Venant  torsion  constant  for  the 
uniform  shear  flow  of  the  closed  cell  is  expressed  by 


K 


,  (2A')Z 

,  ds 

i  T 


(3.4) 


in  which  s  =  perimeter  length  coordinate  for  the  closed  cell  cross-section, 
and 


K 


si 


-  kl 


(3.5) 


in  which  b  =  cross-sectional  segment  length  of  thickness  t. 


Therefore,  for  the  closed  single  cell  shown  in  Fig.  3.2(a), 
the  St.  Venant  torsion  resultant  for  the  equivalent  open  cell  is  given  by 


si 


K  .  T 
si 

K 


(3.6a) 


Substituting  for  Kg^  and  Kg, 


j  Et3  ds  j,  ~  T 
Tsl  =  ( 2A^ 

t2  (b+d)2  T 

3b2d2 

in  which  b  and  d  are  the  cross-sectional  dimensions  as  shown  in  Fig. 


(3.6b) 

3.2(a). 


With  reference  to  Fig.  3.2(b),  the  torque  per  unit  length 
developed  by  the  equivalent  open  cell  shear  stress  distribution  is 


given  by 
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Ay  t2 

6 


(3.8) 


in  which  Ay  =  variation  in  wall  shear  stress  from  mid-thickness  to 
surface. 


For  the  complete  cell. 


1 

6 


<J>  Ayt2  ds 


=  j  Ayt2  (b+d) 


(3.9) 


Equating  the  right  hand  sides  of  Eq.  3.7  and  Eq.  3.9  yields 


Ay 


(b+d)  T 

"b2d2" 


From  Eq.  3.2, 


y 


JT _ 

2bd  t 


(3.10) 


(3.11) 


Therefore,  the  ratio  of  the  variation  in  the  wall  shear  stress 
distribution  from  the  mid-thickness  to  the  surface,  to  the  uniform  shear 
stress  is  given  by 


Ay  =  2  (b+d)  t 

y  bd 

A 

-  TT  (3.12) 


where  A  =  cross-sectional  area  of  cell, 
c 

Thus,  the  definition  of  a  "thin-walled'*  box  girder  hinges  on 
the  degree  of  variation  in  the  shear  stress  distribution  that  is 
acceptable.  Many  authorities  postulate  a  threshold  value  of  10%  that 
distinguishes  a  "thick"  from  a  "thin"  box  girder  wall. 
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In  employing  plane  stress  rather  than  three-dimensional  finite 
elements,  the  total  number  of  degrees  of  freedom,  and  thus  equations  of 
equilibrium,  is  reduced  considerably.  As  the  computer  model  will  be 
used  to  analyze  only  uniform  beams  of  rectangular  and  trapezoidal  cross- 
sections  in  this  particular  study,  a  rectangular  finite  element  was 
chosen  as  the  basic  concrete  wall  element.  The  higher  order  rectangular 
finite  element  has  three  degrees  of  freedom  per  node,  namely  two  in-plane 
translational  degrees  of  freedom  and  a  rotational  degree  of  freedom 
orthogonal  to  the  element  plane. 

To  reduce  computing  costs,  the  behaviour  of  each  rectangular 
concrete  finite  element  is  characterized  by  the  stress-strain  condition 
at  the  centroid  of  the  element.  Although  this  simplification  approximates 
real  behaviour,  the  degree  of  the  approximation  will  be  demonstrated  to 
be  of  minor  significance. 

A  detailed  description  of  the  rectangular  finite  element  chosen 
for  this  analytical  model  is  given  in  Section  3.3. 

3.2.2  Reinforcement  Element 

The  three  distinct  types  of  reinforcement  represented  in  the 
computer  model  are: 

1.  Prestress  reinforcement. 

2.  Conventional  bar  reinforcement. 

3.  Bond  spring  linkage. 

All  three  are  modelled  by  a  one-dimensional  constant  strain 
finite  element  whose  equations  of  equilibrium  are  given  below: 
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R, 

K'  -K' 

r . 

i 

= 

i 

R. 

-K'  K’_ 

r . 

r 

v  J' 

(3.13) 


where  R^,  -  forces  at  nodes  i  and  j  respectively,  K'  =  stiffness 

term,  and  r^,  r^  =  displacements  at  nodes  i  and  j  respectively. 

The  nature  and  derivation  of  the  reinforcement  element 
stiffness  elements  K'  will  be  treated  in  detail  in  Section  3. 4. 2.1. 


3.2.3  Diaphragm  Elements 

Since  the  inception  of  box  girder  design,  the  need  was 
recognized  for  the  provision  of  transverse  plates  to  strengthen  the 
girders’  cross-sectional  rigidity,  distribute  shears  from  webs  to 
bearings  over  supports,  and  to  prevent  excessive  cross-sectional 
distortion.  In  practice,  all  box  girders  have  diaphragms  over  supports, 
while  the  number  of  intermediary  diaphragms  varies  considerably  depending 
upon  the  span  and  cross-sectional  shape. 

Since  the  prime  function  of  the  diaphragm  finite  element  in 
the  analytical  model  is  to  provide  cross-sectional  shear  rigidity,  a 
lower-order  element  is  quite  satisfactory.  To  accommodate  both  the 
rectangular  and  trapezoidal  cross-sectional  shapes,  the  diaphragm  element 
chosen  is  a  bi-linear  isoparametric  serendipity  element,  as  shown  in 
Fig.  3.3.  This  element  possesses  two  translational  degrees  of  freedom 
per  node,  eight  degrees  of  freedom  in  all  per  element,  with  the  displace¬ 
ment  vector  defined  by: 


u 

<<{)> 

0 

r 

U_ 

V 

k  4 

0 

<(j)> 

V 
k  J 

(3.14) 
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where  shape  function  (f^  =  ^  (1+n  r|) (1+^C)  (3.15) 

where  i  is  the  node  number  and  (£.,  r|.)  are  the  non-dimensional  nodal 

i  1 

coordinates . 

The  geometric  coordinates  (x,  y)  are  mapped  in  the  x-y  plane 
using  the  same  shape  functions  as  above. 

<x  y>  =  «f>>g  [x  y]  (3.16) 

where 

4>lg  -  f-  a+n.nHi+5^)  (3.17) 

Thus,  the  analytical  model  is  capable  of  representing  the  stiffness  for 
a  diaphragm  of  any  quadrilateral  geometry.  Derivation  of  the  stiffness 
of  a  bi-linear  isoparametric  serendipity  element  is  detailed  by 
Zienkiewicz3 5.  In  the  course  of  the  computer  program  logic,  the  stress- 
strain  condition  of  the  diaphragm  elements  is  not  monitored  as  all 
diaphragm  elements  are  assumed  to  remain  elastic.  A  listing  of  the 
coding  of  the  appropriate  subroutines  is  given  in  Appendix  D. 

In  accurately  representing  diaphragm  action,  three  distinct 
behavioural  aspects  are  addressed.  Firstly,  the  in-plane  action  of 
actual  diaphragms  is  incorporated  through  the  use  of  the  bilinear 
finite  element  with  a  full  concrete  constitutive  matrix.  Secondly,  the 
transverse  rigidity  of  a  box  girder  length  without  diaphragms  is 
represented  by  equivalent  diaphragm  elements  similar  to  the  above,  the 
distinguishing  feature  being  that  only  a  shear  stiffness  term  is  present 
in  the  constitutive  matrix.  This  aspect  is  treated  in  detail  in 
Section  3.4.2. 9.  The  third  and  last  aspect  of  diaphragm  action  introduced 
into  the  analytical  model  is  the  out-of-plane  warping  restraint  developed 
by  actual  diaphragms.  Within  the  computer  model,  those  serendipity 
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elements  representing  actual  diaphragm  segments  contribute  predetermined 
longitudinal  stiffness  terms  to  the  total  beam  stiffness  matrix  at  each 
of  the  element's  four  corner  nodes.  However,  derivation  of  the  warping 
restraint  at  the  corner  nodes  is  dependent  upon  the  classification  of 
the  diaphragm  thickness  as  being  "thick"  or  "thin".  The  two  respective 
warping  derivations  are  fully  treated  in  Sections  3. 4. 2. 7  and  3. 4. 2. 8. 

3.3  Choice  of  Plane  Stress  Concrete  Wall  Element 

Three  plane-stress  rectangular  finite  elements  were 
conspicuous  in  a  literature  search  for  an  appropriate  concrete  finite 
element  to  be  utilized  in  the  analytical  model;  the  element's 

i 

prerequisities  were  established  in  3.2.1.  The  three  elements  were: 

1.  McCleod  element36 

2.  Scordelis  element37 

3.  Sisodiya  and  Ghali  element38 

Of  the  three  plane-stress  elements  represented  in  the  Sisodiya  and 
Ghali  paper38,  parallelogram  element  PQC3  was  chosen  for  this 
comparative  investigation. 

The  method  of  comparison  adopted  in  choosing  the  most  suitable 
element  simply  involved  the  modelling  of  identical  single-cell  concrete 
girders  using  each  of  the  above  elements,  the  box  girders  being  subjected 
to  loading  patterns  of  torque,  bending  moment,  and  shear.  The  analysis 
was  strictly  elastic,  with  material  behaviour  simplified  by  omitting 
reinforcement  and  approximating  concrete  stiffness  as  being  linear  and 
isotropic.  To  assess  convergence  and  accuracy,  three  finite  element 
meshes  of  varying  degrees  of  refinement,  as  shown  in  Fig.  3.4,  were  used. 
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As  a  basis  of  comparison,  both  stress-strain  and  deformation 
predictions  were  considered.  The  stress  predictions  for  the  application 
of  a  simple  bending  moment  and  a  pure  St.  Venant  torque,  both  uniform 
over  the  central  length  of  the  beam,  are  plotted  in  Figs.  3.5(a)  and 
3.5(b)  respectively.  All  three  elements  perform  equally  well  in  both 
estimating  direct  stress  and  shear  stress  levels,  and  converging  to  the 
respective  theoretical  values.  The  vertical  deflection  of  the  central 
beam  cross-section  arising  from  the  application  of  bending  moment,  and 
the  differential  rotation  of  the  central  beam  length  subjected  to  a 
uniform  St.  Venant  torque,  are  given  in  Table  3.1  for  each  of  the  three 
element  types.  As  in  the  first  basis  of  comparison,  no  particular 
element  was  significantly  superior  to  another  in  its  accuracy. 

In  a  closer  scrutiny  of  the  results,  the  Sisodiya-Ghali  element 
was  observed  to  have  yielded  an  upper  bound  estimate  to  the  theoretical 
bending  deformation.  Since  deformation  characteristics  are  of  prime 
concern  in  the  principal  analytical  model,  this  element  was  considered 
the  least  preferable.  In  choosing  between  the  Scordelis  and  McCleod 
elements,  inter-element  displacement  compatibility  was  the  criterion 
selected  as  a  basis  of  comparison.  Whereas  the  Scordelis  element  does 
not  consistently  exhibit  complete  displacement  compatibility37,  the  choice 
of  the  nodal  rotational  degrees  of  freedom  and  the  displacement  functions 
for  the  McCleod  element  is  such  that  full  boundary  compatibility  is 
achieved.  Thus,  the  McCleod  element  was  chosen  as  the  marginally  prefer¬ 
able  plane  stress  finite  element  for  this  application. 

To  satisfy  the  requirement  of  complete  displacement  compatibility, 
two  McCleod  finite  element  types  are  employed.  For  each  element  type, 
the  rotational  degrees  of  freedom  are  defined  alternately  as  either 


. 
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3v  3u 

9  .  =  (t— ).  or  0  .  =  (-r — ).,  the  rotation  at  node  1  of  element  type  1 
zi  ax  i  zi  ay  i  JV 

3v 

being  defined  as  0^  =  (g^)^,  an<^  t^ie  rotation  at  node  1  of  element  type 

3u 

2  defined  as  0^  =  (-^p-)^.  By  allocating  element  types  such  that  adjacent 
elements  have  different  type  designations,  each  element  node  will  have 
a  unique  rotation.  Figures  3.6(a)  and  3.6(b)  illustrate  element  types  1 
and  2  respectively,  and  an  element  assemblage  for  a  two  dimensional 
beam  is  displayed  in  Fig.  3.6(c). 


The  McCleod  displacement  functions  chosen  to  satisfy  boundary 
compatibility  are 

u  =  +  A^x  +  A6y  +  A12xy  +  A^2  +  A1()xy2  (3.18) 

v  =  A,.  +  AgX  +  A2y  +  AgXy  +  A^x2  +  A^x2y  (3.19) 

in  which  A_^  =  coefficients  numbered  in  such  an  order  to  facilitate 
inversion  of  [A]  matrices. 


Maintenance  of  boundary  compatibility  is  demonstrated  as 
follows.  For  y  =  constant  =  k,  the  displacements  will  have  the  form 

u  =  (A4  +  kA6  +  k2Ax)  +  (A?  +  kA12  +  k2AlQ)x  (3.20) 

v  =  (A^-  +  kA2)  +  (Ag  +  kAg)x  +  (Ag  +  kA-^x2  (3.21) 

In  this  example,  u  is  a  linear  function  of  x  and  can  be 

defined  by  2  nodal  translations  in  the  x  direction,  and  v  is  a  quadratic 

function  of  x  and  can  be  defined  by  two  nodal  translations  in  the  y 

3v 

direction  together  with  the  single  edge  rotation  0  =  .  Similarly, 

nodal  displacements  are  uniquely  defined  along  x  =  constant  boundaries. 

Thus,  displacement  compatibility  is  maintained  across  element  boundaries 
and  throughout  the  structure. 
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The  stiffness  formulation  of  the  McCleod  finite  element 
chosen  differs  from  the  original  derivation  through  shifting  the  origin 
of  the  element's  local  axes  to  the  centroid,  and  changing  the  node 
numbering  sequence.  In  adopting  these  two  changes,  the  matrices 

[A] -1  for  element  types  1  and  2  (Table  l)36  have  had  to  be  reformulated. 

The  redefined  matrices  are  given  in  Table  3.2.  It  should  be  noted  that 
the  element  dimensions  have  been  redefined  as  (2A  x  2B)  as  illustrated 
in  Figs.  3.6(a)  and  3.6(b). 

The  plane  stress  constitutive  relationship  is  defined  as 

{a}  =  [D] {e}  (3.22) 

Thus,  the  element  stiffness  matrix  [k]  is  given  by  the 
relationship 

[k]  =  [A-1]T[k]  [A’1]  (3.23) 

in  which 

[k]  =  J  [B]T  [D]  [B]  dxdydz  (3.24) 

J  vol 

[B]  is  the  matrix  relating  element  strains  to  element  nodal  displacements. 

3.4  Computer  Program  Description 

3.4.1  General  Characteristics 

This  finite  element  computer  program  has  been  developed  as  an 
analytical  model  to  predict  the  behaviour  of  prestressed  or  reinforced 
concrete  box  girders  subjected  to  torsion,  bending,  and  shear.  Complete 
stress-strain  and  deformation  information  is  provided  at  any  specified 
load  level  in  both  the  uncracked  and  cracked  states  up  to  the  point  of 
girder  failure. 


' 
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Finite  elements  employed  in  the  program  consist  of  rectangular 
plane  stress  concrete  elements,  one-dimensional  reinforcement  elements, 
and  plane  stress  quadrilateral  diaphragm  elements.  Provision  is  included 
within  the  program  to  represent  the  presence  of  a  reinforcing  steel 
mesh  within  any  of  the  concrete  elements.  Those  reinforcement  elements 
whose  bond  with  the  adjoining  concrete  elements  is  suspect  to  deterioration 
during  loading,  are  connected  to  adjacent  concrete  nodes  through  bond 
spring  linkages  whose  stiffnesses  are  modified  as  loading  progresses. 

At  every  beam  cross-section  delineated  by  concrete  element  nodes, 
equivalent  diaphragm  elements  are  introduced  to  simulate  the  box  beam's 
transverse  cross-sectional  rigidity. 

Any  loading  pattern  can  be  superimposed  upon  the  analytical 
model  through  the  judicious  choice  of  nodal  force  combinations,  the 
only  restriction  being  that  force  directions  cannot  be  orthogonal  to 
plane  stress  concrete  elements.  The  presence  of  prestress  forces  is 
readily  accommodated. 

To  accurately  reproduce  concrete  behaviour,  the  program  is 
capable  of  modelling  non-linear  material  characteristics.  Naturally, 
non-linearity  is  not  restricted  solely  to  concrete,  as  the  prestress 
and  conventional  reinforcement  and  bond  spring  linkages  all  exhibit 
non-linear  behaviour.  The  loading  sequence  is  an  incremental  one 
consisting  of  the  superposition  of  successive  load  increments.  Following 
the  application  of  each  load  increment,  all  material  elements  are 
checked  for  deviation  from  their  pre-defined  behavioural  paths,  the 
probability  of  deviation  occurring  being  reduced  through  the  use  of  the 
Runge-Kutta  method.  If  significant  deviation  is  detected,  equilibrium 
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is  restored  using  the  modified  Newton-Rapson  iterative  method. 
Probability  of  material  deviation  within  a  subsequent  load  increment 
is  reduced  by  modifying  all  element  stiffnesses  at  the  end  of  the  current 
load  increment  using  a  tangent  stiffness  formulation.  A  detailed 
description  of  the  latter  method  of  stiffness  adjustment  and  the  Runge- 
Kutta  method  is  given  in  Section  3.4.4. 

To  solve  the  large  set  of  equilibrium  equations  efficiently, 
a  banded-block  solution  process  has  been  adopted.  At  any  instant  of 
the  solution  process,  only  two  stiffness  blocks  of  half-band  width 
are  stored  in  core.  Efficient  transmission  of  the  stiffness  blocks  in 
and  out  of  core  storage  is  achieved  through  the  use  of  public  library 
system  subroutines  provided  by  the  University  of  Alberta  Computing 
Services. 

The  non-linear  behaviour  of  a  rectangular  concrete  element 
under  bi-axial  stresses  is  characterized  by  the  stress-strain  condition 
at  the  centroid  of  the  element.  Once  the  principal  tensile  centroidal 
stress  exceeds  the  concrete  tensile  strength,  the  concrete  element  is 
designated  as  cracked,  the  crack  direction  thereafter  remaining  fixed. 

The  crack  inclination  is  derived  from  the  centroidal  stress  conditions 
at  cracking.  If  the  monitored  crack  width  should  close  immediately 
following  prestress  transfer,  the  element  is  once  again  designated  as  an 
uncracked  element.  In  reconstituting  the  stiffness  of  the  concrete 
element  following  cracking,  the  aggregate  interlock  and  dowel  effects 
are  taken  into  account  in  determining  the  shear  rigidity  across  the 
crack.  At  all  but  the  lowest  stress  levels,  concrete  behaves  as  an 


anisotropic  material. 
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At  the  conclusion  of  each  load  increment,  all  elements  are 
checked  to  detect  local  failure.  Structural  failure  occurs  when  one  of 
the  principal  prestress  strand  elements  fails  in  tension  or  a  concrete 
element  crushes.  In  reality,  the  modelled  structure  might  well  not  have 
failed  under  the  above  circumstances,  but  its  full  load  carrying  capacity 
will  have  been  attained  and  its  subsequent  highly  inelastic  behaviour 
will  indicate  imminent  collapse.  Should  failure  not  have  occurred  following 
the  application  of  the  load  increment,  a  comprehensive  summary  of  stress, 
strain,  and  deformation  values  will  be  printed  out. 

To  enable  large  program  runs  to  be  monitored  during  their 
execution,  several  duplicate  print  statements  were  inserted  in  the 
output  subroutine.  The  fully  comprehensive  output  for  each  load 
increment  is  assigned  to  a  high  speed  printing  device,  whereas  simulta¬ 
neously,  a  dramatically  smaller  representative  output  is  viewed  on  a 
terminal  screen  that  can  subsequently  assume  the  role  of  the  controlling 
device.  Termination  of  execution  of  the  program  can  be  prompted  once 
irregular  behaviour  is  observed. 

3.4.2  Simulation  of  Material  Behaviour 

3. 4. 2.1  Concrete  Stiffness  Under  Bi-axial  Stresses: 

Derivation  of  the  concrete  constitutive  matrix  under  bi-axial  stress 
conditions  closely  follows  the  research  of  Liu,  Nilson,  and  Slate39. 

The  form  of  the  stress-strain  equations  for  concrete  under  uniaxial  or 
biaxial  stress  conditions,  and  the  resulting  constitutive  matrix  proposed 
by  the  above  authors  are  easily  incorporated  in  an  incremental  iterative 
finite  element  analysis  as  developed  in  this  thesis. 
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The  stress-strain  relationship  for  concrete  under  biaxial 
stresses  is  given  by  the  equation 
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(1-va) [l+( 


(t-v  a) 


f  -  2><f> +  <f>2i 

s  p  p 


(3.25) 


in  which  a  =  stress  in  direction  considered,  e  =  strain  in  direction 
considered,  E  =  initial  modulus,  V  =  Poisson's  ratio,  a  =  ratio  of 
orthogonal  stress  to  stress  in  direction  considered,  Op  =  ultimate 
compressive  strength,  £p  =  strain  at  point  of  ultimate  strength,  and 


Liu,  Nilson,  and  Slate  state  that  the  equation  is  applicable 
to  concrete  in  biaxial  compression  only,  but  in  this  analysis,  the 
equation  is  used  in  biaxial  compression  and  compression-tension  conditions, 
where  in  the  latter  state  the  orthogonal  principal  stress  is  tensile. 


For  uniaxial  loading,  the  above  equation  simplifies  to: 


a  = 
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(3.26) 


This  theoretical  equation  was  compared  to  the  experimental 
results  of  a  concrete  cylinder  tested  during  the  experimental  program. 
The  two  curves  exhibited  in  Fig.  3.7  correspond  closely,  demonstrating 
the  reasonably  accurate  approximation  of  the  above  equation  to  real 
behaviour . 


The  two  distinct  load  conditions  encountered  in  the  analytical 
model  are  the  incremental  and  total  load  cases.  In  the  formulation  of 
the  concrete  constitutive  matrix,  the  stiffnesses  in  the  two  orthogonal 
principal  directions  are  derived  on  a  tangent  moduli  basis  for  incremental 
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loading,  and  a  secant  moduli  approach  during  the  iterative  process  when 
the  total  loading  condition  prevails.  Thus,  the  concrete  constitutive 
matrix  is  expressed  in  the  form: 
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and  where  in  the  total  loading  condition, 
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for  biaxial  compression  and  compression  tension.  However,  in  a  uniaxial 
compression  state 
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(3 . 27d) 


while  for  biaxial  tension,  tension  compression,  and  uniaxial  tension 
conditions, 

E»b  =  E  (3 . 27e) 


During  the  application  of  load  increments, 
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for  biaxial  compression  and  compression  tension.  For  the  uniaxial 
compression  state 
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while  for  biaxial  tension,  tension  compression,  and  uniaxial  tension 
conditions. 
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(3 . 27h) 


El,  is  defined  in  a  similar  manner  to  E ’  . 
2b  lb 


In  defining  concrete  moduli  under  biaxial  stress  conditions, 
the  nature  of  the  two  associated  orthogonal  principal  stresses  is  the 
criterion  used  in  establishing  the  biaxial  state.  However,  such  a 
definition  is  misleading  if  one  principal  compressive  stress  is  much 
larger  than  its  orthogonal  compressive  stress  such  that  the  strain  in 
the  orthogonal  direction  is  tensile.  Since  the  behaviour  of  concrete  is 
characterized  in  such  a  conflicting  situation  by  its  strain  state,  not 
stress  state,  a  condition  of  biaxial  compression  only  prevails  when 
the  strain  in  the  two  principal  strain  directions  is  compressive.  In 
its  application  to  the  analytical  model,  use  of  a  principal  strain  rather 
than  principal  stress  criterion  changes  the  concrete  moduli  little,  as 
the  tangent  stiffness  of  concrete  in  compression  at  low  stress  levels 
is  close  to  the  initial  tangent  modulus. 
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3. 4. 2. 2  Concrete  Crack  Width:  Janjua  and  Welch40  propose 
that  the  concrete  crack  width  be  given  by: 
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where  L  =  average  crack  spacing,  R  =  ratio  of  extreme  fibre  distance 
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from  neutral  axis  to  distance  of  steel  centroid  to  neutral  axis,  f  = 
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steel  stress,  and  E  =  steel  elastic  modulus. 
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(3.29a) 


in  which  t  =  concrete  cover  and  D  =  bar  diameter. 


All  the  above  variables  are  in  kip  and  inch  units. 

Unfortunately,  the  above  formulation  does  not  lend  itself 
readily  to  inclusion  within  the  analytical  model  logic  as  the  term  R  is 
not  easily  determined.  Moreover,  the  average  crack  spacing  expression 
shown  above  does  not  accurately  predict  the  test  beam  observations. 
Consequently,  elaborate  formulations  were  discarded  in  preference  for 
the  simple  expression: 


W  =  e  •  L' 
av .  c  av . 


(3.29b) 


in  which  -  finite  element  centroidal  direct  strain  perpendicular  to 
crack  direction  and  L*  =  observed  average  crack  spacing. 

3. 4. 2. 3  Aggregate  Interlock:  Once  plain  concrete  cracks, 
shear  is  still  able  to  be  transmitted  across  the  crack  through  interlock 
of  the  two  adjacent  rough  surfaces.  The  level  of  shear  that  can  be 
transferred,  however,  has  been  a  subject  of  constant  conjecture  and 
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research  in  the  development  of  analytical  models.  The  wide  diversity 
of  opinion  is  reflected  in  the  two  opposing  schools  of  thought  that 
support  the  Space  Truss  and  Skew  Bending  Analogies. 

Whereas  a  common  approach4 1  has  been  to  assume  that  a  constant 
percentage  of  the  concrete  shear  strength  is  retained  after  cracking, 
the  treatment  of  the  aggregate  interlock  effect  in  this  analytical  model 
is  based  on  the  research  conducted  by  Houde  and  Mirza4.  The  influence 
of  the  three  parameters  of  crack  width,  concrete  strength,  and  maximum 
aggregate  size  on  the  shear  rigidity  modulus  was  examined,  and  the 
results  lead  to  the  development  of  the  relationship: 
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where  A  =  shear  rigidity  modulus  and  c  =  crack  width. 


(3.30) 


All  of  the  above  variables  are  expressed  in  inch  pound  units.  As 
observed  in  the  above  expression,  the  effect  of  aggregate  size  was  found 
to  be  negligible. 


The  expression  above  was  derived  from  experimental  measurements 

made  in  the  range  of  crack  widths  of  2  x  10“ 3  to  20  x  10” 3  inches.  Crack 

widths  smaller  than  2  x  10” 3  were  difficult  to  accurately  control.  In 

the  uncertain  region  beyond  =  500,  the  authors  have  suggested  that  the 

curve  for  the  rigidity  modulus  is  asymptotic  to  the  line  A  =  G  (shear 

modulus  of  uncracked  concrete)  for  large  values  of  Such  a  supposition 

does  seem  severe,  however,  and  thus  the  equation  for  A  has  been  applied 

to  the  range  of  ^  >  500  in  the  absence  of  research  that  indicates 

otherwise.  The  curve  for  the  rigidity  modulus  intersects  the  line 

A  =  G  close  to  —  =  1000. 

c 
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3. 4. 2. 4  Dowel  Effect:  In  bridging  across  the  concrete 
crack,  reinforcement  not  only  restricts  the  widening  of  the  crack  such 
that  substantial  aggregate  interlock  can  develop,  but  also  offers  shear 
resistance  normal  to  its  axis.  This  shear  resistive  force  developed  in 
the  reinforcement  is  termed  the  "dowel  force" .  Since  the  dowel  effect 
is  only  significant  across  cracks  that  have  experienced  considerable 
shear  displacement,  reinforcement  must  be  highly  stressed  in  tension  to 
develop  dowel  action.  Principal  longitudinal  tension  reinforcement  in 
an  underreinf orced  beam  is  such  an  example. 


The  research  of  Houde  and  Mirza4  is  used  in  quantitatively 
defining  dowel  stiffness.  The  dowel  load-displacement  relationship  is 
given  by: 

Vd  =  2000.  Df.  Acrack  (3.31) 


where  -  dowel  force,  ^crac|c  =  shear  displacement  across  crack,  and 
Df  =  dowel  failure  force. 


D.  =  40  b  (f')1/3 

t  n  c 


(3.32) 


in  which  b  =  net  beam  width, 
n 

All  units  are  in  pounds  and  inches.  Embedment  length,  bar  size  or 
arrangement,  and  axial  stress  in  reinforcement  below  yield  do  not  have 
a  pronounced  effect  upon  dowel  action. 


Appraisal  of  the  effect  of  parameters  in  addition  to  those 
already  mentioned  is  given  in  other  publications.  The  effect  of 
inclination  of  reinforcement  to  crack  direction  on  dowel  strength  was 
investigated  by  Dulacska5,  and  although  an  expression  was  derived  for 
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the  dowel  failure  force,  a  simple  theoretical  relationship  could  not  be 
found  for  deformations.  Bauman's42  research  demonstrated  that  positioning 
of  stirrups  between  a  diagonal  crack  and  the  support  did  not  increase 
the  dowel  strength  if  the  distance  between  crack  and  stirrup  exceeded 
2.5  cms .  Also,  the  same  author  stated  that  dowel  action  was  not 
influenced  by  crack  width  or  concrete  cover. 


Upon  the  commencement  of  splitting  along  the  reinforcement, 
dowel  action  deteriorates.  The  level  of  residual  dowel  action  in  such 
circumstances  is  very  much  dependent  on  stirrup  spacing,  but  the 
precise  behaviour  is  difficult  to  define.  In  this  analytical  model, 
dowel  strength  after  splitting  is  considered  negligible.  Dowel  failure 
occurs  when  the  dowel  displacement  ^crack  exceeds  5  x  10~4  inches. 


3.4. 2.5  Cracked  Concrete  Stiffness;  The  concrete  constitutive 
matrix  for  a  cracked  finite  element  is  of  the  form: 
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where  is  the  principal  tensile  stress  that  acts  in  the  direction 
normal  to  the  crack.  Since  the  concrete  stiffness  in  this  direction  has 
been  set  to  zero,  the  stress  in  the  direction  of  will  consequently 
be  zero.  In  the  direction  of  the  principal  compressive  stress,  the 
stress-strain  relationship  is  given  by  Equations  3.27d  and  3.27g  in 
3. 4. 2.1  as  a  uniaxial  loading  condition  now  prevails. 


The  cracked  shear  modulus  G'  is  comprised  of  two  contributions; 
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the  aggregate  interlock  and  dowel  stiffnesses.  The  former  is  readily 
evaluated  at  the  centroid  of  a  concrete  element  as  described  in  Sections 
3. 4. 2. 2  and  3. 4. 2. 3.  In  its  form  presented  in  the  preceding  section, 
the  dowel  force  stiffness  developed  in  the  reinforcement  is  not  in  the 
units  of  shear  modulus.  Therefore,  the  following  procedure  has  been 
adopted.  Upon  calculation  of  the  crack  width  and  centroidal  shear 
strain  parallel  to  the  crack,  the  dowel  displacement  A  is  evaluated. 

ClaC  K. 

The  subsequently  calculated  dowel  force  is  then  considered  uniformly 
distributed  over  the  element’s  cracked  concrete  surface  whose  area  is 
given  by  the  product  of  the  element  thickness  and  the  length  of  the 
inclined  crack  that  passes  through  the  centroid  and  extends  from 
element  boundary  to  boundary.  The  equivalent  dowel  shear  modulus  is 
added  to  the  aggregate  interlock  rigidity  modulus  to  define  the  equivalent 
cracked  concrete  shear  modulus  G' . 

In  the  formulation  of  G’,  the  presence  of  stirrups  is  neglected, 
and  all  dowel  stiffness  contributions  are  calculated  disregarding 
deviation  from  perpendicular  inclination  of  reinforcement  to  the  crack 
direction. 


3. 4. 2. 6  Warping  Resistance  of  Thin  Concrete  Diaphragms:  The 
out-of-plane  warping  resistance  at  the  corners  of  a  thin  concrete  plate 
is  derived  in  most  theory  of  elasticity  texts43.  For  a  square  diaphragm 
element,  the  out-of-plane  warping  stiffness  matrix  is  of  the  form  below: 
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where  =  force  at  node  1  in  a  direction  perpendicular  to  the  element’s 
plane,  X  =  corner  warping  stiffness,  and  =  displacement  at  node  1 
corresponding  to  force  R^. 

Et 3 

X  =  3(1  +  v)  kz  (3.35) 

in  which  E  =  concrete  elastic  modulus,  t  =  thickness,  and  k  =  length  of 
element ' s  diagonal . 

The  off-diagonal  stiffness  terms  in  Eq.  3.34  have  the  same 
magnitude  as  the  main  diagonal  terms  since  the  plate  was  visualized  as 
being  simply  supported  at  each  of  its  four  corners  in  the  derivation  of 
the  warping  restraint  forces.  Rectangular  or  quadrilateral  diaphragm 
elements  are  treated  as  square  elements  of  the  same  area  in  computing 
warping  restraint. 

3. 4. 2. 7  Warping  Restraint  of  Thick  Concrete  Diaphragms: 

Since  the  classical  theory  of  plates  formulation  for  diaphragm  warping 
resistance  given  in  the  preceding  section  is  valid  only  for  "thin" 
plates,  the  distinction  between  "thick”  and  "thin"  plate  thicknesses 
must  be  made.  Theoretically,  no  method  of  distinction  is  currently 
available.  An  additional  qualification  of  the  classical  approach  is 
that  the  derivation  is  developed  for  a  square  plate.  Thus,  an  approxi¬ 
mation  is  immediately  introduced  if  the  formulation  is  applied  to  non¬ 
square  diaphragm  shapes. 

To  analytically  simulate  actual  diaphragm  action,  a  finite 
element  model  was  developed  that  permitted  complete  generality  of 
diaphragm  shape  and  thickness.  To  accommodate  a  complete  range  of 
possible  diaphragm  thicknesses,  an  assemblage  of  three  dimensional 
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bi-quadratic  serendipity  finite  elements  was  used  to  represent  the 
presence  of  diaphragms.  Upon  comparing  the  corner  end  warping  displace¬ 
ments  of  two  identical  double  cell  box  beams  (illustrated  in  Fig.  3.8), 
one  beam  being  restrained  longitudinally  by  end  diaphragms,  the  diaphragm 
warping  stiffness  at  its  four  corners  was  derived  as  a  function  of  the 
box  beam  warping  resistance.  The  warping  resistance  equation  is  of  the 
form: 
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where  R^  =  corner  warping  stiffness  of  diaphragm,  wfc  =  warping 
displacement  of  unrestrained  box  beam  corner,  w,  =  corresponding  warping 


displacement  of  box  beam  corner  restrained  by  diaphragm,  and  R 


wb 


corner  warping  stiffness  of  unrestrained  box  beam. 


The  analytical  model  results  and  the  corresponding  classical 
plate  derivations  for  a  particular  test  beam  are  illustrated  in  Fig.  3.9. 
As  expected,  the  classical  approach  dramatically  overestimates  the 
warping  resistance  for  thick  diaphragms,  but  even  for  smaller  thicknesses, 
there  is  not  good  agreement  between  the  two  formulations.  For 
thicknesses  smaller  than  2  inches,  the  dramatic  divergence  of  the 
approaches  is  not  significant  as  both  predict  negligible  diaphragm  warping 
restraint  compared  with  that  of  the  box  beam.  For  thicknesses  exceeding 
2  inches,  it  is  apparent  that  the  classical  approach  does  not  give 
accurate  correspondence  with  actual  diaphragm  behaviour  as  represented 
with  a  reasonable  degree  of  accuracy  by  the  finite  element  method. 
Consequently,  for  the  beam  cross-sectional  geometry  used  in  this 
comparison,  the  geometry  being  that  of  the  five  rectangular  beams  tested 
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in  the  experimental  program,  the  finite  element  approach  is  the  more 
preferred  method  for  modelling  diaphragm  behaviour. 

The  diaphragm  corner  out-of-plane  warping  stiffnesses  derived 
by  the  above  approach  are  readily  incorporated  in  the  principal  analytical 
computer  model.  Throughout  the  analysis,  it  is  assumed  that  the  diaphragms 
behave  elastically. 

3. 4. 2. 8  Shear  Rigidity  of  Equivalent  Diaphragms:  In  modelling 
a  rectilinear  box  girder  cross-section  as  an  assemblage  of  plane  stress 
finite  elements,  no  account  has  been  taken  of  the  girder's  intrinsic 
cross-sectional  distortion  rigidity.  At  a  cross-section  in  the  analytical 
model  where  an  actual  diaphragm  is  not  provided,  externally  applied  loads 
will  not  be  distributed  correctly  unless  an  equivalent  diaphragm  is 
introduced.  Evaluation  of  the  shear  rigidity  of  an  equivalent  diaphragm 
is  treated  in  detail  by  Sawko  and  Cope41*.  In  the  formulation  of  the 
element  stiffness  matrix,  the  only  non-zero  term  in  the  constitutive 
matrix  is  the  shear  modulus.  Thus,  the  function  of  the  equivalent 
diaphragm  is  twofold  in  preserving  both  cross-sectional  geometry  and 

real  structural  performance. 

3.4.2. 9  Concrete-Reinforcement  Bond ;  Assumption  of  perfect 
bond  between  reinforcement  and  concrete  can  result  in  significant  error45 
in  analytical  modelling.  Thus,  the  concept  of  bond  spring  linkages 
connecting  reinforcement  and  concrete  elements  was  developed,  wherein 
the  loss  of  adhesion  between  concrete  and  steel  is  represented  by  a 
softening  of  the  spring  stiffness.  Nilson45  was  the  first  to  introduce 

a  non-linear  bond-slip  relationship  into  a  finite  element  analysis, 
and  subsequent  extensive  research  investigations  in  this  field,  as  that 
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conducted  by  McCutcheon,  Mirza  and  Mufti49,  used  similar  modelling 
systems  to  that  of  Nilson.  The  adopted  equation  relating  bond  stress 
to  bond  spring  elongation  used  in  stiffness  formulation  is  that  given 
by  Houde  and  Mirza4.  In  both  the  incremental  and  total  load  cases, 
the  bond  spring  linkage  stiffness  is  the  product  of  the  bond  stress  curve 
slope4  and  the  contributing  circumferencial  area  of  the  reinforcement 
bar  to  which  the  bond  linkage  is  attached.  Failure  of  the  bond  spring 
occurs  when  the  spring  elongation  exceeds  .0012  inches. 

3.4.2.10  Reinforcement  Stiffness:  To  simplify  the  analytical 
representation  of  reinforcement  bars,  bi-linear  stress-strain  curves 
have  been  assumed  for  both  conventional  and  prestress  reinforcement, 
as  shown  in  Fig.  3.10. 

Below  yield,  the  stiffness  of  both  reinforcement  types  in  the 
incremental  and  total  load  cases  is  given  by  the  respective  moduli  of 
elasticity.  Beyond  yielding,  the  slope  of  the  strain  hardening  segment 
of  the  stress-strain  curve  defines  the  incremental  reinforcement  stiffness. 
However,  in  the  total  load  condition  that  prevails  in  the  iterative 
process,  the  method  of  stiffness  derivation  for  the  conventional  and 
prestress  reinforcement  differs.  Since  the  slope  of  the  strain  hardening 
section  of  the  conventional  reinforcement  stress-strain  curve  is  highly 
inelastic,  an  excessive  number  of  modified  Newton-Rapson  iterations  are 
often  required  to  restore  equilibrium.  If  the  set  of  equations  were 
large,  this  approach  could  be  highly  impractical.  Therefore,  once  a 
conventional  reinforcement  bar  has  yielded  in  the  total  load  condition, 
it  is  represented  in  the  analytical  model  by  a  bar  of  zero  stiffness, 
the  load  carried  by  the  bar  being  represented  by  equivalent  external 
loads  applied  at  the  bar  nodes.  Following  such  a  stiffness  change. 
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The  strain  in  the  yielding  bar  will  increase  as  the  load  is  maintained, 
but  the  magnitude  of  the  equivalent  bar  loads  will  increase  only  slightly. 
Thus,  at  a  particular  load  level,  no  iteration  is  required  to  restore  an 
acceptable  degree  of  equilibrium  to  the  yielding  bars.  This  approach 
cannot  be  applied,  however,  to  the  yielding  of  the  prestress  reinforcement. 
Prestress  reinforcement  constitutes  the  underreinforced  beam’s  last 
reserve  of  strength,  and  setting  of  its  stiffness  to  zero  will  immediately 
produce  instability.  Thus,  the  modified  Newton-Rapson  method  is  retained 
for  modelling  yielding  prestress  reinforcement  in  the  total  load 
condition.  Use  of  a  relaxation  factor  has  been  introduced  to  improve 
the  rate  of  convergence  of  prestress  strand  deviation. 

3.4.3  Failure  Criteria 

Of  the  two  material  constituents,  concrete  and  reinforcing 
steel,  the  failure  characteristics  of  concrete  are  complex  and  warrant 
clarif icat ion.  In  establishing  the  tensile  and  compressive  strengths 
of  concrete  under  biaxial  stress  conditions,  Kupfer,  Hilsdorf,  and 
Rusch46  have  postulated  a  biaxial  stress  envelope  illustrated  in  Fig. 

3.11.  In  addition  to  specifying  biaxial  stress  failure  combinations, 
the  authors  also  established  that  the  Poisson  ratio  in  biaxial  tension, 
biaxial  compression,  and  tension-compression  were  .18,  .2,  and  an 
average  value  of  .19  respectively.  The  simplified  envelope  used  in  the 
analytical  model  is  shown  in  Fig.  3.12.  Beyond  cracking,  the  shear 
rigidity  across  a  concrete  crack  is  developed  by  the  aggregate  interlock 
and  dowel  mechanisms.  In  the  absence  of  relevant  information,  the 
failure  stress  level  for  the  aggregate  interlock  mechanism  has  been  set 
equal  to  the  uniaxial  concrete  compressive  strength.  Under  actual  test 
conditions,  the  evaluated  aggregate  interlock  stress  has  not  been  found 
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to  rise  beyond  750  p.s.i.4.  Dowel  linkage  and  the  concrete-related 
bond  spring  linkage  failure  criteria  are  specified  in  preceding  Sections 
3. 4. 2. 6  and  3. 4. 2. 9  respectively. 

In  a  total  structure  context,  failure  in  the  analytical  model 
is  defined  as  having  occurred  when  a  concrete  element  has  crushed  or  a 
major  prestress  reinforcement  element  has  failed  in  tension. 

3.4.4  Numerical  Methods 

The  three  numerical  methods  that  perform  vital  functions  in 
the  analytical  model  are  the  block-by-block  gaussian  elimination  solution 
process,  Runge-Kutta  and  modified  Newton-Rapson  methods. 

To  reduce  the  size  of  the  allocated  core  storage  area  used  in 
the  equation  solution  process,  only  two  stiffness  blocks  are  retained 
in  core  at  any  one  time.  After  a  stiffness  block  has  been  reduced  by 
gaussian  elimination,  it  is  written  onto  auxiliary  disc  storage,  and 
the  following  unreduced  stiffness  block  is  subsequently  read  into  core. 

If  an  efficient  method  of  transfer  of  the  large  data  sets  in  and  out  of 
core  storage  is  used,  this  block-by-block  solution  process  is  ideally 
suited  to  the  solution  of  large  equation  systems.  Logic  details  are 
given  in  the  listing  of  subroutine  SOLVE  in  Appendix  D. 

The  area  of  application  of  the  Runge-Kutta  method  is  the 
progressive  adjustment  of  the  structural  stiffness  as  the  external 
loads  are  applied  incrementally.  Figure  3.13(a)  illustrates  the  procedure. 
Using  the  stiffness  of  the  material  element  in  the  previous  (n-l)th. 
load  increment,  defined  by  the  slope  of  line  AB,  the  stress-strain  state 
upon  application  of  half  of  the  current  nth.  load  increment  is  calculated. 
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denoted  by  point  C.  The  tangent  stiffness  at  the  corresponding  point  D 
on  the  predefined  stress-strain  curve  for  the  material  is  subsequently 
adopted  as  a  reasonable  prediction  of  the  average  material  element 
stiffness  for  the  nth.  load  increment.  Upon  superposition  of  the  full 
nth.  load,  the  stress-strain  state  at  the  end  of  the  nth.  increment  is 
established,  denoted  by  point  E.  Thus,  the  Runge-Kutta  method  of 
stiffness  prediction  improves  program  efficiency  through  minimizing  the 
need  for  the  time  consuming  iterative  process  that  must  be  invoked  if 
significant  material  deviation  occurs. 

Should  deviation  at  the  end  of  a  load  increment  be  significant, 
equilibrium  is  restored  by  employing  the  iterative  modified  Newton- 
Rapson  method.  After  reducing  the  deviation  of  point  Q  in  Fig.  3.13(b) 
to  an  acceptable  level,  denoted  by  point  R,  the  tangent  stiffness  at 
point  S  is  the  initial  estimate  of  increment  stiffness  used  in  the 
Runge-Kutta  stiffness  adjustment  for  the  nth.  load  increment. 

Both  the  Runge-Kutta  and  modified  Newton-Rapson  techniques 
are  treated  in  detail  in  most  applied  mechanics  texts47  that  address 
non-linear  behaviour. 

3.4.5  Flow  Chart  for  Main  Program 

In  the  flow  chart  that  follows,  the  logic  skeleton  of  the 
main  program  is  presented,  the  symbolic  names  in  brackets  designating 
the  subroutines  that  accomplish  the  defined  logic  steps.  Further  logic 
clarification  is  supplied  by  the  numerous  comment  cards  distributed 
through  the  main  program,  whose  listing  is  given  in  Appendix  C. 
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FLOW  CHART  FOR  MAIN  PROGRAM 
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49. 


3.4.6  Derivation  of  Subroutine  Logic 

Logic  development  notes  of  the  more  complex  subroutines 
comprise  Appendix  E. 

3.5  Program  Usage 

3.5.1  Capabilities  and  Restrictions 

(A)  Beam  Geometry 

The  program  has  been  developed  for  the  analysis  of  concrete 
box  girder  structures  with  modest  wall  thicknesses,  but  any  wall  thickness 
can  be  accommodated  with  the  qualification  that  accuracy  of  solution 
will  be  prejudiced  for  "thick"  walls,  as  detailed  in  3.2.1.  Cross- 
sectional  geometry  can  assume  any  form  that  can  be  represented  by  linear 
segments,  but  in  the  program's  present  form,  no  geometry  change  is 
permitted  along  the  beam's  length.  Such  a  shortcoming  can  be  remedied 
by  replacing  the  rectangular  concrete  finite  element  by  a  general 
quadrilateral  element38. 

(B)  Material  Behaviour 

To  reduce  computational  effort,  the  behaviour  of  a  concrete 
element  is  characterized  by  the  stress-strain  condition  at  the  element's 
centroid.  For  modestly  fine  meshes  of  high  order  elements,  this 
approach  is  an  acceptable  approximation.  Uncracked  concrete  is  modelled 
as  a  heterogeneous,  anisotropic,  non-linear  material,  but  cracking 
introduces  two  phenomena  associated  with  the  definition  of  shear  rigidity 
across  a  crack,  namely  aggregate  interlock  effect  and  dowel  action. 

As  cracking  progresses,  the  stiffness  of  these  two  contributions  are 
modified  in  addition  to  the  concrete  stiffness  in  the  parallel  crack 
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direction.  The  closing  of  cracks  produced  by  dramatic  load  reversal 
in  the  load  increment  following  prestress  transfer,  and  the  superposition 
of  a  different  loading  type  (torque  upon  bending  load)  is  also  within 
the  program's  capability. 

The  two  principal  reinforcement  types,  conventional  steel 
and  prestress  strand,  are  modelled  by  one-dimensional  finite  elements. 
Consequently,  a  continuous  reinforcement  bar  is  represented  by  a  series 
of  segments  of  constant  stress,  producing  a  step-like  variation  in  bar 
force.  In  direct  contrast,  steel  mesh  reinforcement  is  considered  to 
be  uniformly  distributed  throughout  a  concrete  element,  and  its  stiffness 
is  formulated  in  an  identical  manner  to  that  for  a  rectangular  concrete 
finite  element.  The  stress-strain  curves  for  all  three  reinforcement 
types  are  approximated  by  multi-linear  functions. 

In  the  development  of  the  finite  element  mesh,  the  mechanism 
of  force  transfer  at  the  steel-concrete  interface  is  simulated  by  bond 
spring  linkages  connecting  adjacent  steel  and  concrete  nodes.  As  loading 
progresses,  the  stiffness  of  the  bond  spring  is  adjusted  to  reflect 
the  gradual  deterioration  in  concrete  bond. 

Should  significant  material  behavioural  deviation  be  detected 
within  a  load  increment,  an  iterative  process  is  automatically  initiated, 
and  proceeds  until  equilibrium  is  restored. 

(C)  Loading  and  Boundary  Conditions 

Any  loading  or  boundary  condition  combination  can  be 
imposed  upon  the  analytical  model  provided  that  the  appropriate  degrees 
of  freedom  are  present.  In  specifying  load  magnitudes,  the  subsequent 
number  of  load  increments  should  not  be  small  such  that  the  computationally 


51. 


expensive  iterative  process  is  invoked  frequently.  Provision  for  load 
application  is  not  restricted  to  external  test  loads,  since  both 
prestressing  and  post-stressing  techniques  can  be  represented 
by  equivalent  model  forces. 

(D)  Diaphragms 

The  presence  of  diaphragms  of  any  thickness  and  the  intrinsic 
cross-sectional  deformation  resistance  of  the  hollow  beam  are  two 
important  components  of  box  girder  behaviour,  and  both  characteristics 
are  incorporated  in  the  program  logic.  Throughout  the  entire  load  range, 
diaphragm  action  is  assumed  to  be  elastic. 

(E)  Printout  Information 

Through  the  use  of  output  control  parameters,  the  precise 
nature  of  the  printed  output  can  be  specified  in  a  selective  qualitative 
and  quantitative  manner,  as  detailed  in  3.5.4. 

(F)  Program  Monitoring 

Provision  is  made  within  the  program  for  transmission  of 
output  information  to  a  monitoring  device  such  as  a  CRT  display  terminal. 
In  the  running  of  large  problems,  the  user  has  the  option  of  observing 
program  progress  and  controlling  the  rate  of  execution  in  the  event 
that  erroneous  results  may  prompt  run  termination,  thus  reducing 
unnecessary  complete  run  expense. 

3.5.2  Structure  Discretization 

Figure  3.13  illustrates  the  finite  element  mesh  chosen  to 
model  three  double-cell  prestressed  concrete  rectangular  box  beams 
tested  in  the  experimental  program.  In  specifying  structural  geometry, 
the  global  x  axis  must  be  parallel  to  the  beam's  longitudinal  centroidal 
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3.5.3  Input  Specifications 

In  the  preparation  of  the  input  data  file,  reference  should 
be  made  to  Appendix  H  in  which  a  detailed  card-by-card  description  of 
the  input  data  is  given,  together  with  implicit  logic  assumptions  that 
have  a  direct  bearing  on  data  derivation.  For  ease  of  preparation,  the 
input  format  is  semi-free  field,  as  illustrated  in  Appendix  F  which 
lists  a  sample  input  data  file. 

3.5.4  Output  Description  and  Interpretation 

The  form  of  the  printed  output  is  governed  by  the  choice  of 
the  output  control  variables  described  in  Appendix  H.  Independent  of 
specified  controls,  increment  headings  are  printed,  together  with  local 
material  failure  messages  and  progressive  CPU  and  program  cost  estimates. 
For  each  load  increment,  structure  deformations,  concrete  centroidal 
stresses  and  strains  in  the  global  and  principal  axes  directions,  steel 
mesh  and  bar  reinforcement  stresses  and  strains,  and  increment  and  total 
load  levels  can  be  selectively  printed  on  an  element  or  category  basis. 
Units  of  pounds  wt.,  inches,  and  radians  are  used  throughout. 

In  the  interpretation  of  output  deformations,  all  values 
correspond  to  displacements  in  structural  degree  of  freedom  directions. 

The  results  must  be  processed  further  to  determine  cross-sectional 
rotation,  beam  curvature,  and  other  such  descriptive  deformations. 

For  both  concrete  and  steel  mesh  elements,  all  stress  and  strain 
estimates  are  derived  for  the  centroidal  element  location.  Since  the 
reinforcement  bar  element  is  one-dimensional,  the  respective  stress-strain 
values  are  constant  for  the  element's  length.  Structural  failure  is 
defined  as  occurring  when  a  concrete  element  crushes  or  a  prestress 
strand  fails  in  tension.  Thus,  the  failure  statement  that  immediately 
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precedes  run  termination  is  only  of  two  possible  forms.  However, 
structural  failure  can  develop  through  instability,  such  a  failure  mode 
resulting  in  the  formation  of  an  ill-conditioned  set  of  equilibrium 
equations.  Under  such  a  situation  program  progress  ceases  when  a 
negative  term  is  encountered  on  the  main  diagonal  of  the  total  stiffness 
matrix  during  the  reduction  process. 

Since  the  success  of  the  finite  element  method  of  analysis 
is  reflected  in  the  quality  of  the  program  input,  careful  consideration 
should  be  given  to  the  evaluation  of  the  numerous  strength  of  materials 
parameters,  in  particular  those  that  have  a  strong  influence  on  post¬ 
cracking  behaviour. 
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TABLE  3.2(B)  [A]-1  FOR  ELEMENT  TYPE  2 


FIG.  3.1  SIMPLIFIED  REPRESENTATIVE  FINITE  ELEMENT  MESH 


(b)  Force  couple  due  to  open-cell 
shear  stress  distribution 


FIG.  3.2  NOTATION  FOR  DERIVATION  OF  ST.  VENANT  SHEAR 
STRESS  DISTRIBUTION 


FIG.  3.3  BI-LINEAR  ISOPARAMETRIC  SERENDIPITY  FINITE 
ELEMENT 
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(d)  Location  of  loading  cross-sections 


FIG.  3.4  THREE  FINITE  ELEMENT  MESHES  USED  IN  CONCRETE  SELECTION 
PROCESS 
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FIG.  3.5(A)  DIRECT  BENDING  STRESSES  FOR  THREE  CONCRETE  FINITE 
ELEMENTS 


- - * 


SCORDELIS 

McCLEOD 

SISODIYA  AND 
GHALI 


2  4  6  8  |0 

NUMBER  OF  ELEMENTS  IN  LENGTH 


12 


FIG.  3.5(B)  ST.  VENANT  TORSION  SHEAR  STRESSES  FOR  THREE  CONCRETE 
FINITE  ELEMENTS 
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FIG.  3.6(B)  ELEMENT  TYPE  2 


FIG.  3.6(C)  TWO-DIMENSIONAL  BEAM  ASSEMBLAGE  OF  McCLEOD 
FINITE  ELEMENTS 


FIG.  3.7  THEORETICAL  AND  EXPERIMENTAL  CONCRETE  COMPRESSION 
STRESS-STRAIN  CURVES 


/classical*  diaphragm  warping  resistance 

FINITE  ELEMENT  DIAPHRAGM  WARPING  RESISTANCE 
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FIG. 3.9  DIAPHRAGM  WARPING  RESTRAINT  CURVES 
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FIG.  3.10(A)  CONVENTIONAL  REINFORCEMENT  STRESS-STRAIN  CURVE 


FIG.  3.10(B)  PRESTRESS  REINFORCEMENT  STRESS-STRAIN  CURVE 
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FIG. 3. 11  KUPFER,  HILSDORF,  RUSH  BIAXIAL  STRESS  ENVELOPE 


FIG.  3.12  SIMPLIFIED  ANALYTICAL  MODEL  BIAXIAL  STRESS 
ENVELOPE 
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STRAIN 


FIG.  3.13(A)  RUNGE-KUTTA  METHOD 


STRAIN 


FIG.  3.13(B)  MODIFIED  NEWTON  RAPSON  METHOD 
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FIG.  3.14  FINITE  ELEMENT  MESH  DISCRETIZATION 


CHAPTER  4 


EXPERIMENTAL  PROGRAM 


4 . 1  Introduction 

Consistent  with  the  defined  scope,  the  scale  of  the 
experimental  program  was  modest  in  the  number  of  beam  specimens  tested 
and  the  degree  of  parameter  variation.  In  all,  seven  prestressed  box 
beams  were  cast  and  tested  to  failure,  with  each  beam  subjected  to 
different  ratios  of  torque  to  bending  moment  to  shear.  All  facets  of 
the  experimental  program  are  dealt  with  in  this  chapter,  and  the  results 
are  presented  in  the  following  chapter. 

4 . 2  Definition  of  Basic  Experimental  Parameters 

In  essence,  all  the  test  beam  specimens  were  hollow,  precast, 
prestressed  concrete  members.  At  the  preliminary  design  stage,  it  was 
recognized  that,  if  the  analytical  model  was  to  be  tested  rigorously, 
complex  beam  cross-sections  of  more  than  one  cell  should  be  tested. 

An  extensive  literature  survey  revealed  that  little  if  any  experimental 
research  had  been  conducted  on  the  post-cracking  behaviour  of  multi-celled 
members  under  the  action  of  torsion,  bending,  and  shear.  From  an 
analytical  stance  of  generality,  little  benefit  would  have  been  derived 
in  having  more  than  two  cells  per  beam,  and  practical  expedience  in  the 
casting  and  testing  sequences  dictated  against  an  excessive  number  of 
voids.  Thus,  the  decision  was  made  that  all  beams  would 'be  cast  with 
two  closely  spaced  voids. 
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Since  the  experimental  tests  were  to  be  used  strictly  for 
comparative  purposes,  the  number  of  beams  cast  for  the  testing  program 
was  small.  The  only  principal  experimental  parameter  whose  variation 
was  investigated  in  the  series  of  tests  was  the  load  parameter.  The 
more  significant  aspect  of  this  parameter  was  not  the  level  of  individual 
loads,  but  the  ratio  of  bending  moment  to  torque  to  shear.  In  all, 
seven  beams  were  cast,  with  five  beams  being  subjected  to  varying  ratios 
of  torque  to  bending  moment  with  little  shear  present.  The  remaining 
two  beams  were  subjected  to  significant  levels  of  shear,  in  conjunction 
with  the  more  predominant  torque-bending  moment  loading  combination. 

The  primary  region  of  concern  in  the  full  range  of  test 
beam  behaviour  was  the  post-cracking  region.  Thus,  the  initial  design 
of  the  test  specimens  had  to  meet  the  important  prerequisite  that  the 
inelastic  behavioural  path  beyond  cracking  was  of  reasonable  length 
compared  with  the  respective  elastic  path.  This  condition  was  assured 
by  not  permitting  the  ratio  of  ultimate  strength  to  cracking  strength 
to  fall  much  below  the  value  of  two. 

To  introduce  a  degree  of  generality  of  cross-sectional  shape 
in  the  experimental  study,  two  geometrical  shapes  were  adopted  in  the 
beam  design;  rectangular  and  trapezoidal.  Of  the  total  of  seven  beams, 
five  were  cast  with  rectangular  sections. 

The  one  overriding  consideration  in  the  selection  of  void 
sizes  was  that  all  wall  segments  must  be  classified  as  "thin".  The 
need  for  thin  wall  segments  arose  from  the  use  of  plane  stress  finite 
elements  to  model  the  beam  walls.  The  theoretical  classification  as  to 
whether  a  beam  cross-section  is  "thick"  or  "thin",  is  treated  in  detail 


in  3.2.1. 
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4.3  Design  of  Test  Specimens 

All  beams  in  the  test  program  were  precast,  prestressed 
concrete  members  with  modest  conventional  reinforcement  present. 

Beams  were  categorized  as  belonging  to  either  of  the  two  test  beam 
series  according  to  their  cross-section  shape,  those  rectangular  beams 

belonging  to  the  "R"  series  and  those  trapezoidal  assigned  to  the  "T" 
series. 

In  all  seven  beams,  a  design  concrete  strength  of  5,000  p.s.i. 
was  adopted.  To  achieve  the  appropriate  workability  that  would  allow 
the  small  wall  thicknesses  to  be  cast  free  of  air  voids,  the  mix  design 
was  made  considerably  richer  than  was  generally  required,  the  actual 
concrete  component  proportions  being  given  in  Table  4.1.  The  respective 
compressive  and  tensile  concrete  strengths  for  the  seven  beams  are  given 
Table  4.2.  The  load  versus  elongation  curves  for  the  remaining  two 
material  components,  the  prestress  strand  and  the  conventional  reinforce¬ 
ment,  are  displayed  in  Figures  4.1  and  4.2  respectively.  As  a  point 
of  clarification,  the  prestress  strand  was  250K  grade,  1/4"  diameter 
seven  strand  stress  relieved  cable,  whereas  the  conventional  reinforce¬ 
ment  used  throughout  was  #3  deformed  mild  steel. 

In  defining  full  beam  geometry,  beam  length,  wall  thickness, 
number  of  cells,  cross-sectional  shape  and  dimensions  have  to  be 
specified.  Of  the  three  variables  that  are  dimensionally  functions  of 
length,  the  independent  variable  was  the  wall  thickness.  This  beam 
dimension  had  to  be  sufficiently  large  such  that  the  prestress  and 
conventional  reinforcement  could  be  accommodated  without  drastically 
impeding  the  flow  of  concrete  during  casting.  As  mentioned  previously. 
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a  wall  thickness  of  2"  was  chosen  for  all  beams.  Having  specified  the 
wall  thickness,  the  dependent  variables  of  cross-section  dimensions  and 
beam  length  were  determined  upon  recognition  of  their  respective 
constraints.  The  web  height  and  flange  width  had  to  be  a  considerable 
order  of  magnitude  larger  than  the  wall  thickness  to  ensure  a  reasonable 
degree  of  plane-stress  action  of  the  concrete  beam  walls,  and  the  beam 
length  such  that  a  central  test  section  for  deformation  measurement  and 
the  appropriate  loading  and  support  apparatus  could  be  accommodated. 

The  remaining  two  geometrical  variables  of  cross-sectional  shape  and 
number  of  cells  have  been  treated  in  Section  4.2.  Complete  beam 
dimensioning  is  given  in  Figures  4.3  and  4.4. 

In  design,  the  appropriate  level  of  prestress  was  chosen 
such  that,  at  transfer  four  days  after  casting,  the  maximum  allowable 
tensile  stresses  were  developed  in  the  top  flange.  In  all  beam  designs, 
the  bottom  flange  compressive  stresses  did  not  approach  the  allowable 
limits. 

In  the  calculation  of  the  ultimate  bending  moment  capacity, 
the  level  of  reinforcement  was  such  that  the  concrete  crushed  just 
prior  to  the  prestress  strand  developing  its  ultimate  strength.  In 
establishing  the  ultimate  torsional  capacity,  the  space  truss  analogy 
equations  as  presented  by  Collins  and  Lampert25  were  used  to  estimate 
beam  strength. 

As  mentioned  in  the  previous  section,  the  three  distinct 
loading  types  represented  in  this  experimental  program  were  bending  moment, 
torque,  and  shear.  Of  the  five  rectangular  beams  cast,  three  were 
subjected  to  combined  bending  and  torsion  only,  the  remaining  two  acted 
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upon  by  complete  bending-torsion-shear  loading  combinations.  The  two 
trapezoidal  beams  were  restricted  to  bending  and  torsion  loads.  In 
predicting  beam  strength  under  combined  loading,  the  interactive  equations 
given  in  the  Collins  and  Lampert  paper  were  utilized,  such  equations 
given  in  Fig.  4.5  together  with  the  design  curves  for  the  rectangular 
and  trapezoidal  beams.  The  predicted  beam  strengths  and  failure  loads 
are  given  in  Table  4.3. 

Premature  shear  failure  was  prevented  in  the  beam  length 
outside  the  central  test  section  by  the  provision  of  substantial  shear 
reinforcement,  achieved  in  most  beams  by  extending  torsion  hoop  rein¬ 
forcement  beyond  the  central  test  section.  Reinforcement  details  are 
provided  in  Figures  4.3  and  4.4.  Also  shown  in  the  latter  figures  are 
the  locations  of  the  points  of  application  of  the  bending,  torsion 
and  shear  loads. 

4 . 4  Specimen  Fabrication 

To  enable  both  rectangular  and  trapezoidal  beams  to  be  cast 
with  the  same  forming  unit,  the  timber  formwork  consisted  of  two 
identical  reversible  form  segments  placed  adjacent  to  one  another,  as 
shown  in  Fig.  4.6. 

Provision  of  the  two  adjacent  cells  within  each  beam  was 
simply  achieved  by  the  use  of  styrofoam  blocks  held  in  place  by  the 
reinforcement  cage.  In  previous  research  conducted  at  this  university, 
the  styrofoam  was  removed  by  piping  acetone  into  the  interior  of  the 
cast  beams.  This  added  precaution,  however,  was  not  deemed  necessary 
as  the  styrofoam  strength  contribution  was  negligible  compared  with 
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that  of  the  prestressed  beam.  In  all  beams,  the  void  lengths  did  not 
extend  as  far  as  the  points  of  support,  as  exhibited  in  Figure 
4.4. 


As  indicated  above,  the  presence  of  the  reinforcement  cage 
was  utilized  to  maintain  the  positions  of  the  styrofoam  blocks.  To 
prevent  the  styrofoam  blocks  from  floating  during  cast,  restraining 
beams  were  placed  across  the  top  of  the  forms  to  inhibit  upward  cage 
movement.  Since  the  cage  tolerances  were  quite  small  arising  from  the 
selection  of  narrow  wall  thicknesses,  the  reinforcement  bars  and  hoops 
were  lightly  tack-welded  in  position  to  minimize  movement  and  facilitate 
handling.  As  a  matter  of  necessity,  the  two  top  corner  longitudinal 
bars  had  to  be  moved  1"  from  their  corner  positions  to  allow  a  small 
vibrating  needle  to  pass  down  the  outside  web  walls  during  casting. 

Despite  the  small  wall  thickness,  the  rich  concrete  mix  whose 
proportioning  is  given  in  Table  4.1  did  exhibit  good  workability  with 
the  result  that  very  few  air  voids  were  evident  upon  examination  after 
formwork  stripping. 

The  complete  casting  bed  before  concrete  pouring  is  shown 
in  Plate  4.1. 

4 . 5  Loading  Apparatus  Design 

4.5.1  Beam  Supports 

In  addition  to  their  customary  role  of  simple  supports,  the 
beam  reactions  permitted  the  development  of  a  uniform  St.  Venant  torque 
applied  across  the  central  test  section.  Consequently,  each  support  had 
to  allow  the  beam  to  rotate  freely  about  its  centre  of  rotation.  Lateral 
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stability  of  the  beam  was  ensured  by  the  very  nature  of  the  torsion 
loading  arms,  as  described  in  Section  4.5.3.  Figure  4.7  illustrates 
all  end  support  details.  The  roller  housing  shown  in  this  figure  was 
bolted  to  the  top  of  the  conventional  roller  and  hinge  supports  to 
achieve  the  desired  beam  support  conditions. 

4*5.2  Application  of  Bending  Moment  and  Shear 

Both  bending  moment  and  shear  force  were  developed  by  the 
application  of  downward  vertical  concentrated  loads  applied  by  50  ton 
Amsler  jacks.  If  a  vertical  concentrated  load  was  not  to  create  secondary 
torsion  during  the  course  of  a  beam  test,  the  centre  line  of  the  Amsler 
jack  had  to  pass  through  the  shear  centre  of  the  beam.  This  alignment 
was  achieved  through  incorporating  a  design  that  was  identical  to  that 
of  the  rotational  end  support,  with  the  apparatus  simply  inverted.  The 
one  additional  design  provision  not  present  in  the  end  support  design 
was  the  roller  housing  bracing  whose  sole  function  was  to  provide 
horizontal  stability  to  the  roller  housing  under  all  loading  extremes. 

As  the  beam  deflected  under  load,  the  bracing  arms  attached  to  the  housing 
by  a  central  pin  maintained  their  horizontal  posture  upon  adjustment 
of  the  bracing  jack.  The  apparatus  is  displayed  in  Fig.  4.8. 

4.5.3  Application  of  Torque 

Only  the  central  test  section  was  subjected  to  torque,  more 
specifically  a  uniform  St.  Venant  torque  developed  by  a  pair  of  equal 
and  opposite  forces  whose  placement  along  the  test  beam  delineated 
the  central  test  section.  Each  torsion  force  was  applied  by  a  cable 
that  draped  over  a  curved  torsion  arm  (as  in  Fig.  4.9)  and  passed 
through  the  testing  floor  where  it  was  anchored  by  a  hydraulic  jack 
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loading  system.  The  curved  arm  geometry  ensured  that  the  torsion  lever 
arm  dimension  remained  constant  as  the  beam  rotated.  Stability  of  a 
test  beam,  once  the  beam  was  placed  in  the  testing  arrangement,  was 
provided  by  minimal  tension  in  the  torsion  cables. 

4 . 6  Instrumentation 

To  plot  beam  behaviour  throughout  the  test  program,  both  stress 
and  deformation  data  were  recorded.  As  the  seven  strand  prestress 
cable  was  of  small  diameter  (1/4”) ,  strain  gages  were  difficult  to 
attach  securely,  and  thus  could  not  be  relied  upon  to  accurately 
monitor  prestress  strain  levels.  However,  the  strain  states  of  the 
conventional  longitudinal  and  central  transverse  hoop  reinforcement 
were  monitored  by  strain  gages  to  yield  a  representative  record  of  stress 
versus  load  behaviour. 

More  emphasis  was  placed  in  the  instrumentation  phase  on  the 
accurate  measurement  of  cross-section  deformation  within  the  central 
test  section  rather  than  extensive  monitoring  of  reinforcement  strain. 

By  surveying  the  vertical  and  horizontal  deflections  of  three  cross- 
sections  within  the  central  test  section,  the  central  vertical  deflection, 
beam  curvature,  and  rate  of  twist  were  able  to  be  evaluated.  To 
facilitate  fast  test  measurement  recording,  linear  displacement  transducers 
were  used  in  preference  to  conventional  dial  gages.  Plate  4.2  displays 
several  linear  displacement  transducers  operating  under  test  conditions. 

In  recognition  of  the  fact  that  the  number  of  both  load 
increments  and  individual  instrument  readings  was  large,  an  automatic 
data-recording  system  was  used.  The  electrical  signals  from  the  linear 
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displacement  transducers  and  the  strain  gages  were  monitored,  trans¬ 
lated,  and  stored  on  disc  by  a  Data  General  Corporation  Nova  210E 
Computer.  During  testing,  the  depression  of  a  console  key  prompted 
the  complete  data  set  of  load  increment  measurements  to  be  recorded  and 
written  on  disc. 

All  locations  of  strain  gages  and  linear  displacement  transducers 
are  described  in  Chapter  5. 

4 . 7  Testing  Procedure 

Upon  application  of  each  load  increment,  the  entire  regime  of 
instrument  measurement  was  recorded  automatically  as  described  in  the 
preceding  section.  If  cracking  had  commenced,  crack  propagations  were 
clearly  marked  and  the  load  intensity  indicated  at  the  furthest  point 
of  crack  propagation,  as  seen  in  the  following  chapter's  photographic 
plates  of  tested  beams.  As  a  precautionary  measure,  the  transverse 
alignment  of  the  Amsler  jacks  was  checked  at  regular  intervals  and 
corrected  when  necessary  by  adjustment  of  the  bracing  turnbuckles  and 
jacks.  All  cracking  patterns  and  mode  of  failure  were  photographed 
upon  completion  of  the  test. 
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Ingredient 

Wt.  in  lbs. /batch 

Cement 

211 

Sand 

420 

Coarse  Aggregate 
(3/8") 

540 

Water 

120 

Approximate  batch  volume  =  10  ft.3 

slump  =  3"  ■+  4" 

TABLE  4.1  CONCRETE  MIX  PROPORTIONS 


Beam 

Designation 

Concrete 
Compressive 
Strength  (psi) 

Concrete 
Tensile 
Strength  (psi) 

Concrete 
Young  * s 
Modulus  (psi) 

R1 

4816 

396 

2.36 

x  106 

R2 

4580 

462 

2.133 

x  10G 

R3 

4545 

479 

2.6 

x  106 

R4 

4362 

288 

2.69 

x  106 

R5 

4562 

329 

1.84 

x  10G 

T1 

3985 

354 

3.0 

x  10G 

T2 

4262 

338 

2.79 

x  10G 

TABLE  4.2  CONCRETE  STRENGTH 
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TABLE  4.3  BEAM  DESIGN  TABULATION 
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FIG.  4.1  LOAD-ELONGATION  CURVE  FOR  PRESTRESS  STRAND 
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(!8M)  SS3U1S 


FIG.  4.2  STRESS-STRAIN  CURVE  FOR  #3  DEFORMED  BAR  REINFORCEMENT 
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FIG.  4.3(A)  RECTANGULAR  BEAM  CROSS-SECTION 


FIG.  4.3(B)  TRAPEZOIDAL  BEAM  CROSS-SECTION 

Notes:  (1)  All  wall  thicknesses  =  2" 

(2)  O  Conventional  #3  bar  •. 25 "  Prestress  strand 

(3)  All  beams  are  symmetrical  w.r.t.  vertical  centreline 


79. 


B  I'.il" 
'  '  2 


LONGITUDINAL  #  3  BARS 


-  I  M 
2-0 


— . ■■■■  V - - - - - - - 

r 

1 

1 

S’-S" 

• 

VOID 

1 1 

L 

MCK' 


18 


STIRRUPS-*- 


S  =4 


S=  8 


■2  at  2 


i 


HOOPS 


Note:  Total  void  length  for  beam  R1  =  10'-  0" 


FIG.  4.4(A)  REINFORCEMENT  FOR  BEAMS  R1  AND  R2 


Notes:  (1)  All  bars  are  #3 

(2)  All  half  beam  lengths  =  7'-  7.5" 

(3)  R  =  support  B  =  jack  T  =  torsion 

(4)  All  void  lengths  =11'-  0"  unless  noted  otherwise 

(5)  All  longitudinal  steel  is  1.5"  short  of  beam  ends 

(6)  Drawings  are  not  to  scale 
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FIG.  4.5  INTERACTION  EQUATIONS  FOR  TEST  BEAMS  UNDER  COMBINED 
BENDING  AND  TORSION 
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Juxtiposition  for  casting  rectangular  beams 


FIG.  4.6  FORMWORK  SECTION 
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FIG.  4.7  TEST  BEAM  SUPPORTS 
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FIG.  4.9  TORSION  LOAD  ARM 
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PLATE  4.1  FORMWORK  BEFORE  CASTING  PLATE  4.2  LINEAR  DISPLACEMENT  TRANSDUCERS 
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PLATE  4.3  COMPLETE  EXPERIMENTAL  TESTING  SETUP 


CHAPTER  5 


EXPERIMENTAL  RESULTS 

5-1  Format  of  Presentation  of  Results 

The  complete  spectrum  of  experimental  results  for  the  seven 
beams  tested  is  presented  in  this  chapter.  The  two  distinct  aspects 
of  the  testing  program,  test  measurements  and  observed  behaviour,  are 
treated  separately  in  the  following  Sections  5.2  and  5.3. 

In  the  measurement  of  beam  response,  greater  emphasis  was 
placed  on  the  recording  of  member  deformation  than  reinforcement  stress 
levels.  For  each  beam,  graphs  of  torque  versus  differential  rotation 
and  bending  moment  versus  central  deflection  describe  complete 
deformation  characteristics.  In  the  only  representation  of  monitored 
stresses,  stress  levels  of  the  longitudinal  and  transverse  conventional 
reinforcement  at  the  central  cross-section  are  presented  in  the 
respective  plots.  Beam  strength  is  tabulated  fully  in  the  specification 
of  prestress  levels,  elastic  stiffness,  cracking  and  failure  loads. 

Experimental  observations  are  confined  to  modes  of  failure 
and  distinctive  cracking  patterns. 

In  the  concluding  section  of  this  chapter,  the  potential 
sources  of  experimental  error  and  irregularities  in  the  experimental 
program  are  identified  and  evaluated. 

A  complete  detailed  tabulation  of  experimental  results  is 
provided  in  Appendix  I. 


-  85  - 


86. 


5 . 2  Test  Measurements 

The  recorded  numerical  data  in  its  processed  form  is  presented 
under  the  three  categories  of  beam  strength,  beam  stiffness,  and 
reinforcement  stress. 

The  beam  strength  components  of  initial  level  of  prestress, 
final  level  of  prestress,  cracking  load,  and  loading  combination  at 
failure  are  summarized  in  Table  5.1.  At  both  the  cracking  and  failure 
values  the  bending  moment  specified  is  the  total  bending  moment  arising 
from  the  torsion  arm  and  Amsler  jack  loads. 

To  provide  a  suitable  basis  of  comparison  with  analytical 
model  results,  the  load-deformation  relationships  for  each  test  specimen 
have  been  plotted  for  the  complete  loading  range.  Beam  response  in 
the  pre-cracked  state  is  elastic  to  within  a  reasonable  degree  of 
accuracy,  and  the  respective  torsion  and  bending  stiffness  constants 
for  the  seven  beams  are  given  in  Table  5.2.  Beyond  cracking,  beam 
behaviour  is  highly  inelastic  as  exhibited  in  Figures  5.1  to  5.5.  The 
torque  versus  rotation  curve  for  beam  R4  is  not  presented  as  it  closely 
retraces  the  appropriate  segment  of  the  curve  for  beam  R3.  In  Figures 
5.1,  5.2,  and  5.3,  the  abscissa  is  the  differential  rotation  over  the 
central  30  inch  length  of  the  beam.  The  deflection  ordinate  of  the 
moment-deflection  curves  is  the  vertical  displacement  of  the  central 
cross-section  due  to  bending  action  only.  Thus,  in  all  tests  where 
torque  was  present,  the  initial  deflection  measurements  had  to  be 
appropriately  adjusted  to  reflect  pure  bending  behaviour. 

In  monitoring  reinforcement  stress  levels,  strain  gages  were 
attached  only  to  the  longitudinal  and  transverse  conventional  rein- 
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forcement  at  the  beam's  central  cross-section.  In  most  cases,  the 
bottom  and  top  longitudinal  reinforcement  stress  values  are  averages 
of  three  readings,  whereas  the  hoop  stress  levels  are  derived  from  single 
gage  measurements.  In  the  test  setup,  the  beams  were  aligned  longitu¬ 
dinally  close  to  a  north-south  direction.  Thus,  differentiation  of  the 
beam's  east  and  west  sides  in  elevation  was  established.  The  respective 
stress  plots  are  illustrated  in  Figures  5.6  to  5.12,  with  identifying 
curve  nomenclature  given  in  Table  5.3.  All  stress  levels  are  tensile 
with  the  exception  of  those  of  the  top  longitudinal  bars. 

5.3  Observed  Behaviour 

In  almost  all  cases,  beams  were  not  tested  to  complete 
destruction,  but  were  loaded  just  beyond  their  observed  maximum  load 
carrying  capacity.  Upon  approaching  the  load  capacity  of  all  beams, 
primary  cracks  widened  markedly,  especially  in  those  tests  where 
there  was  a  high  ratio  of  bending  moment  to  torque.  Thus,  failure 
appeared  to  be  precipitated  by  excessive  yielding  of  conventional 
reinforcement,  yielding  being  most  pronounced  in  the  bottom  tension 
flange.  In  the  testing  of  beam  R1  to  destruction,  failure  occurred 
upon  the  crushing  of  a  compression  diagonal  across  the  width  of  the 
top  flange.  As  indicated  by  gaping  crack  widths,  the  strands  had 
undergone  considerable  inelastic  strain,  but  complete  disintegration 
of  the  beam  was  prevented  by  the  ductility  of  the  bottom  prestress 
strands.  Close  to  the  failure  of  beam  R4,  the  test  in  which  the  highest 
ratio  of  bending  moment  to  torque  was  applied,  a  similar  mode  of  failure 
was  observed,  but  well  defined  splitting  cracks  at  the  level  of  the 
bottom  longitudinal  reinforcement  were  also  apparent,  as  shown  in  Plate 
5.1.  The  only  beam  failure  that  differed  from  this  general  mode  of 
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failure  was  that  of  beam  Tl.  Excessive  local  crushing  and  a  punching 
shear  effect  became  apparent  beneath  one  torsion  arm  as  displayed  in 

Plate  5.2.  Initiation  of  the  local  failure  was  primarily  due  to  the 
absence  of  the  plaster  of  paris  pad  between  beam  and  torsion  arm,  thus 
resulting  in  extremely  severe  load  stress  concentrations.  However,  the 
load  at  failure  was  close  to  the  predicted  level. 

Representative  cracking  patterns  for  two  of  the  beams  are 
illustrated  in  Fig.  5.13.  Beams  R4  and  R5  were  selected  as  they 
represented  the  extremes  of  highest  and  lowest  ratio  of  bending  moment 
to  torque  respectively.  The  numbering  adjacent  to  the  crack  paths  in 
Fig.  5.13  corresponds  to  the  torsion  load  at  which  the  particular  crack 
was  formed.  The  suffix  B  designates  a  bending  load.  In  the  testing 
of  beam  R4,  the  ratio  of  bending  moment  to  torque  at  failure  was  close 
to  five,  and  thus  the  resultant  cracking  pattern  is  much  as  expected. 

At  the  other  end  of  the  testing  spectrum,  the  ratio  of  bending  moment 
to  torque  for  beam  R5  was  slightly  greater  than  one.  Consequently,  the 
formation  of  parallel  compression  diagonals  is  well  defined.  Generally, 
most  primary  cracks  were  7  to  8  inches  apart,  with  secondary  cracks  more 
closely  spaced  at  2  to  4  inches. 

5.4  Potential  Sources  of  Anomalies 


In  the  casting  of  test  specimens,  dimensional  tolerances  are 
inevitably  introduced,  and  their  severity  must  be  accounted  for  in  the 
estimation  of  accuracy  of  strength  and  deformation  predictions.  The 
most  significant  source  of  potential  dimensional  inaccuracy  unique  to 
the  seven  beams  cast  was  the  presence  of  the  styrofoam  voids.  Any 
substantial  movement  of  the  styrofoam  blocks  during  casting  would  have 
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introduced  dramatic  variation  in  the  thickness  of  the  thin  concrete 
webs.  Presence  of  the  closed  torsion  hoop  reinforcement 
was  used  to  advantage  to  secure  the  position  of  the  voids,  and  post¬ 
test  examination  revealed  that  little  variation  in  wall  thickness  was 
apparent. 


The  nature  of  the  beam  supports  was  such  that  a  uniform 
St.  Venant  torque  could  be  accommodated  along  any  length  of  the  beam. 
However,  the  presence  of  the  18  inch  long  solid  beam  ends  beyond  the 
supports  offered  longitudinal  warping  restraint  to  the  small  out-of¬ 
plane  warping  displacements  that  were  generated  at  the  beam  ends  during 
torsion  tests.  As  a  result,  the  beam  length  along  which  deformation 
measurements  were  taken  was  centrally  located  such  that  the  influence 
of  the  solid  ends  would  have  diminished  to  a  negligible  level. 

For  the  sake  of  expediency  in  ease  of  handling  and  positioning, 
the  hoops  and  bars  were  lightly  tack-welded  during  construction  of  the 
reinforcement  cages.  The  welding  was  sufficiently  light  such  that  no 
brittle  joints  would  be  formed.  The  rigidity  thus  introduced  into  the 
cages  had  no  effect  on  the  pre-cracked  beam  stiffness,  and  contributed 
little  to  the  post-cracked  stiffness.  Of  greater  significance  concerning 
the  reinforcement  cage  design  was  the  location  of  the  top  longitudinal 
bars  one  inch  in  from  the  hoop  corners.  This  feature  was  necessary  to 
facilitate  casting.  Although  Collins  and  Lampert25  state  that  the 
positions  of  the  corner  longitudinal  bars  define  cross-section  geometry 
of  the  cracked  concrete  beam  subjected  to  torsion,  the  movement  of  the 
top  corner  bars  did  not  have  a  marked  effect  on  the  torsion  failure 
loads  for  the  seven  beams  tested.  This  was  primarily  due  to  the  ratio 
of  bending  moment  to  torque  remaining  sufficiently  high  such  that 
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the  uncracked  state  of  the  top  flange  was  preserved  up  to  failure  in 
most  instances. 

The  presence  of  the  clamping  bolts  that  secured  the  torsion 
arms,  Amsler  point  loading  and  beam  support  apparatus  did  restrain  the 
propagation  of  cracking  beyond  the  central  test  section  in  each  beam. 

Individual  bolts  were  3/4"  in  diameter  and  effectively  acted  as  oversized 
stirrups . 

In  the  lower  range  of  elastic  behaviour  where  beam  deformations 
were  small,  the  linear  displacement  transducers  did  not  yield  consistently 
accurate  results  in  the  calculation  of  beam  twist.  This  inaccuracy 
arose  partly  through  the  method  of  attaching  to  the  concrete  the 
smooth  metal  plates  on  which  the  transducer  needles  impinged,  and 
partly  because  the  measuring  instrument  was  being  used  at  the  lower  limit 
of  its  range  of  accuracy. 
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© - © 

• 

Experimental  Bottom  Longitudinal  Steel 

Model  Bottom  Longitudinal  Steel 

A - A 

A 

Experimental*  Top  Longitudinal  Steel 

Model*  Top  Longitudinal  Steel 

b - 0 

■ 

Experimental  West  Hoop  Leg 

Model  West  Hoop  Leg 

O - €> 

♦ 

Experimental  East  Hoop  Leg 

Model  East  Hoop  Leg 

V - V 

▼ 

Experimental  Bottom  Hoop  Leg 

Model  Bottom  Hoop  Leg 

*  Note:  In  Figures  5.6  and  5.12,  the  designation  for  the 
ordinates  of  the  experimental  top  longitudinal  steel 
curves  is  the  symbol  A 


TABLE  5.3  NOMENCLATURE  FOR  CONVENTIONAL  REINFORCEMENT 
STRESS-LOAD  CURVES 
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FIG.  5.3  TORQUE-ROTATION  CURVES  FOR  BEAMS  Tl  AND  T2 


96. 


(sdj>|  ipui)  J.N3IAI0IAI  9NIQN38 


97. 


(sdp|  qouj)  1N3IMOW  9NIQN38 


STRESS  (psi) 


98. 


FIG.  5.6  STRESS-LOAD  CURVES  FOR  BEAM  R1  REINFORCEMENT 


FAILURE  LOAD  *  21.83  k 


54.7  ksi 
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FIG.  5.7  STRESS-LOAD  CURVES  FOR  BEAM  R2  REINFORCEMENT 
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TOTAL  LOAD  (kips) 

FIG.  5.8  STRESS-LOAD  CURVES  FOR  BEAM  R3  REINFORCEMENT 
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FIG.  5.9  STRESS-LOAD  CURVES  FOR  BEAM  R4  REINFORCEMENT 
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FIG.  5.10  STRESS-LOAD  CURVES  FOR  BEAM  R5  REINFORCEMENT 
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ORIGIN  LOAD  (kips) 


FIG.  5.11  STRESS-LOAD  CURVES  FOR  BEAM  T1  REINFORCEMENT 
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FIG.  5.12  STRESS-LOAD  CURVES  FOR  BEAM  T2  REINFORCEMENT 
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FIG.  5.13  CRACKING  PATTERNS  FOR  BEAMS  R4  AND  R5 
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PLATE  5.1  FAILURE  MODE  OF  BEAM  R4 


PLATE  5.2  LOCAL  FAILURE  OF  BEAM  T1 


CHAPTER  VI 


EVALUATION  OF  COMPUTER  MODEL  RESULTS 


6.1  Introduction 


Motivation  to  develop  this  finite  element  computer  model  arose 
through  the  need  for  a  flexible,  precise  analytical  method  of  analysis 
of  concrete  box  girders  subjected  to  a  general  load  condition.  The 
principles  of  the  mechanics  incorporated  in  the  computer  model  have  been 
described  in  detail  in  Chapter  3  but  performance  of  the  assembled  model 
was  not  addressed.  This  chapter  focuses  on  the  evaluation  of  the  computer 
model  results  through  comparison  with  experimental  tests  and  related 
current  theory. 

Initially ,  those  aspects  that  strongly  influence  the  computer 
model  response  and  do  not  represent  common  finite  element  modelling 
considerations,  are  examined  to  establish  an  insight  into  model 
behavioural  characteristics  revealed  in  the  subsequent  comparisons. 

The  subject  of  comparison  is  the  spectrum  in  behaviour  of  the  seven 
prestressed  concrete  box  girders  that  were  tested  in  the  experimental 
program  described  in  Chapter  4.  Following  the  presentation  of  the 
experimental  test  results  and  the  corresponding  computer  model  results 
in  Section  6.3,  theoretical  estimates  of  ultimate  beam  strength  and 
interactive  response  are  in  turn  compared  with  the  analytical  model 
predictions  in  Section  6.4.  Performance  of  the  finite  element  model 
is  assessed  on  the  basis  of  the  two  comparative  presentations,  and  is 
pursued  in  Section  6.5. 
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6*2  Extraordinary  Influences  on  Computer  Model  Response 

In  assessing  the  performance  of  the  computer  model,  considera¬ 
tion  must  be  given  to  the  following  important  aspects  that  both  individu¬ 
ally  and  collectively  have  a  strong  bearing  on  the  model  results. 

Several  aspects  are  unique  in  the  application  of  the  analytical  model 
to  simulate  the  response  of  the  seven  prestressed  concrete  box  beams 
tested  in  the  experimental  program.  All  aspects  are  of  either  a  material 
or  structural  nature. 

6.2.1  Material  Behaviour  Aspects 

In  simulating  the  stiffness  of  an  uncracked  reinforced 
concrete  beam,  the  material  property  of  paramount  concern  is  the  initial 
elastic  modulus  for  the  concrete.  Since  the  stiffness  contribution  of 
the  reinforcement  in  an  underreinforced  concrete  beam  is  small  compared 
with  that  of  the  concrete,  an  inaccurate  determination  of  the  initial 
concrete  modulus  will  result  in  a  correspondingly  inaccurate  prediction 
of  beam  deformations.  Since  concrete  cyclinder  tests  did  not  yield 
consistent  results  and  empirical  formulae  were  not  sufficiently  accurate, 
the  initial  deflection  measurements  of  the  test  beams  were  used  in 
precisely  calculating  the  respective  initial  concrete  moduli. 

A  significant  discrepancy  was  observed  between  the  average 
crack  spacing  value  derived  using  Eq.  3.29  and  the  corresponding 
experimental  value  determined  from  the  beam  cracking  patterns  at  failure. 
Since  aggregate  interlock  stiffness  is  a  function  of  crack  width  and 
thus  of  average  crack  spacing,  adoption  of  the  considerably  smaller 
theoretical  value  would  have  resulted  in  an  overestimation  of  post¬ 
cracking  concrete  stiffness. 
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Reconstitution  of  the  stiffness  of  a  reinforced  concrete 
element  beyond  cracking  is  difficult  in  that  the  contribution  of  "dowel 
action"  cannot  be  determined  accurately.  Research  has  not  conclusively 
established  the  quantitative  effect  upon  dowel  stiffness  of  several 
layers  of  web  reinforcement  or  the  inclination  of  the  crack  to  the 
reinforcement  axis  other  than  at  ninety  degrees.  However,  consequences 
of  the  lack  of  complete  understanding  and  definition  of  this  phenomenon 
are  not  of  major  significance  since  dowel  action  accounts  for  less  than 
20  percent  of  the  shear  strength  of  a  cracked  reinforced  concrete  beam. 

6.2.2  Structural  Aspects 

Although  both  conventional  and  prestress  longitudinal 
reinforcement  are  represented  in  the  analytical  model  by  one  dimensional 
bar  finite  elements,  the  same  method  cannot  be  extended  to  the  modelling 
of  stirrup  and  hoop  reinforcement  in  this  computer  model.  This 
restriction  has  arisen  through  characterizing  the  concrete  finite 
element  behaviour  by  the  stress-strain  condition  at  the  centroid  of 
the  element.  When  the  principal  tensile  stress  at  the  concrete  element 
centroid  exceeds  the  specified  tensile  strength  of  concrete,  the  element 
is  designated  as  a  cracked  element.  However,  the  element  is  not  fractured 
by  several  parallel  cracks  distributed  across  the  element,  but  by  one 
crack  passing  through  the  centroid.  Thus,  when  the  vertical  web  rein¬ 
forcement  is  concentrated  into  bars  located  at  the  vertical  boundaries 
of  an  element  whose  aspect  ratio  exceeds  unity,  it  is  highly  probable 
t-hat  the  sole  centroidal  crack  will  not  intersect  the  side  boundaries 
of  the  element.  Consequently,  the  cracked  element  cannot  utilize  its 
web  reinforcement,  and  a  premature  shear  failure  will  occur.  To  resolve 
this  modelling  difficulty,  the  web  reinforcement  is  distributed  throughout 
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the  element  as  a  uni-directional  steel  mesh  of  closely  spaced  bars. 

When  cracking  occurs,  the  complete  mesh  stiffness  for  the  concrete 
element  is  engaged.  Thus,  modelling  of  the  mechanics  of  stirrup  and 
hoop  reinforcement  in  the  cracked  concrete  condition  closely  reflects 
real  behaviour.  Figures  6.1(a)  and  6.1(b)  illustrate  the  problem  and 
solution  of  this  vitally  important  aspect  of  analytical  modelling. 

In  Section  3.2.1,  an  expression  was  derived  for  the  variation 
in  shear  stress  across  a  "thick"  wall  as  a  function  of  the  uniform 
shear  flow,  in  the  form: 

_  Ac 

V  A'  (3.12) 

Clearly,  for  those  tubular  members  whpse  height  and  width  dimensions 
are  only  moderately  large  compared  with  the  wall  thickness,  the  maximum 
shear  stress  at  the  wall  surface  can  be  considerably  larger  than  that 
at  the  wall  mid- thickness;  ie.  the  uniform  shear  flow.  Considering  that 
plane  stress  finite  elements  do  not  permit  shear  stress  variation 
across  their  thickness,  and  that  such  elements  are  located  at  the  plates' 
mid-thickness,  only  uniform  shear  flows  can  be  adequately  represented. 

The  computer  model's  insensitivity  to  shear  stress  variation  is  of 
concern  since  surface  cracking  in  concrete  will  immediately  propagate 
through  the  entire  thickness  of  a  tubular  member's  wall.  Thus,  if  the 
torsional  shear  stresses  are  high  at  the  critical  cross-section  of  a 
concrete  box  girder,  the  computer  model  will  yield  under conservative 
estimates  of  cracking  load  and  to  a  lesser  extent  post-cracking  member 
stiffness . 


A  further  possible  consequence  of  the  location  of  plane  stress 
elements  at  wall  mid-thickness  is  the  inaccurate  prediction  of 
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the  ultimate  bending  moment  capacity  of  an  underreinforced  concrete 
box  girder.  In  the  analytical  model,  the  lever  arm  length  from  the 
force  resultant  in  the  concrete  compression  flange  to  the  tension 
reinforcement  at  ultimate  load  conditions  is  fixed  by  geometry.  If  the 
compression  flange  is  thick  with  respect  to  the  girder’s  depth,  the 
concrete  compressive  force  resultant  will  be  located  between  the  wall 
mid-thickness  and  the  upper  flange  surface,  the  height  of  the  resultant 
force  above  the  compression  flange  mid-thickness  constituting  the 
difference  between  the  experimental  and  model  lever  arms.  Also,  since 
the  thick  concrete  compression  flange  of  the  computer  model  is  too  stiff 
close  to  failure,  the  tension  reinforcement  will  be  forced  to  carry  a 
disproportionately  higher  load,  further  reducing  the  ultimate  bending 
moment  capacity  of  the  girder  below  its  actual  strength. 

Before  a  reinforced  concrete  beam  is  tested,  the  concrete 
and  reinforcement  are  already  stressed  through  shrinkage  and  creep 
member  deformations,  the  effects  of  these  phenomena  being  most  pronounced 
in  precast  prestressed  concrete  beams.  Failure  to  include  the  presence 
of  the  tensile  concrete  shrinkage  stresses  can  result  in  significant 
overestimation  of  the  cracking  load,  and  a  slight  distortion  of  pre¬ 
cracked  deformations. 

Although  the  longitudinal  warping  restraint  of  beam  end 
diaphragms  is  not  normally  of  real  structural  significance,  the  modelling 
of  the  seven  prestressed  box  beams  tested  in  the  experimental  program 
was  complicated  by  the  presence  of  the  18  inch  long  solid  beam  ends 
extending  beyond  each  beam  support.  Since  the  warping  resistance  of 
the  solid  ends  was  comparable  to  that  of  the  box  beams,  accurate 
representation  of  their  influence  on  beam  behaviour  was  essential. 
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Consequently,  the  approach  outlined  in  Section  3. 4. 2. 7  was  adopted, 
and  the  warping  stiffnesses  established  in  the  auxiliary  finite  element 
program  were  read  in  as  input  in  the  principal  analytical  model. 

As  shown  in  Figures  4.7,  4.8,  and  4.9,  clamping  bolts  were 
used  to  maintain  the  position  of  testing  equipment  at  several  locations 
along  the  length  of  the  test  specimens.  Verified  by  observations  during 
testing,  the  clamping  bolts  and  the  radial  torsion  loading  arms  effectively 
acted  as  very  stiff  stirrups,  greatly  restricting  the  propagation  of 
web  cracking  in  their  vicinity.  Thus,  their  presence  in  the  analytical 
model  is  essential  from  both  deformation  and  strength  aspects.  In  a 
similar  manner  to  that  adopted  for  conventional  stirrup  and  hoop 
reinforcement,  equivalent  steel  meshes  were  included  in  adjacent  concrete 
elements  to  represent  their  influence. 


In  each  of  the  seven  test  beams,  the  solid  cross-section 
infringed  upon  the  test  span,  the  distance  from  the  support  to  the 
commencement  of  the  styrofoam  voids  varying  from  7.5  to  13.5  inches 
(Fig.  4.4).  Since  the  stress  levels  at  the  beam  ends  were  not  sufficient 
to  produce  cracking,  the  primary  modelling  concern  was  deformational 
accuracy,  especially  with  respect  to  bending  deflections.  Deformation 
characteristics  were  modelled  closely  through  the  use  of  plane  stress 

finite  elements  of  an  appropriate,  uniform  thickness  such  that  the  moment 
of  inertia  remained  unchanged. 

In  the  casting  of  prestressed  concrete  beams,  small  flexural 
cracks  will  often  appear  in  the  top  flange  after  transfer.  Upon  the 
application  of  positive  moment,  these  cracks  will  immediately  close, 
and  the  subsequent  beam  response  will  not  reflect  their  initial 


presence. 


113. 


Provision  is  made  within  the  model  to  reproduce  such  behaviour.  If  a 
model  crack  closes  in  the  load  increments  immediately  following  transfer, 
the  full  concrete  constitutive  matrix  is  recovered.  Should  a  crack  close 
at  a  later  stage  in  the  loading  sequence,  the  aggregate  interlock 

stiffness  does  not  increase,  but  is  maintained  at  the  level  prior  to 
crack  closing. 

Comparison  of  Computer  Model  and  Experimental  Results 

In  this  Section,  the  computer  model  and  experimental  results 
for  the  seven  prestressed  box  girders  are  presented  in  a  form  that 
facilitates  a  thorough  comparison.  Assessment  of  the  analytical  model 
performance  on  the  basis  of  this  comparison  is  treated  in  Section  6.5.1. 

Since  the  complete  spectrum  of  beam  response  from  initial 
loading  to  failure  is  to  be  examined,  plots  of  torque  versus  rotation, 
bending  moment  versus  central  beam  deflection,  and  reinforcement  stresses 
versus  load  have  been  prepared  for  every  test  beam,  each  plot  displaying 
the  corresponding  model  and  experimental  coordinates.  Monitoring  of 
reinforcement  stresses  embraces  the  behaviour  of  the  top  and  bottom 
flange  conventional  reinforcement  at  the  beam  span  centreline  and  the 
four  legs  of  the  centrally-located  hoop.  The  legend  for  identification 
of  the  reinforcement  stress  plots  is  provided  in  Table  5.3. 

The  twenty  figures  displaying  comparative  beam  plots,  Figures 
6.2  to  6.21,  are  not  completely  comprehensive  as  two  plots  have  not  been 
presented.  Reference  to  the  Figure  List  indicates  that  the  torque 
versus  rotation  curve  for  beam  R4  and  the  hoop  reinforcement  stresses 
versus  load  graphs  for  beam  T2  are  the  two  absent  plots.  For  beam 


f 


114. 


M‘  ltS  tor<3ue-rotation  curve  closely  follows  the  respective  plot  for 
beam  R3  since  the  loading  of  both  beams  only  differed  significantly 
beyond  the  scope  of  the  torsion  plot  of  beam  R4.  During  the  test  of 
beam  T2,  the  central  hoop  strain  gages  exhibited  erratic  readings 
which  prejudiced  their  experimental  value.  In  the  torque-rotation  graph 
of  the  experimental  results  of  beam  T2,  shown  in  Fig.  6.20,  the  initial 
coordinates  are  not  plotted  as  they  displayed  considerable  inconsistency. 
The  two  computer  model  curves  for  the  bending  moment-deflection  response 
of  beam  R4  represent  the  two  loading  extremes  of  the  heavy  Amsler  jack 
load  applied  at  the  beam  centreline.  If  the  beam  failure  mechanism 
comprised  a  centreline  hinge  in  the  compression  flange,  and  the  loading 
Plate  was  sufficiently  stiff  to  bridge  the  hinge,  the  loading  pattern 
of  Fig.  6.22(b)  is  feasible,  in  contrast  to  that  illustrated  in  Fig. 
6.22(a)  where  complete  contact  is  assumed. 

Two  computer  model  graph  characteristics  exhibited  in  the 
majority  of  bending  moment-deflection  and  torque-rotation  plots  require 
clarification.  For  those  beams  where  failure  of  the  finite  element 
model  occurred  within  a  load  increment,  the  model  strength  results  are 
described  by  a  band  delineating  the  lower  and  upper  load  bounds.  The 
second  point  of  clarification  is  the  significance  of  the  "translated" 
computer  model  curve  shown  in  those  figures  where  there  is  considerable 
discrepancy  between  experimental  and  model  cracking  loads.  If  the 
computer  model  behaviour  was  adjusted  such  that  its  cracking  load  agreed 
closely  with  the  corresponding  experimental  load,  the  projected  analytical 
model  results  are  represented  by  the  "translated"  computer  model  curve. 

Most  importantly,  the  translation  procedure  does  not  entail  pre-cracking 

or  post-cracking  stiffness  modification.  Development  of  the  translated 
curve  is  detailed  in  Section  6.5.1. 
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The  cracking  and  ultimate  strength  combined  loading  conditions 
for  the  experimental  and  the  corresponding  computer  model  beams  are 
presented  in  Table  6.1.  Invariably,  the  deformation  value  at  the  point 
of  ultimate  strength  was  indeterminable  from  both  the  experimental  and 
computer  model  results.  Consequently,  there  is  no  objective  method  of 
comparing  experimental  and  computer  model  postcracking  stiffnesses, 

and  thus  a  subjective  evaluation  of  the  corresponding  curves  is  the 
best  recourse. 

6*4  Comparison  of  Computer  Model  Results  with  Current  Theory 

Of  the  vast  number  of  research  publications  in  the  field  of 
torsion,  bending,  and  shear  in  reinforced  concrete,  reference  has  been 
made  to  Collins  and  Lampert25  in  estimating  pure  ultimate  torsional 
strength,  Thurliman50  in  the  evaluation  of  bending  moment-shear  interaction 
and  pure  shear  strength,  and  Elfgren51  concerning  torsion-bending  and 
torsion-bending-shear  interaction.  Since  current  theory  cannot  predict 
post-cracking  member  deformations  under  combined  loading,  only  ultimate 
strength  predictions  will  be  considered.  Evaluation  of  cracking  loads 
will  not  be  treated  as  this  area  of  research  has  been  exhaustively 
examined  in  the  past,  and  resolved  to  a  satisfactory  degree.  Similar 
to  the  previous  Section,  this  Section  will  only  present  a  comparison 
of  theoretical  and  model  predictions,  with  the  assessment  and  implications 
of  the  comparison  addressed  in  Section  6.5.2.  Since  geometry,  concrete 
strength,  and  reinforcement  levels  are  almost  identical  within  the  two 
test  beam  classifications,  average  ultimate  strengths  will  be  presented 
for  the  rectangular  and  trapezoidal  beam  categories. 
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6*4.1  Ultimate  Strength 

In  the  application  of  the  Space  Truss  Theory  proposed  by 
Collins  and  Lampert25  to  estimate  the  ultimate  torsional  capacity  of  a 
prestressed  concrete  box  girder,  the  form  of  their  equations  is  not 
modified  to  account  for  initial  prestress.  The  prestress  strand  is 
simply  considered  as  reinforcement  of  strength  f*  A*,  where  f*  is 

S  S  s 

the  yield  stress  of  the  prestress  steel  and  A*  the  area  of  the  individual 
strand.  For  both  the  rectangular  and  trapezoidal  beams,  the  weaker 
top  reinforcement  determines  the  ultimate  torsional  capacity  of  the 
two  beam  classifications.  When  considering  the  trapezoidal  beams,  the 
web  prestress  strand  was  distributed  to  the  top  and  bottom  flange 
stringers  such  that  the  ultimate  bending  moment  capacity  remained 
unchanged.  The  ultimate  torsional  capacity  is  given  by 

£Z  S 

To  =  2Ao/  — ^  •  -Z  „ 

u  t  (6.1) 

where  Ao  -  area  enclosed  by  corner  longitudinal  bars,  EZ  =■  twice  the 
sum  of  yield  forces  of  longitudinal  bars  in  the  weaker  flange,  S  =  hoop 
yield  force,  u  =  corner  longitudinal  bar  perimeter,  and  t  =  hoop 
spacing.  The  average  ultimate  torsional  capacity  of  the  five  rectangular 
beams  using  the  Space  Truss  Theory  is  379  inch  kips,  and  the  corresponding 
strength  for  the  two  trapezoidal  beams  is  239  inch  kips. 

Since  the  Space  Truss  Theory  cannot  accurately  estimate  the 
compression  zone  depth,  and  Skew  Bending  Theory  cannot  accommodate  the 
presence  of  longitudinal  web  reinforcement  or  longitudinal  reinforcement 
of  different  yield  strengths,  classical  bending  theory52  as  currently 
incorporated  in  prestressed  concrete  design  has  been  used  to  evaluate 
the  ultimate  bending  moment  capacity  of  the  two  beam  types.  For  the 
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five  rectangular  beams,  the  average  ultimate  bending  moment  capacity  is 
1515  inch  kips,  and  985  inch  kips  for  the  two  trapezoidal  beams. 


;  Calculation  of  the  pure  shear  strength  of  the  two  box  beam 
types  follows  the  approach  advocated  by  Thurliman50.  On  the  basis  of 


a  simplified,  generalized  space  truss  model,  Thurliman  devised  the 

following  expression  for  the  "plastic  shear  force"  V  • 

po* 


V 

po 


/  2  F  •  S  •  - 
yi  y  t 


(6.2) 


where  Fyl  =  yield  force  of  bottom  flange  longitudinal  stringers,  Sy  = 
yield  force  of  stirrups  at  cross-section  considered,  h  =  ultimate  bending 
moment  lever  arm,  or,  in  the  absence  of  an  accurate  estimation  of  this 
value,  the  centre  to  centre  distance  between  the  top  and  bottom  stringers, 
and  t  =  stirrup  spacing. 


However,  since  both  the  transverse  and  longitudinal  reinforce¬ 
ment  must  yield  simultaneously  before  failure  for  the  above  derivation 
to  be  valid,  the  variable  angle  of  inclination  of  the  compression  concrete 
diagonals  must  lie  within  the  range: 

•  5  1  tan  ag  <  2  (6.3) 

where  ag  =  angle  of  inclination  of  compression  struts  to  horizontal. 

For  such  an  inclination  range,  the  maximum  shear  resistance 
of  a  concrete  beam  is  modified  to 

Vpmax  /  k  Vpo  (6.4) 

where  Vpmax  =  maximum  shear  resistance,  and 
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In  the  application  of  the  above  equations  to  prestressed 
concrete  (Thurliman’s  equations  were  developed  for  reinforced  concrete), 
the  axial  force  N  in  the  equation  below  is  not  zero,  and  corresponds  to 
the  initial  prestress  force. 

v  _  M  N  ,  V 

Fyl  "  h  +  2  +  2  '  cot0ls  (6.6) 

The  form  of  the  developed  equations  remains  unaltered  if  the  term  F 

yl 

is  replaced  by  the  term  (F  ^  4  ^— )  where  N  is  the  initial  prestress 

force,  and  h’  is  the  distance  from  the  prestress  strand  to  the  concrete 
force  resultant.  Mostly,  h'  equals  h.  Therefore,  Eq.  6.6  is  valid  in 
its  original  form  if  F^  is  defined  as  the  total  bottom  stringer  yield 
strength  including  the  initial  prestress  force. 

For  the  rectangular  box  beams,  the  plastic  shear  force  is 
70.7  kips  and  the  maximum  shear  resistance  is  49.5  kips.  Similarly, 
for  the  trapezoidal  box  beams,  the  plastic  shear  force  is  57.4  kips 
and  the  maximum  shear  resistance  is  40.5  kips. 

Since  a  computer  model  loading  combination  could  not  be 
devised  that  produced  a  pure  shear  loading  similar  to  that  considered  in 
the  preceding  theory,  an  analytical  estimate  of  the  ultimate  shear 
capacity  could  not  be  determined. 

The  experimental  and  corresponding  computer  model  ultimate 
strengths  are  summarized  in  Table  6.3. 
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6.4.2  Torsion  -  Bending  Interaction 

Of  the  seven  box  beams  tested  in  the  experimental  program, 
five  were  subjected  to  torque  and  bending  moment  only,  those  beams  being 
Rl,  R2 ,  R5,  Tl,  and  T2. 


The  interaction  equations  of  Elfgren51,  in  an  identical  form 
to  those  developed  by  Collins  and  Lampert25,  are  utilized  to  describe 
beam  response  under  combined  torque-bending  moment  load  conditions. 

Under  such  combined  loading,  beam  failure  under  the  action  of  a  positive 
bending  moment  can  be  either  of  two  modes,  mode  t  when  failure  is  initiated 
by  yielding  of  the  bottom  flange  longitudinal  reinforcement,  or  mode  c 
when  yielding  of  the  top  flange  stringers  precipitates  beam  failure. 


Mode  t 


2L  +  (i)2r 

M  'T  ' 
o  o 


=  1 


Mode  c 
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where  r  =  ratio  of  yield  forces  of  top  and  bottom  flange  longitudinal 
reinforcement . 


In  determining  the  value  of  term  r  in  the  above  two  equations, 
the  full  yield  strength  of  the  prestress  strand  is  used.  For  the  two 
trapezoidal  beams  that  contain  web  prestress  strand  reinforcement, 
distribution  of  the  longitudinal  web  steel  to  the  top  and  bottom  flange 
stringers  is  such  that  the  ultimate  bending  moment  capacity  of  the  beams 
remains  unchanged.  The  yield  force  ratio  "r"  is  .2377  and  .448  for 
the  rectangular  and  trapezoidal  beam  categories  respectively. 


As  a  basis  of  comparison,  the  interaction  Equations  6.7  and 
6.8  are  plotted  in  Figure  6.23  for  both  the  rectangular  and  trapezoidal 
beam  types,  together  with  the  corresponding  computer  model  results  for 
the  five  beams  that  failed  under  torque  and  bending  moment  loading  only. 
The  computer  model  results  are  displayed  in  their  dimensional  and  non- 
dimensional  form  in  Table  6.3,  all  results  being  average  values  for 
the  load  increment  in  which  failure  of  the  analytical  model  occurred. 

6*4.3  Torsion  -  Bending  -  Shear  Interaction 

Only  two  of  the  seven  test  beams  were  subjected  to  shear  in 
addition  to  torque  and  bending  moment  at  their  critical  cross-sections, 
the  two  beams  being  rectangular  box  beams  R3  and  R4. 

The  interaction  equations  adopted  to  define  beam  response 
under  the  combined  loading  of  torque,  bending  moment,  and  shear  are 
those  proposed  by  Elfgren52.  In  conjunction  with  the  two  previously 
defined  modes  of  failure  the  presence  of  shear  introduces  an  additional 
mode  of  failure,  mode  s,  where  the  compression  zone  is  formed  on  the 
side  of  the  beam.  The  corresponding  interaction  equations  are  as  follows: 

Mode  t 


(6.9) 


Mode  c 
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Mode  s 

2  _2r  .  V  2  2r  ,  TV  2r  2 

To  r+1  r+1  ToVo  *  r+1  *  j/  l  +  b'/h'  "  1  (6.11) 

where  b'  =  horizontal  centre  to  centre  distance  of  corner  stringers 
in  top  or  bottom  flanges,  and  h'  =  vertical  centre  to  centre  distance 
of  corner  stringers  in  top  and  bottom  flanges. 


In  conjunction  with  the  above  equations,  Thurliman's  inter¬ 
action  equations  for  bending  moment  and  shear  is  also  considered  as  an 
additional  interactive  constraint  that  must  not  be  violated: 


M  V 

—E_  +  /_£_)  2  =  ± 

po  po 


(6.12) 


where  applied  moment,  and  -  "plastic  moment"  or  ultimate  bending 

moment  capacity.  From  earlier  discussion,  the  applied  shear  V  must 

P 

not  exceed  the  maximum  shear  capacity  V  established  by  Eq.  6.4. 

pUldiX 

Similarly,  as  a  consequence  of  the  inclination  of  the  compression  struts 


being  restricted  to  .5  £  tan  CXg  £  2,  the  maximum  applied  bending  moment 
must  satisfy  the  condition: 
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The  interaction  Eq.  6.12,  together  with  the  imposed  limits  of  Equations 
6.4  and  6.13,  is  shown  in  Fig.  6.24. 


Using  Eq .  6.2  to  evaluate  the  plastic  shear  force  V  appearing 
in  Equations  6.9,  6.10  and  6.11,  the  interaction  equations  for  the 


three  modes  of  failure  of  beams  R3  and  R4  are  illustrated  in  Fig.  6.25. 
Only  one  set  of  interaction  equations  is  shown  as  the  respective  equations 
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for  the  two  beams  are  very  similar.  All  three  interaction  equations  have 
been  simplified  through  the  evaluation  of  the  (~)  terms  which  are  then 
transferred  to  the  numerical  right  hand  sides  of°their  respective 
equations.  Thus,  three  dimensional  interaction  is  reduced  to  a  two 
dimensional  torque-bending  moment  interaction. 

In  Table  6.3  that  displays  the  computer  model  results  for 
beams  R3  and  R4,  the  bending  moment  failure  loads  have  been  adjusted 
to  allow  for  the  presence  of  the  central  downward  vertical  concentrated 
load  illustrated  in  Fig.  1.1(c)  in  Appendix  I.  Under  such  a  loading 
system,  Thurliman50  has  established  that  the  principal  design  cross- 
section  for  shear  is  not  located  beneath  the  central  load,  but  at  a 
distance  h  either  side  of  the  concentrated  load.  Elfgren's51  examination 
of  the  effect  of  a  concentrated  load  yields  a  similar  result.  Consequently, 
the  failure  bending  moments  for  beams  R3  and  R4  have  been  altered 

accordingly.  The  adjusted  computer  model  results  for  the  latter  two 
beams  are  plotted  in  Fig.  6.25. 

6.5  Assessment  of  Computer  Model  Results 

6,5,1  Computer  Model  Assessment  in  Light  of  Experimentation 

6* 5.1.1  Prominent  Aspects  of  Model’s  Performance:  In  the 
uncracked  state,  the  computer  model  yielded  an  accurate  assessment  of 
beam  stiffness  since  the  initial  modulus  of  elasticity  for  concrete 
was  evaluated  directly  from  beam  deformations  at  low  load  levels.  This 
procedure  in  determining  the  concrete  modulus  was  adopted  as  cylinder 
tests  conducted  for  the  sole  purpose  of  modulus  measurement  produced  a 
wide  scatter  of  results.  In  this  study  the  need  for  estimating  the 
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modulus  of  concrete  as  accurately  as  possible  is  paramount  as  both 
deformation  and  strength  characteristics  are  examined.  The  few  instances 
of  disparity  in  agreement  between  experimental  and  model  behaviour  in 
the  pre-cracked  condition  appear  in  the  torque-rotation  relationship 
since  the  experimental  deflection  monitoring  equipment  was  incapable  of 
consistently  accurate  measurements  of  very  small  differential  deflections 
of  adjacent  beam  locations. 

Although  corresponding  experimental  and  model  elastic  stiffnesses 
are  in  close  agreement,  a  very  discernible  discrepancy  is  evident  at 
the  onset  of  cracking.  For  all  beams  where  cracking  occurred  under  a 
combined  torque-bending  moment  load  condition:  i.e.  all  beams  with  the 
exception  of  R3  and  R4;  the  model  beam  cracked  at  a  considerably  higher 
load  than  the  experimental  beam.  The  three  potential  sources  of  this 
deviation  in  behaviour  are:  concrete  tensile  strength,  concrete  shrinkage 
stresses,  and  variation  in  shear  stresses  across  the  wall  thickness.  Of 
the  three  potential  sources,  the  effect  of  the  variation  in  wall  shear 
stresses  is  the  most  dominating  influence  on  behavioural  discrepancy. 

In  Section  3.2.1,  the  variation  in  the  St.  Venant  torsional  shear  stress 
from  the  wall  surface  to  wall  mid-thickness  was  shown  to  be  equal  to  the 
ratio  of  the  cross-sectional  area  to  the  area  enclosed  by  the  corner 
longitudinal  reinforcement  stringers.  For  the  rectangular  and  trapezoidal 
beams,  the  ratio  equals  .61  and  .87  respectively.  Since  the  analytical 
model  evaluates  the  torsional  shear  stresses  at  the  wall  mid-thickness, 
the  maximum  torsional  shear  stresses  at  the  surface  of  the  experimental 
beams  are  61%  and  87%  higher  than  the  corresponding  model  stresses  for 
the  rectangular  and  trapezoidal  beams.  This  substantial  discrepancy 
between  mid-thickness  and  surface  torsional  shear  stresses  is  critical  as 
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cracks  in  the  outer  surface  fibres  immediately  propagate  through  the 
entire  wall  thickness.  If  the  shear  stress  variation  were  taken  into 
account  in  the  calculation  of  the  model  concrete  principal  tensile 
stresses,  the  discrepancy  in  cracking  loads  would  be  substantially 
reduced.  However,  difficulty  in  accurately  measuring  the  concrete 
tensile  strength  and  uncertainty  in  establishing  concrete  shrinkage 
stresses  can  collectively  contribute  to  inaccurate  cracking  estimates. 
Thus,  although  the  variation  in  torsional  shear  stresses  can  be  taken 

into  account,  accurate  prediction  of  cracking  loads  may  not  be  achieved 
consistently. 


As  a  result  of  the  model's  higher  cracking  load,  a  response 
delay  is  exhibited  in  the  bending  moment-deflection,  torque-rotation, 
and  reinforcement  stresses-load  graphs.  In  those  moment -deformat ion 
figures  where  the  response  delay  is  considerable,  a  "translated  computer 
model  curve"  has  been  plotted  in  addition  to  the  actual  model  curve  to 
illustrate  the  projected  computer  model  response  if  the  cracking  load 
discrepancy  did  not  occur.  To  define  the  translated  curve,  the  post- 
cracking  portion  of  the  actual  model  curve  is  moved  horizontally  to  the 
right  until  the  tangent  projection  of  the  inelastic  curve  below  its 
local  origin  intersects  the  elastic  slope  at  the  experimental  cracking 
load.  The  prophetic  significance  of  the  translation  procedure  illu¬ 
strated  in  Fig.  6.1(b)  is  as  follows.  Should  a  concrete  box  beam  be 
modelled  where  no  account  is  taken  of  the  variation  of  torsional  shear 
stresses  and  the  presence  of  shrinkage  stresses,  the  analytical  load- 
deformation  results  will  be  of  the  form  of  curve  OBC  in  the  latter 
figure,  where  a  considerable  discrepancy  is  apparent  between  experimental 
and  model  cracking  loads.  However,  if  both  the  two  previously  specified 


125. 


sources  of  error  in  cracking  prediction  are  taken  into  account,  the 
modified  analytical  results  are  defined  by  the  curve  OAB'C'.  The 
stiffnesses  at  corresponding  intermediary  points  along  BC  and  B'C', 
points  D  and  D'  for  example,  correspond  closely  as  the  effects  of 
shrinkage  and  torsional  shear  stress  variation  are  minimal  in  comparison 
to  the  dramatic  redistribution  of  stresses  that  follows  cracking.  In 
beam  tests  that  exhibited  considerable  ductility  beyond  cracking,  the 
inelastic  slope  up  to  moderate  load  levels  exhibited  little  decay,  and 
thus  construction  of  the  segment  AB'  is  a  reasonable  approximation  of 
the  initial  modified  post-cracking  response.  The  principal  motivation 
m  introducing  the  translated  curves  is  that  a  closer  subjective  comparison 
can  be  made  of  experimental  and  model  post-cracking  stiffnesses,  and 
projected  member  deformations  are  predicted  more  accurately. 

With  few  exceptions  model  post-cracking  stiffness  closely 
follows  experimental  behaviour,  and  deformation  predictions  of  the 
translated  model  curves  are  consistently  accurate  before  entering  the 
highly  inelastic  deformation  regions  close  to  ultimate  failure.  The 
less  pronounced  effect  of  the  torsional  shear  stress  variation  after 
beam  cracking  arises  since  redistribution  of  forces  through  concrete 
cracking  and  reinforcement  yielding  produces  large  stress  variations 

that  greatly  exceed  the  influence  of  surface  to  mid-thickness  shear 
stress  variation. 


At  ultimate  failure  load  conditions,  the  computer  model  bending 
moments  are  consistently  8%  to  10%  below  the  corresponding  experimental 
moments  for  those  beams  that  were  subjected  to  a  high  ratio  of  bending 
moment  to  torque.  This  modelling  inaccuracy  is  largely  due  to  the  method 
of  delineation  of  the  beam  cross-section  in  the  finite  element  model. 
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In  representing  a  concrete  wall  by  a  two-dimensional  plane  stress  finite 
element  located  at  the  wall  mid-thickness,  the  moment  of  inertia  of  the 
uncracked  beam  is  accurately  represented  and  actual  reinforcement 
locations  can  usually  be  accommodated.  However,  the  thickness  of  the 
concrete  compression  zone  and  the  location  of  the  compressive  force 
resultant  are  fixed,  effectively  predefining  the  model  bending  moment 
lever  arm  independently  of  reinforcement  levels  and  critical  beam  cross- 
sectional  parameters.  For  both  the  rectangular  and  trapezoidal  finite 
element  meshes,  the  bending  moment  lever  arms  are  4.5%  shorter  than 
theoretical  estimates  at  ultimate  load  conditions,  resulting  in  almost  a 
5%  reduction  in  the  ultimate  bending  moment  capacity.  The  higher  stress 
levels  in  the  longitudinal  tension  reinforcement  of  the  computer  model 
have  little  effect  on  post-cracking  stiffness  except  in  the  region  close 

to  failure  where  model  stiffness  is  invariably  less  than  experimentation 
ind icated. 


Several  less  significant  sources  of  error  are  introduced  through 
selection  of  the  geometry  of  the  finite  element  mesh.  Defining  cross- 
sectional  geometry  by  mid-thickness  dimensions  dictates  the  location  of 
all  reinforcement.  However,  imposed  reinforcement  bar  movements  do  not 
exceed  the  magnitude  of  their  respective  diameters,  thus  limiting  error 
to  a  minimal  degree.  Of  a  more  difficult  nature  to  evaluate  objectively, 
the  fineness  of  a  finite  element  mesh  can  influence  convergence  and 
modelling  accuracy. 

The  mesh  size  was  selected  in  this  instance  as  a  compromise 
between  realistic  accuracy  expectation  and  computer  execution  costs. 
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6. 5.1. 2  Review  of  Individual  Beam  Results:  Two  areas  of 
discrepancy  in  the  comparison  of  model  and  experimental  behaviour  are 
consistent  for  all  beams.  Model  bending  moments  at  failure  under¬ 
estimate  experimental  values,  thus  affecting  the  corresponding  torques 
in  a  similar  manner.  At  low  load  levels  when  beam  deformations  are  small, 
the  experimental  differential  rotation  measurements  are  not  reliable.  In 
the  following  examination  of  Figures  6.2  to  6.21,  the  legend  of  graph 

coordinates  for  the  reinforcement  stresses-load  plots  is  given  in  Table 
5.3. 


Beam_Rl.:  The  relevant  plotted  relationships  are  illustrated 
in  Figures  6.2,  6.3  and  6. A. 

Since  the  ratio  of  torque  to  bending  moment  is  moderate  (.5) 
throughout  the  loading  sequence,  the  model  cracking  load  is  significantly 
higher  than  the  corresponding  experimental  value  as  anticipated  in  the 
discussion  of  Section  6. 5. 1.1.  The  "translated"  computer  model  curves 
illustrated  in  Figures  6.2  and  6.3  show  a  reasonably  close  post-cracking 
stiffness  correspondence,  the  exception  being  in  the  highly  inelastic 
region  of  the  torque-rotation  relationship.  During  the  experimental 
testing  of  beam  Rl,  the  central  beam  length  over  which  the  differential 
rotation  was  calculated,  was  further  removed  from  the  influence  of  the 
substantial  stirrup-like  torsion  load  arms  than  in  the  computer  model. 
Consequently,  average  model  differential  rotation  measurements  are 
significantly  smaller  than  their  experimental  counterparts  at  high  load 
levels  when  the  restraint  of  the  torsion  arms  on  crack  widening  is  most 
pronounced.  Difference  in  cracking  load  predictions  is  also  exhibited 
in  the  reinforcement  stress  plots  of  Figures  6.4(a)  and  6.4(b). 
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Beam_R2:  The  relevant  plotted  relationships  are  illustrated 
in  Figures  6.5,  6.6,  and  6.7. 

As  the  bending  moment-deflection  ratio  at  failure  is  higher 
for  beam  R2  than  Rl,  the  degree  to  which  the  model  ultimate  strength 
predictions  underestimate  beam  capacity  is  more  pronounced.  Correspondence 
of  the  elastic  and  post-cracking  stiffnesses  is  accurate  in  both  the 
bending  moment  and  torque  graphs,  with  modest  divergence  occurring  close 
to  failure.  The  only  significant  difference  in  the  reinforcement 
stress  plots  is  that  the  eastern  leg  of  the  central  hoop  did  not  yield 
in  the  model  beam,  but  the  discrepancy  is  not  significant. 

Beam_R3:  The  relevant  plotted  relationships  are  illustrated 
in  Figures  6.8,  6.9,  and  6.10. 

For  beam  R3,  the  disproportionate  underestimation  of  the 
ultimate  torsional  capacity  compared  with  that  of  the  ultimate  bending 
moment  capacity  is  a  direct  result  of  the  loading  sequence.  In  the  load 
increments  preceding  failure,  the  ratio  of  bending  moment  to  torque  was 
equal  to  1.06.  Consequently,  premature  failure  of  the  test  specimen  had 
an  exaggerated  influence  on  the  accuracy  of  the  ultimate  torsional 
estimate.  Accurate  prediction  of  the  model  cracking  load  is  expected  as 
the  beam  was  subjected  to  bending  moment  only  at  the  formation  of  initial 
cracks.  The  reinforcement  stress  plots  also  reflect  the  close  estimate 
of  the  cracking  load. 

Beam_RA:  The  relevant  plotted  relationships  are  illustrated 
in  Figures  6.11  and  6.12. 
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During  testing,  an  unexpected  increase  in  beam  stiffness 
occurred  beyond  the  applied  bending  moment  of  1100  inch  kips.  In  retro¬ 
spect,  it  seems  that  the  nature  of  application  of  the  central  load  altered, 
reducing  the  central  bending  moment  as  illustrated  in  Fig.  6.22.  In  Fig. 
6.11,  the  "adjusted"  experimental  curve  represents  the  bending  moment- 
deflection  curve  based  on  the  assumption  that  the  central  loading  plate 
ceased  to  maintain  uniform  contact  with  the  concrete  beam  at  centre 
span,  effectively  maintaining  only  edge  contact  as  shown  in  Fig.  6.22(b). 
The  smooth  stiffness  change,  absence  of  a  marked  post-cracking  stiffness 
increase,  and  a  more  realistic  ultimate  bending  moment  capacity  suggest 
that  the  adjusted  curve  is  a  more  reasonable  representation  of  the  actual 
experimental  behaviour.  Correspondence  in  post-cracking  stiffness 
between  the  c  .uputer  model  and  adjusted  experimental  curve  is  close. 

The  disagreement  between  the  corresponding  east  hoop  stress  curves  in 
Fig.  6.12(b)  is  the  result  of  the  hoop  strain  gage  being  located  above 
its  assumed  mid-depth  position. 

Beam_R_5:  The  relevant  plotted  relationships  are  illustrated 
in  Figures  6.13,  6.14,  and  6.15. 

illustrated  in  Table  4.2,  the  initial  concrete  modulus 
for  beam  R5  was  considerably  below  the  average  modulus  of  the  other  six 
beams,  despite  the  comparable  concrete  strengths  obtained  from  cylinder 
tests.  The  low  concrete  modulus  was  undoubtedly  due  to  inadequate 
vibration  in  the  narrow  web  walls  and  bottom  flange,  whereas  the  top 
flange,  being  completely  exposed,  was  compacted  sufficiently.  Thus, 
the  modulus  value  in  Table  4.2  slightly  over-estimates  the  web  and  bottom 
flange  stiffness,  and  substantially  underestimates  the  top  flange 
stiffness.  Consequently,  the  post-cracking  stiffness  of  the  analytical 
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model  is  less  than  the  corresponding  experimental  stiffness  in  the  bending 
moment-deflection  relationship  where  the  top  flange  stiffness  is  critical, 
and  greater  in  the  torque-rotation  relationship  where  web  and  bottom 
flange  stiffness  is  influential.  The  top  longitudinal  steel  strain 
gage  was  malfunctioning  during  the  experimental  test. 

BeamVTl:  The  relevant  plotted  relationships  are  illustrated 
in  Figures  6.16,  6.17,  and  6.18. 

In  addition  to  the  expected  prediction  of  a  higher  cracking 
load,  the  only  region  of  significant  deviation  is  the  torque-rotation 
curve  segment  close  to  failure.  During  the  latter  stages  of  the  experi¬ 
mental  test,  localized  crushing  commenced  on  a  top  flange  corner  directly 
beneath  a  torsion  loading  arm.  The  crushing  was  confined  to  the  corner 
only,  but  the  trapezoidal  beam's  torsional  stiffness  was  undoubedly 
diminished,  thus  contributing  to  the  increasingly  inelastic  response  of 
the  experimental  beam. 

Beam_T^:  The  relevant  plotted  relationships  are  illustrated 
in  Figures  6.19,  6.20,  and  6.21. 

Beyond  550  inch  kips  in  the  bending  moment— deflection  curve, 
considerable  discrepancy  is  apparent  that  defies  explanation.  Deviation 
of  model  and  experimental  behaviour  displayed  in  the  torque-rotation 
graph  is  exaggerated  since  both  curves  are  highly  inelastic  immediately 
prior  to  failure  and  the  ultimate  torque  at  failure  is  small  in  comparison 
to  the  corresponding  bending  moment.  Experimental  coordinates  were  not 
plotted  in  the  latter  figure  below  a  torque  of  109  inch  kips  as  the 
results  were  inconsistent. 
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6*5. 1.3  Summary:  The  computer  model  simulation  of  the  seven 
prestressed  box  girders  tested  In  the  experimental  program  satisfactorily 
described  member  deformations  and  evaluated  ultimate  strengths  to  within 
an  acceptable  degree  of  accuracy.  To  assess  the  analytical  model's 
performance  in  its  proper  perspective,  two  important  sources  of  model 
behavioural  inadequacy  must  be  recognized.  Firstly,  the  computer  model 
consistently  overestimates  the  cracking  strength  of  "thick"  walled 
concrete  box  girders  that  are  subjected  to  any  loading  combination  that 
includes  torque.  However,  this  shortcoming  can  be  overcome  as  the  actual 
cracking  load  can  be  calculated  accurately,  and  a  "translated  computer 
model  curve"  can  be  drawn  to  illustrate  model  behaviour  if  model  and 
experimental  cracking  loads  coincide.  Secondly,  upon  development  of  a 
finite  element  mesh,  the  ultimate  bending  moment  lever  arm  is  fixed  by 
geometry,  resulting  in  an  inaccurate  estimate  of  the  ultimate  bending 
moment  capacity.  Beam  strength  under  combined  loading  that  includes 
bending  moment  is  influenced  in  a  similar  manner.  No  corrective  measures 
can  be  made  in  the  computer  model  to  negate  the  influence  of  the  latter 
anomaly . 

6.5.2  Computer  Model  Assessment  in  Light  of  Current  Theory 

6. 5. 2.1  Ultimate  Strength 

Ul^ima_te_T£rs^ional_S_treng^th:  Calculation  of  the  theoretical 
ultimate  torsional  capacity  of  both  rectangular  and  trapezoidal  beams 
has  been  made  on  the  basis  of  the  "compressive  stress  field  theory" 
or  space  truss  theory.  The  theory  is  based  on  a  kinematic  approach  in 
the  theory  of  plasticity  where  an  external  loading  is  sought  at  which 
a  mechanism  of  deformations  is  formed.  However,  the  formation  of  the 
mechanism  implicitly  assumes  that  both  stirrups  and  longitudinal  steel 
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yield  before  failure.  Although  the  space  truss  model  has  several  other 
inherent  limitations  in  its  application,  the  requirement  that  reinforcement 
proportions  be  such  that  stirrups  and  stringers  yield  prior  to  failure 
is  crucial  in  establishing  the  theory’s  validity  in  this  instance. 

Elfgren51  reports  that  several  researchers  have  confirmed 
that  the  yield  prerequisite  is  met  if  the  ratio  of  longitudinal  to 
transverse  reinforcement  is  such  that 

0.5  <  cot  aT  <  2.0  (6.14) 

where 

2A  ay 

COt  “T  =  /  FTV  •  (6.15) 

A  a  y 
w  w 

in  which  A^  =  area  of  longitudinal  reinforcement,  =  yield  stress  of 
longitudinal  reinforcement,  b’  =  width  of  beam  measured  between  corner 
stringers,  h’  =  beam  depth  measured  from  top  to  bottom  stringer,  t  = 
stirrup  spacing,  A^  =  stirrup  area,  and  oj  =  stirrup  yield  stress. 

For  the  five  rectangular  beams,  cot  aT  equals  1.96,  and  2.05 
for  the  two  trapezoidal  beams.  The  value  of  cot  aT  for  the  trapezoidal 
beams  is  of  questionable  significance  as  Eq.  6.15  was  developed  for 
rectangular  beams  only. 

On  the  basis  of  the  theoretical  verification,  a  valid  comparison 
can  be  made  of  the  computer  model  and  theoretical  ultimate  torsional 
capacity  estimates  for  the  rectangular  beams,  and  the  respective  values 
in  Table  6.2  illustrate  good  agreement.  However,  Equations  6.14  and 
6.15  indicate  that  a  comparison  of  trapezoidal  results  is  highly  question¬ 
able  and  furthermore,  experimental  and  computer  model  results  verify  that 
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the  stirrup  reinforcement  in  both  the  trapezoidal  beams  was  not  yielding 
at  failure. 

The  slightly  higher  computer  model  estimate  for  the  rectangular 
beams  is  as  anticipated,  since  only  the  yield  force  of  the  reinforce¬ 
ment  is  used  in  the  theoretical  calculation.  Unlike  the  computer  model, 
space  truss  theory  does  not  take  into  account  reinforcement  strain 
hardening  or  dowel  action.  The  additional  reserve  of  reinforcement 
strength  beyond  yielding  does  not  increase  the  ultimate  torsional  capacity 
substantially,  however,  as  the  transfer  of  load  from  concrete  to  steel 
increases  reinforcement  stresses  dramatically  as  failure  is  approached. 

Ul_t ima_t e_B^nd_ing_Momen t_  Strength :  Both  the  space  truss  and 
computer  models  cannot  accurately  estimate  ultimate  bending  moment 
capacities  as  their  respective  ultimate  bending  moment  lever  arms  are 
fixed  by  geometry.  Although  not  afflicted  by  the  latter  shortcoming, 
the  skew  bending  model  cannot  accommodate  longitudinal  web  bars  and 
tension  reinforcement  of  different  yield  strengths.  Thus,  the  equivalent 
stress  block  theory  incorporated  in  prestressed  concrete  design  is  the 
most  accurate  method  of  evaluating  the  ultimate  bending  moment  capacity. 
The  shortcoming  of  the  computer  model  in  this  respect  is  treated  in 
detail  in  Section  6. 5. 1.1. 

Ul_timaJte_Sh_ear_St^ren^th_!!  The  approach  adopted  by  Thiir liman  in 
formulating  the  theoretical  shear  capacity  of  a  reinforced  or  prestressed 
concrete  beam  is  only  applicable  to  beams  underreinforced  for  shear. 
Definition  of  an  underreinforced  beam  in  the  context  of  pure  shear  is 
identical  to  that  for  pure  torsion,  in  that  both  the  longitudinal 
stringers  and  the  transverse  reinforcement  must  yield  at  failure.  Should 
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either  reinforcement  type  not  yield  at  failure,  the  yield  criteria  are 
violated  and  the  theory  becomes  invalid. 

Unfortunately,  a  loading  pattern  could  not  be  devised  to 
enable  the  computer  model  to  predict  the  maximum  shear  resistance  of 
rectangular  beams  R3  and  R4.  The  choice  of  loading  pattern  was 
restricted  since  the  beam  span  could  not  be  reduced  without  modifying 
behaviour,  and  concentrated  loads  could  not  be  located  close  to  the 
critical  cross-section  where  shear  failure  was  anticipated  as  their 
presence  altered  the  local  shear  stress  distribution50’51.  Recognizing 
these  two  restrictions,  the  only  feasible  load  combination  produced  a 
design  region  of  maximum  shear  accompanied  by  a  large  bending  moment. 

6,5*2‘3  Combined  Loading  Interaction:  Under  combined  loading 
conditions,  the  correspondence  between  computer  model  and  theoretical 
results  is  difficult  to  analyze  objectively,  expecially  since  the 
number  of  beam  specimens  is  small  and  no  noticeable  trends  of  deviation 
are  apparent.  Strength  comparisons  for  each  individual  load  type  have 
been  made,  and  qualifications  that  arose  through  those  comparisons  are 
applicable  to  the  combined  loading  evaluation. 

Correspondence  in  the  torque-bending  moment  interaction, 
illustrated  in  Fig.  6.23,  is  good,  although  it  should  be  recognized 
that  non-dimensional  presentation  of  results  can  disguise  fundamental 
differences  in  performance.  For  example,  beams  T1  and  T2  are  both 
overreinforced  for  torsion,  but  their  interaction  coordinates  lie  close 
to  the  corresponding  theoretical  underreinforced  interaction  curve. 

Displayed  in  Fig.  6.25,  the  interaction  equations  for  beams 
R3  anc^  R4  d°  n°t  significantly  reflect  the  influence  of  the  modest  amount 
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of  shear  present  in  both  beams.  Since  the  degree  of  error  that  is 
inevitably  encountered  in  experimentation  and  computer  modelling  is  of 
the  same  order  of  magnitude  as  the  adjustment  of  the  two  interaction 
equations  for  the  presence  of  shear,  little  can  be  deduced  from  the 
proximity  of  the  computer  model  coordinates  to  the  theoretical  torque¬ 
bending  moment-shear  interaction  curve  illustrated  in  the  latter  figure. 
In  the  testing  of  beam  R3,  the  ratio  of  torque  to  bending  moment  was 
close  to  unity  in  the  load  increments  preceding  failure.  Thus,  experi¬ 
mental  and  computer  modelling  error  would  have  resulted  in  disproportion¬ 
ately  large  changes  in  the  torque  loading  at  failure. 

The  most  serious  reservation  concerning  interactive  behaviour 
is  the  obscure  definition  of  an  underreinforced  beam,  especially  when 
shear  is  present.  Current  provisions  proposed  to  prevent  premature 
failure  by  crushing  of  concrete  are  conservative  in  the  absence  of  a 
more  thorough  understanding  of  load  interaction.  In  contrast,  the 
computer  model  is  a  generalized  analytical  tool  whose  applicability  is 
not  restricted  by  loading  or  reinforcement  limitations. 

6. 5. 2. 3  Limitations  in  Application  of  Theory:  The  following 
limitations  apply  to  the  calculation  of  the  ultimate  torsional  capacity 
of  a  reinforced  concrete  beam  using  space  truss  theory: 

1.  The  beam  must  be  underreinforced:  ie.  longitudinal  and 
transverse  reinforcement  must  yield  prior  to  failure. 

2.  St.  Venant  torsion  must  be  dominant. 

3.  Along  the  beam  length,  the  cross-section  must  be  uniform. 

4.  Dowel  action  of  reinforcement  is  neglected. 
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The  theory  adopted  to  calculate  the  ultimate  bending  moment 
strength,  as  described  in  Section  6.4.1  and  6. 5. 2.1,  is  not  subjected 
to  limitation  in  its  applicability. 

Thurliman  s  approach  in  estimating  the  ultimate  shear  capacity 
of  a  reinforced  or  prestressed  beam  is  only  subject  to  two  limitations: 

1.  The  beam  must  be  underreinforced  for  shear:  ie.  longitudinal 
and  transverse  reinforcement  yield  simultaneously  at  failure. 

2.  The  top  uncracked  flange  does  not  carry  shear. 

The  interaction  equations  presented  by  Elfgren51  are  subject 
to  the  following  qualifications  in  their  application: 

1.  The  beam  cross-section  must  be  consistent  along  its  length. 

2.  Interaction  equations  for  polygonal  cross-sections  are 
not  verified. 

3.  Criteria  for  preventing  premature  failure  through  concrete 
crushing  are  not  explicit. 

The  only  limitation  or  qualification  equally  applicable  to  the 
computer  model  is  that  the  uncracked  top  flange  does  not  carry  shear. 

One  serious  qualification  not  attributed  to  theory  is  the  inaccurate 
length  of  the  bending  moment  lever  arm  at  failure.  In  conclusion,  the 
flexibility  of  the  computer  model  is  characterized  by  several  important 
areas  of  application  that  are  beyond  theoretical  capabilities: 

1.  The  member  cross-section  can  change  along  its  length 
through  the  use  of  a  general  quadrilateral  concrete  finite  element. 

2.  Any  reinforcement  arrangement  can  be  accommodated  for 
both  underreinforced  and  overreinforced  conditions. 

3.  Member  deformations  are  described  comprehensively. 
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4.  All  stress-deformation  information  is  available  at  any 
loading  stage  between  cracking  and  failure. 

5.  Indeterminate  structural  analysis  under  any  loading 
combination  is  possible. 

6. 5. 2. 4  Summary ;  In  pure  torsion,  agreement  between  the 
computer  model  and  theoretical  results  was  satisfactory  for  the  under¬ 
reinforced  rectangular  beams.  However,  since  the  two  trapezoidal  beams 
were  overreinforced  for  torsion,  a  meaningful  comparison  could  not  be 
made.  Analytical  model  estimates  of  the  ultimate  bending  moment  capacities 
of  all  beams  were  inaccurate  since  the  bending  moment  lever  arms,  fixed 
by  finite  element  mesh  geometry,  were  incorrect. 

Although  difficult  to  assess  objectively,  torque-bending 
moment  interaction  corresponded  closely.  Such  close  correspondence 
was  not  apparent  in  the  torque-bending  moment-shear  interactive  behaviour 
of  beams  R3  and  R4. 

Objective  assessment  of  the  computer  model  results  through 
comparison  with  the  corresponding  theoretical  predictions  was  not 
conclusive  as  the  theory  is  limited  in  its  range  of  application,  and 
complex  interactive  behaviour  is  presented  theoretically  in  a  general 
form  that  does  not  permit  explicit  comparison. 
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WEB  REINFORCEMENT 
CONCENTRATED  AT 


(a)  Unstable  cracked 

reinforced  concrete  element 


WEB  REINFORCEMENT 
DISTRIBUTED  ACROSS  CRACK 


(b)  Stable  cracked 

reinforced  concrete  element 


SHEAR 

FORCES 


i 


r 


*  Note:  Web  reinforcement  is  uniformly  distributed  across  element, 
but  when  crack  forms  through  centroid,  the  reinforcement 
is  effectively  distributed  across  the  crack  only. 


FIG.  6.1(A)  METHOD  OF  MODELLING  WEB  REINFORCEMENT 


FIG.  6.1(B)  TRANSLATION  PROCEDURE 
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FIG.  6.2  MODEL  AND  TEST  BENDING  MOMENT-DEFLECTION  CURVES 
FOR  BEAM  R1 
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FIG.  6.4(A)  TEST  AND  MODEL  LONGITUDINAL  CONVENTIONAL 
REINFORCEMENT  STRESSES  FOR  BEAM  R1 


FIG.  6.4(B)  TEST  AND  MODEL  HOOP  REINFORCEMENT  STRESSES 
FOR  BEAM  R1 
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FIG.  6.5  MODEL  AND  TEST  BENDING  MOMENT-DEFLECTION  CURVES 
FOR  BEAM  R2 
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Experimental  ultimate  torque  =  372"k 

Computer  model  failure  torque  =  336”k 

./•r 
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(•>  1 Q  Experimental  curve 


Computer  model  curve 

Translated  computer  model  curve 
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:L  AND  TEST  TORQUE-ROTATION  CURVES  FOR  BEAM  R2 
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FIG.  6.7(A)  TEST  AND  MODEL  LONGITUDINAL  CONVENTIONAL 
REINFORCEMENT  STRESSES  FOR  BEAM  R2 


FIG.  6.7(B)  TEST  AND  MODEL  HOOP  REINFORCEMENT  STRESSES 
FOR  BEAM  R2 


1400  Experimental  ultimate  moment  =  1357"k 
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FIG.  6.8  MODEL  AND  TEST  BENDING  MOMENT-DEFLECTION  CURVES 
FOR  BEAM  R3 


Experimental  ultimate  torque  *  532"k 
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FIG.  6.10(A)  TEST  AND  MODEL  LONGITUDINAL  CONVENTIONAL 
REINFORCEMENT  STRESSES  FOR  BEAM  R3 


FALSE 

ORIGIN 


FIG.  6.10(B)  TEST  AND  MODEL  HOOP  REINFORCEMENT  STRESSES 
FOR  BEAM  R3 


Experimental  ultimate  moment  =  1624"k 
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FIG.  6.11  MODEL  AND  TEST  BENDING  MOMENT  DEFLECTION  CURVES  FOR  BEAM  R4 
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FIG.  6.12(A)  TEST  AND  MODEL  LONGITUDINAL  CONVENTIONAL 
REINFORCEMENT  STRESSES  FOR  BEAM  R4 
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FIG.  6.12(B)  TEST  AND  MODEL  HOOP  REINFORCEMENT  STRESSES 
FOR  BEAM  R4 


FAILURE  LOAD  -  50  k 


BENDING  MOMENT  (inch  kips) 
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FIG.  6.13  MODEL  AND  TEST  BENDING  MOMENT-DEFLECTION  CURVES 
FOR  BEAM  R5 


TORQUE  ( inch  kips ) 
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FIG.  6.14  MODEL  AND  TEST  TORQUE-ROTATION  CURVES  FOR  BEAM  R5 
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FIG.  6.15(A)  TEST  AND  MODEL  LONGITUDINAL  CONVENTIONAL 
REINFORCEMENT  STRESSES  FOR  BEAM  R5 


FIG.  6.15(B)  TEST  AND  MODEL  HOOP  REINFORCEMENT  STRESSES 
FOR  BEAM  R5 
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FIG.  6.16  MODEL  AND  TEST  BENDING  MOMENT  DEFLECTION  CURVES 
FOR  BEAM  T1 
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FIG.  6.17  MODEL  AND  TEST  TORQUE-ROTATION  CURVES  FOR  BEAM  T1 
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FIG.  6.18(A)  TEST  AND  MODEL  LONGITUDINAL  CONVENTIONAL 
REINFORCEMENT  STRESSES  FOR  BEAM  T1 


FIG.  6.18(B)  TEST  AND  MODEL  HOOP  REINFORCEMENT  STRESSES 
FOR  BEAM  T1 


BENDING  MOMENT  (inch  kips) 
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FIG.  6.19  MODEL  AND  TEST  BENDING  MOMENT-DEFLECTION  CURVES 
FOR  BEAM  T2 


200  "  Experimental  ultimate  torque  =  196.5"k 
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FIG.  6.21  TEST  AND  MODEL  LONGITUDINAL  CONVENTIONAL 
REINFORCEMENT  STRESSES  FOR  BEAM  T2 
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3" 

PLASTER  OF  PARIS 

' 

' 

LOADING  PLATE 

PAD'S\k 
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r* 

v  TEST  1 

BEAM  , 

f 
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FIG.  6.22  CENTRAL  LOAD  VARIATION  FOR  BEAM  R4 
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FIG.  6.23  TORQUE-BENDING  MOMENT  INTERACTION  DIAGRAM  FOR 
MODEL  EVALUATION 
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FIG.  6.24  INTERACTION  DIAGRAM  FOR  BENDING  MOMENT  AND  SHEAR 


164. 


INTERACTION  EQUATIONS  FOR  BEAM  R3: 

MODE  t :  -fi-  +  0.2377  (^-)2  =  0.994 

MODE  c:  -4.207  -jjfc-  +  (-|~)2  =  0.976 
MODE  *  :  ~  =  | ,6 

'O 


INTERACTION  EQUATIONS  FOR  BEAM  R4  : 

MODE  t:  -!£-  +  0.2377  =  0.985 

MODE  c:  -4.207  =0.936 

MODE  s  :  —I—  =  |.6 

To 

FIG.  6.25  ADJUSTED  TORQUE-BENDING  MOMENT-SHEAR  INTERACTION 
DIAGRAM  FOR  BEAMS  R3  AND  R4 


CHAPTER  VII 


CONCLUSION  AND  SUMMARY 

7 . 1  Principal  Implications  of  Comparison 

Of  the  two  sources  of  comparison  that  were  utilized  in  the 
evaluation  of  the  analytical  model,  the  results  of  the  small,  comprehen¬ 
sive  experimental  program  provided  the  better  basis  for  a  conclusive 
assessment. 

From  experimental  comparisons,  two  potential  sources  of  model 
inaccuracy  were  isolated;  underconservative  estimation  of  the  cracking 
load  under  torsion  loading  conditions,  and  inaccurate  representation  of 
the  bending  moment  lever  arm  at  ultimate  failure.  The  former  qualifi¬ 
cation  of  behaviour  can  be  overcome  if  the  variation  of  torsional  shear 
stress  across  the  wall  thickness  is  taken  into  account.  To  reflect  the 
more  accurate  prediction  of  cracking  load,  the  computer  model  results 
are  transformed  accordingly  through  the  translation  of  the  post-cracking 
sections  of  the  bending  moment-deflection  and  torque-rotation  curves. 
Inaccurate  prediction  of  the  ultimate  bending  moment  capacity  cannot  be 
resolved  as  the  maximum  lever  arm  length  is  fixed  by  geometry,  but  this 
deficiency  is  only  of  significance  in  underreinforced  beams  where  the 
ratio  of  compression  flange  thickness  to  beam  depth  is  unusually  high. 
Since  the  modest  degree  of  discrepancy  between  experimental  and  model 
results  is  principally  derived  from  the  latter  qualification,  accurate 
simulation  of  beam  behaviour  can  be  attained  for  concrete  box  girders 
of  lower  wall  thickness  to  depth  ratios. 
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The  principal  conclusion  that  arose  from  the  model  assessment 
in  light  of  current  theory  was  the  general  applicability  of  the  analytical 
model  achieved  through  freedom  from  restrictive  assumptions.  Where 
valid  comparisons  can  be  made,  agreement  between  analytical  and  theoretical 
results  is  good.  Theoretical  interactive  behaviour,  presented  in  its 
dimensionless  equation  form,  does  not  permit  an  explicit  evaluation  of 
model  results,  and  can  indeed  be  misleading  because  of  the  dimensionless 
form  of  presentation. 

Consideration  of  the  significant  aspects  of  the  assessment 
procedure  affirm  the  value  of  the  computer  model  as  a  capable,  versatile, 
^n^-lytical  tool  in  its  context  of  the  analysis  of  concrete  box  girders 
acted  upon  by  torsion,  bending,  and  shear.  Qualification  of  its  accuracy 
is  isolated  to  one  potential  source  of  error,  and  the  range  of  applica¬ 
tion  is  almost  without  restriction.  Characteristic  of  all  finite  element 
modelling,  the  reliability  of  output  information  is  reflected  in  the 
quality  of  input  specifications. 

7.2  Application  of  the  Analytical  Model 

Within  their  respective  confines  of  validity  determined  by 
explicit  assumptions,  theoretical  strength  predictions  discussed  in 
Chapter  6  are  not  in  conflict  with  analytical  model  results.  The  import¬ 
ant  qualification  concerning  the  analytical  value  of  current  theory  is 
the  restrictive  nature  of  assumptions  that  limit  the  regions  of  theoret¬ 
ical  application.  As  specified  in  Section  6.5.2. 3,  the  following  computer 
model  capabilities  are  beyond  the  range  of  theory  application: 
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1.  Beams  may  contain  any  level  of  reinforcement,  to  the  extent  of 
being  overreinforced. 

2.  Cross-sectional  geometry  may  vary  along  the  member's  length. 

3.  St.  Venant  torsion  need  not  be  dominant. 

4.  Member  deformations  are  described  comprehensively  at  all  load 
levels  from  commencement  of  loading  to  failure. 

5.  The  stress  levels  of  all  component  materials  are  evaluated 
throughout  the  loading  sequence. 

6.  Indeterminate  analysis  under  any  loading  condition  is  possible. 

In  its  present  form,  the  computer  model  cannot  permit  variation  in 
cross-sectional  geometry  along  its  length  as  a  rectangular  concrete 
finite  element  has  been  used.  Replacement  of  the  rectangular  element 
with  a  plane  stress  general  quadrilateral  element38  that  has  the  same 
degrees  of  freedom  will  overcome  this  shortcoming. 

The  flexibility  and  fully  comprehensive  nature  of  the  analyti¬ 
cal  model  can  be  a  valuable  asset  in  the  design  process.  Although  the 
model  is  not  suited  to  direct  incorporation  in  preliminary  design, 
final  design  proposals  can  be  checked  thoroughly  through  the  use  of  the 
analytical  model  after  material  and  geometrical  parameters  have  been 
chosen.  Since  current  practice  is  increasingly  oriented  toward  ultimate 
strength  design,  a  certain  degree  of  concrete  cracking  is  often  tolerated 
at  extreme  service  load  conditions.  Evaluation  of  reinforcement  and 
concrete  stress  levels  under  such  conditions,  though  beyond  theoretical 
capabilities,  is  readily  achieved  by  the  computer  model  which,  in 
addition,  can  provide  accurate  estimation  of  structural  deformations. 

In  statically  indeterminate  structures,  the  analytical  model  may  well 
represent  the  only  reliable  method  of  analysis. 
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Enthusiastic  adoption  of  the  analytical  model  approach  must  be 
tempered  through  recognition  of  the  inherent  dependence  of  model  accuracy 
on  quality  of  input  specifications.  Evaluation  of  failure  load  con¬ 
ditions  is  most  strongly  influenced  by  beam  geometry  and  reinforcement 
strength,  the  effect  of  concrete  compressive  strength  being  more  pro¬ 
nounced  as  the  overreinforced  condition  is  approached.  In  the  accurate 
determination  of  deformations,  the  initial  concrete  modulus  is  the  most 
crucial  parameter,  assuming  that  the  cracking  load  is  estimated  closely 
through  accurate  monitoring  of  concrete  tensile  strength  and  shrinkage 
stresses.  Invariably,  concrete  strength  parameters  can  be  notoriously 
inconsistent  if  care  is  not  exercised  in  their  evaluation.  Therefore, 
suitability  of  the  computer  model  method  must  be  viewed  within  the 
context  of  input  parameter  accuracy,  as  the  highly  sophisticated  and 
comprehensive  nature  of  finite  element  analysis  can  be  completely 
negated  by  erroneous  input. 

7 . 3  Summary 

The  principal  objective  of  this  thesis  topic  is  the  development 
of  an  analytical  computer  model  that  can  analyze  reinforced  or  pre¬ 
stressed  concrete  box  girders  of  arbitrary  cross-section  for  any  loading 
combination  of  bending  moment,  torque,  and  shear.  Of  the  many  refined 
capabilities  of  the  finite  element  model,  prediction  of  deformations  in 
the  post-cracking  region  has  been  chosen  as  the  main  premise  on  which 
performance  of  the  model  is  evaluated. 

Linear  concrete  segments  of  the  box  girder  are  represented  in 
the  model  by  higher-order  plane  stress  rectangular  finite  elements  that 
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are  assigned  twelve  degrees  of  freedom,  two  translational  and  one 
rotational  degree  of  freedom  at  each  element  node.  Reinforcement  is 
represented  by  one-dimensional  bar  elements.  To  enable  the  complete 
load-deformation  path  to  be  described  from  onset  of  loading  through  to 
failure,  the  loading  sequence  is  an  incremental,  iterative  one  where 
the  probability  of  material  behavioural  deviation  in  successive  load 
increments  is  reduced  by  the  Runge-Kutta  method.  Should  significant 
material  deviation  be  detected,  the  modified  Newton-Rapson  method 
restores  equilibrium  in  an  iterative  procedure. 

Concrete  is  modelled  as  a  non-linear,  anisotropic  material, 
and  following  cracking,  its  effective  shear  modulus  in  the  reformulated 
constitutive  relationship  is  redefined  on  aggregate  interlock  and  dowel 
action  considerations.  The  important  aspects  of  diaphragm  action  are 
treated  comprehensively  in  the  analytical  model  in  the  simulation  of 
cross-sectional  distortional  stiffness,  longitudinal  warping  restraint, 
and  intrinsic  transverse  shear  rigidity  of  box  girders  without  diaphragms. 
In  the  absence  of  theoretical  or  experimental  information,  an  auxiliary 
finite  element  model  has  been  employed  to  investigate  the  longitudinal 
warping  resistance  of  "thick"  diaphragms. 

To  test  the  capabilities  of  the  computer  model  rigorously, 
a  literature  survey  was  conducted  to  find  experimental  results  of  multi- 
celled  box  girders  subjected  to  torsion,  bending,  and  shear  where 
testing  was  pursued  to  failure.  No  suitable  reference  was  found,  and 
consequently,  a  small  experimental  program  was  undertaken  in  which  seven 
double-celled,  prestressed  concrete  box  girders,  five  rectangular  and 
two  trapezoidal  in  cross-section,  were  tested  to  failure  under  varying 


. 
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combinations  of  the  latter  three  load  types. 

To  evaluate  the  accuracy  of  the  analytical  model  predictions, 
corresponding  computer  model  and  experimental  results  were  compared  in 
detail  through  examination  of  the  respective  bending  moment-deflection, 
torque-rotation,  and  reinforcement  stresses-load  relationships.  Where 
applicable,  the  model  results  were  also  assessed  in  light  of  current 
theoretical  predictions.  From  an  overall  perspective,  performance  of 
the  analytical  model  was  very  satisfactory. 

7.4  Conclusion 


The  object  of  this  thesis  has  been  achieved  in  the  development 
a  finite  element  analytical  model  that  can  analyze  concrete  box 
girders  of  arbitrary  cross-section  subjected  to  torsion,  bending,  and 
shear.  Being  afflicted  by  only  one  shortcoming  that  is  not  significant 
in  members  of  common  cross-sectional  geometry,  the  computer  model  is  a 
flexible,  comprehensive  mode  of  analysis  whose  range  of  application  is 
limited  by  few  constraints.  Of  the  numerous  material  parameters  that 
govern  the  model  s  behaviour,  only  the  phenomenon  of  dowel  action  is 
ill-defined  as  certain  aspects  of  its  contribution  to  shear  strength 
have  not  as  yet  been  clarified.  In  the  research  of  concrete  box  girder 
behaviour,  exhausive  time-consuming  experimental  test  programs  are  now 
superseded  in  almost  every  respect  by  the  analytical  model.  Considering 
the  restrictive  assumptions  in  current  theory,  the  computer  model  greatly 
increases  the  analytical  scope  and  design  capability  in  this  field  of 
study. 
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APPENDIX  A 


SIGN  CONVENTIONS  AND  SYSTEMS 

(1)  Global  Axes  Directions 

The  global  x  axis  must  be  in  the  beam's  longitudinal  direction, 
and  the  global  y  axis  in  the  upward  vertical  direction.  Having  chosen 

the  x  and  y  directions,  the  global  z  axis  is  defined  by  the  right  hand 
axis  convention. 

(2)  Node  Numbering 

(a)  Rectangular  concrete  elements 
Refer  to  Figure  A-l. 

(b)  Reinforcement  elements 

The  node  numbering  sequence  should  be  in  the  positive 

global  axis  direction.  A  vertical  inclined  bar  should  be  numbered 
in  the  upward  direction. 

(c)  Bond  spring  linkages 

No  predefined  node  numbering  convention 


(3)  Nodal  Forces  and  Displacements 


Positive  in  positive  global  axis  direction. 


(4)  Nodal  Moments  and  Rotations 


Positive  when  acting  clockwise  looking  along  global  axis 


in  positive  direction. 
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(5)  Stresses  and  Strains 

(a)  Direct  stresses  and  strains 

Tension  positive,  compression  negative 

(b)  Shear  stresses  and  strains 

Positive  when  acting  on  a  positive  face  in  a  positive  global 
axis  direction,  or  a  negative  face  in  a  negative  global  axis 
direction. 

(6)  Mohr’s  Circle 

Same  sign  convention  for  direct  stresses  and  strains  as  above 
is  used.  However,  shear  stresses  and  strains  are  positive  when  the 
couple  acts  in  the  clockwise  direction. 

(7)  Dimensional  Parameter  Units 

Length  -  inches 
Force  -  pounds  wt. 

Rotation  -  degrees  (CANGLE(NEL) ) 


180. 


Y 


Z 


FIG.  A-1  NODE  NUMBERING  SEQUENCE  FOR  RECTANGULAR 
CONCRETE  ELEMENTS 
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APPENDIX  E 


SUBROUTINE  LOGIC  NOTES 

The  purpose  of  this  appendix  is  to  clarify  any  subroutine 
logic  development  that  might  not  be  sufficiently  illustrated  by  the 
numerous  comment  cards  inserted  throughout  the  program.  The  comment 
card  heading  in  each  subroutine  states  the  routine's  principal 

function.  Those  routines  warranting  comment  are  listed  in  alphabetical 
order. 

Subroutine  DOWEL 


To  detect  the  failure  of  a  dowel  mechanism  across  a  concrete 
crack,  the  dowel  displacement  of  the  adjacent  crack  surfaces  at  the 
concrete  element's  centroid  is  calculated  and  compared  to  the  critical 
dowel  displacement  value  DF. 

Subroutine  FORCE 

When  a  bi-linear  stress-strain  curve  is  assumed  for  the 
conventional  reinforcement,  considerable  difficulty  is  experienced 
in  modelling  the  reinforcement's  behaviour  beyond  yielding.  Inevitably, 
an  excessive  number  of  modified  Newton-Rapson  iterations  is  required 
to  produce  convergence,  especially  when  the  second  linear  segment  of 
the  stress-strain  curve  is  almost  perfectly  plastic.  To  eliminate  this 
particular  cause  of  time  consuming  and  costly  iteration,  a  yielding 
conventional  reinforcement  bar  is  modelled  in  the  total  stress  condition 
as  a  bar  of  zero  stiffness  with  externally  applied  nodal  forces  at 
each  bar  node  representing  the  influence  of  the  bar  of  the  remainder 
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of  the  structure.  In  the  incremental  stress  condition,  the  bar’s 
stiffness  corresponds  to  the  slope  of  the  strain-hardening  portion  of 
the  stress-strain  curve.  For  the  range  of  strain-hardening  usually 
encountered,  this  approach  simulates  reinforcement  yielding  accurately. 

i 

Subroutine  HORIZ 

Before  the  stiffness  terms  of  a  horizontal  rectangular  element 
are  added  into  the  in-core  total  stiffness  matrix,  the  component  (3x3) 
element  stiffness  blocks  must  be  transformed  as  the  local  element 
axes,  as  shown  in  Figure  A-l,  do  not  correspond  to  the  global  axis 
orientation.  The  sign  of  four  stiffness  terms  is  changed  as  a  result 
of  the  transformation  process. 

Subroutine  NDCONC 

The  formulations  used  in  this  subroutine  are  given  in  detail 
in  Sections  3. 4. 2.1  and  3. 4. 2. 6.  In  deriving  the  stiffness  of  a 
cracked  concrete  element,  the  constitutive  matrix  is  formulated  for 
axes  in  the  crack  and  orthogonal  to  crack  directions.  Consequently, 
transformation  of  the  constitutive  matrix  must  be  undertaken. 

Subroutine  SHRMOD 

The  shear  rigidity  of  a  cracked  concrete  element  in  a  direction 
parallel  to  the  crack  direction  is  the  sum  of  the  aggregate  interlock 
and  dowel  stiffnesses.  The  empirical  derivations  of  the  two  contributing 
stiffnesses  are  given  in  Sections  3. 4. 2. 4  and  3. 4. 2. 5.  In  formulation 
of  the  dowel  stiffness,  the  magnitude  is  not  permitted  to  exceed 
twenty  percent  of  the  aggregate  interlock  value. 
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Subroutine  WARP 

The  warping  resistance  of  "thick"  diaphragms  is  greatly 
overestimated  by  classical  thin  plate  theory,  as  detailed  in  Section 
3. 4. 2. 8.  However,  the  warping  stiffness  formulation  used  in  this 
subroutine  is  the  classical  form  given  in  3. 4. 2. 7,  with  real  behaviour 
reflected  by  the  insertion  of  an  equivalent  diaphragm  thickness. 
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210,  62,  15376 
40,  40, 
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ELEMENT  NUMBER  BOND  SPRING  LINKAGE 
(a)  Exploded  view  of  concrete  web  and  flange  elements 


I  'll 

(b)  Diaphragm  elements 


FIG.  F-1  BEAM  MESH  FOR  INPUT  DATA  FILE  EXAMPLE 


APPENDIX  G 
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The  purpose  of  this  appendix  is  to  assist  the  computer 
program  user  in  the  compilation  of  the  input  data  file.  Units  of  inches 
and  pounds  are  used  throughout,  and  if  large  numbers  are  to  be  read  in, 
the  E  format  should  be  used.  The  user  should  be  especially  conscious 
of  the  sign  conventions  listed  in  Appendix  A.  In  referring  to  several 
input  variables  in  the  description  below,  the  symbolic  name  may  be 
used  for  the  sake  of  brevity.  The  alphabetic  listing  of  all  symbolic 
names  and  their  corresponding  codes,  where  applicable,  are  given  in 
Appendix  B.  Complete  freedom  in  input  data  entry  is  provided  through 
the  use  of  a  semi-f reef ield  format,  where  commas  are  used  to  separate 
adjacent  entries,  as  illustrated  in  the  sample  input  listing  in 
Appendix  F.  A  sufficient  number  of  columns  have  been  allocated  for  the 
entry  on  all  variables  of  realistic  size,  but  reference  to  the  format 
statements  in  subroutine  READIN  listed  in  Appendix  D  will  clarify  the 
maximum  number  of  assigned  columns  to  accommodate  the  full  digit  length. 
Following  the  description  of  each  input  card  or  assemblage  of  associated 
cards,  the  significance  and  qualifications  that  pertain  to  the  important 
variable  entries  will  be  commented  upon. 

The  data  cards  must  occur  in  the  following  order: 

1 .  Echo  Check  Card 

Order  of  Entries  Description 

1  I WRITE 

Comments:  If  IWRITE  =  1,  the  input  data  will  be  printed  out  so  that  it 

can  be  quickly  checked . 
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2.  Number  of  Finite  Elements  Card 


Order  of  Entries 

Description 

1 

NELT 

2 

NELCHK 

Comments:  (a)  If  it  is  recognized  that  all  concrete  elements  will 

not  crack,  the  element  numbering  system  should  be  chosen  such  that 
only  the  first  NELCHK  number  of  elements  will  be  checked  for  cracking. 

(b)  NELT  does  not  include  the  number  of  diaphragm  elements. 

3.  Concrete  Finite  Elements  Description  Cards 

One  card  for  each  concrete  finite  element  in  ascending  order. 


Order  of  Entries 

Description 

1 

INDCEL  (NEL) 

2 

INELSZ  (NEL) 

3 

INELTY  (NEL) 

4 

ELTHN  (NEL) 

5  First  element  node  number 


6 

Second  node  number 

7 

Third  node  number 

8 

Fourth  node  number 

9 

INDOWL  (NEL) 

10 

INMESH  (NEL) 

11 

INSZMS  (NEL) 

12 

WIDTHC  (NEL) 

Comments:  (a)  Concrete  elements  of  the  same  dimensions,  thickness, 

and  percentage  of  reinforcement  are  assigned  the  same  arbitrary  integre 


number  INELSZ  (NEL) . 
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(b)  Only  those  web  elements  that  adjoin  a  tension  flange 
element  will  develop  dowel  shear  resistance. 

(c)  Elements  that  have  the  same  percentages  of  steel  mesh 
reinforcement  are  assigned  the  same  arbitrary  integre  number  INSZMS  (NEL) 

(d)  The  effective  dowel  width  of  the  beam  is  distributed 
proportionately  between  the  web  elements  that  adjoin  the  tension  flange. 


4.  Concrete  Shrinkage  Data  Card 

Order  of  Entries 

Description 

1 

NELTOP 

2 

SIGXT1 

3 

SIGXT2 

4 

SIGXB1 

5 

SIGXB2 

6 

SIGXS1 

7 

SIGXS2 

8 

EXTl 

9 

EXT2 

10 

EXBl 

11 

EXB2 

12 

EXS1 

13 

EXS2 

5.  Alternate  Top  and  Bottom  Flange  Shrinkage  Elements  Cards 

Alternate  cards  for  top  and  bottom  flange  concrete  elements 


that  develop  shrinkage  stresses. 
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The  first  card  will  list  ten  top  flange  concrete  elements, 
the  second  card  ten  bottom  flange  concrete  elements.  This  sequence 
of  cards  is  continued  until  all  flange  elements  subjected  to  shrinkage 
stresses  are  listed.  The  last  two  cards  will  contain  at  least  one 
entry,  but  will  not  exceed  ten  entries. 

6.  Number  of  Elements  with  Steel  Mesh  Card 

The  only  entry  on  the  card  is  NELSM. 

Comments:  It  is  imperative  for  the  proper  functioning  of  the  program 
that  all  transverse  shear  or  torsion  reinforcement  be  represented  by 
an  equivalent  steel  mesh  distributed  throughout  the  element. 

7 .  Number  of  Layers  in  Steel  Mesh  Card 

The  only  entry  is  NDIRNS . 

Comments:  If  NELSM  =  0,  omit  this  card  and  proceed  to  card  9. 

8.  Steel  Mesh  Description  Cards 

Two  cards  for  each  element  containing  a  steel  mesh.  The  element  number 
is  contained  on  the  first  card.  The  second  card  is  as  below: 

Order  of  Entries  Description 

1  Steel  mesh  percentages  for  all  layers 

2  Corresponding  inclinations 
Comments:  A  steel  mesh  layer  is  defined  as  a  discrete  system  of  parallel 
bars. 

9.  Number  of  Diaphragms  Card 


The  only  entry  is  NELD. 
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Comments:  (a)  NELD  is  the  total  sum  of  actual,  equivalent,  and 

warping  diaphragm  elements. 

(b)  If  NELD  =  0,  proceed  directly  to  card  11. 


10.  Diaphragm  Description  Cards 

One  card  for  each  diaphragm. 

Order  of  Entries 
1 
2 

3 

4 

5 

6 

7 

8 

Comments:  Identical  diaphragm 
integre  value. 


Description 

INELTY  (NEL) 

ELTHN  (NEL) 

First  element  node  number 
Second  element  node  number 
Third  element  node  number 
Fourth  element  node  number 
NDREF  (NEL) 

DPMMOD  (NEL) 

ments  have  the  same  NDREF  (NEL) 


11.  Number  of  Reinforcement  Elements  Card 

The  only  entry  is  NREO. 

Comments:  NREO  is  the  sum  of  all  prestress  and  conventional  longitudinal 
reinforcement  elements  and  bond  linkages. 


12.  Reinforcement  Type  Description  Cards 


One  card  for  each  reinforcement  element. 
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Order  of  Entries  Description 

1  First  reinforcement  node  number 

2  Second  reinforcement  node  number 

3  Reinforcement  element  type  indicator 

Comments:  (a)  The  reinforcement  node  numbers  must  be  numbered  in  the 

positive  axis  direction. 

(b)  The  following  individual  reinforcement  cards  must  be 
in  the  same  order. 


13 .  Conventional  Longitudinal  Reinforcement  Card 


This  card  is  read  if  INRTY(NR)  =  -1. 
Order  of  Entries 
1 
2 


Description 

INRDN  (NR) 
RAREA  (NR) 


14.  Prestress  Reinforcement  Card 


This  card  is  read  if  INRTY 
Order  of  Entries 
1 
2 

3 

4 


(NR)  =  0. 

Description 
INRDN  (NR) 
RAREA  (NR) 
TSGPRE  (NR) 
TEPRE  (NR) 


15 .  Bond  Linkage  Card 

This  card  is  read  if  INRTY  (NR)  =  1. 

The  only  entry  is  CAREA  (NR) . 
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16.  Number  of  Finite  Element  Nodes  Card 

The  only  entry  is  NNODES. 

Comments:  Bond  spring  linkage  elements  adjoin  adjacent  concrete  and 
reinforcement  nodes.  However,  no  bond  spring  linkages  are  provided 
at  beam  ends. 


17 .  Finite  Element  Node  Description  Cards 


One  card  for  each  node. 

Order  of  Entries 
1 
2 

3 

4 


Description 
ICNODE  (I) 

X  (I) 

Y  (I) 

Z  (I) 


Comments:  (a)  Nodes  adjoining  two  non-planar  elements  must  be 

designated  as  corner  nodes:  ie.  they  possess  five  degrees  of  freedom. 

(b)  Internal  node  numbers  connecting  web  and  diaphragm 
elements  are  assigned  ICNODE  (I)  =2:  ie.  they  have  four  degrees  of 
freedom. 


(c)  Node  adjoining  co-planar  elements  are  designated  as 
interior  nodes,  with  three  degrees  of  freedom. 

(d)  The  X,  Y,  and  Z  coordinates  are  global  axis  coordinates, 
the  direction  of  the  global  axes  being  defined  in  Appendix  A. 


18 .  Strength  Moduli  Card 

Order  of  Entries 

1 


Description 

CONMOD 


2 


REOMOD 


■ 
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3 

PREMOD 

4 

SLPMOD 

5 

SMSMOD 

19.  Ultimate  Strength  Card 

Order  of  Entries 

Description 

1 

FC 

2 

FT 

3 

FAGG 

4 

FUR 

5 

FUP 

6 

FUMS 

20.  Ultimate  Strain  Cards 

Order  of  Entries 

Description 

1 

EEREO 

2 

EEPRE 

3 

ECULT 

4 

ERULT 

5 

EPULT 

6 

EMSULT 

7 

ESLIP 

8 

DF 
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21.  Poisson  Ratios  Card 


Order  of  Entries 
1 
2 

3 

4 


Description 

PI 

P2 

P3 

PU 


22.  Material  Deviation  Card 


Order  of  Entries 

1 

2 

3 

4 

5 

6 

7 

8 


Description 

CONDEV 

DEVCON 

DEVREO 

PREDEV 

SLPDEV 

IDEV 

AVCSP 

RELAX 


Comments:  (a)  Entries  1  to  5  specify  the  maximum  allowable  deviation 

percentages  in  material  behaviour  that  is  permitted  before  corrective 
measures  are  taken.  Only  CONDEV  and  PREDEV  of  the  five  variables 
mentioned  above  have  a  direct  influence  on  the  decision  to  invoke  the 
iterative  process.  Obviously,  if  the  tolerances  expressed  by  CONDEV 
and  PREDEV  are  too  demanding,  the  expensive  iterative  process  will  be 
undertaken  frequently.  Thus,  a  compromise  has  to  be  struck  between 
the  maximum  allowable  deviation  that  is  acceptable  and  the  cost 
incurred  if  that  maximum  level  of  deviation  is  to  be  enforced.  Values 
used  in  this  application  of  the  computer  model  were  CONDEV  =  8%  and 


■ 
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PREDEV  =  10%.  Severe  tolerances  expressed  by  DEVCON  and  DEVREO  do 
not  Increase  computing  cost  significantly,  and  it  is  recommended 
that  these  values  be  kept  small. 

(b)  IDEV  is  the  critical  weighted  iteration  integre  that 
determines  whether  the  iterative  procedure  will  be  invoked.  Each 
important  aspect  of  structural  deviation  is  given  a  weighting,  as 
declared  in  the  following  card.  If  the  sum  total  of  the  deviation 
counters  exceeds  IDEV,  a  cycle  of  the  iterative  process  will  be  initiated. 

(c)  RELAX  is  the  relaxation  factor  employed  to  improve  the 
rate  of  convergence  in  the  highly  inelastic  segments  of  reinforcement 
stress-strain  curves  during  the  iterative  process.  The  variable  will 
invariably  exceed  unity,  but  its  optimum  value  is  very  much  dependent 
upon  the  nature  of  the  analysis. 

23.  Deviation  Weighting  Card 

Order  of  Entries  Description 

1  MAXIT 

2  NWTCON 

3  NWTREO 

4  NWTPRE 

Comments:  If  allowable  behavioural  deviation  is  exceeded  in  a  concrete 

element,  a  conventional  reinforcement  or  steel  mesh  element,  or  a 
prestress  reinforcement  element,  the  deviation  integre  counter  is 
increased  by  the  addition  of  NWTCON,  NVTREO ,  or  NWTPRE  respectively. 

The  deviation  integre  counter  within  the  program  is  the  variable  ID. 

Thus,  when  ID  exceeds  IDEV  after  the  application  of  a  load  increment, 
iteration  is  commenced  until  deviation  is  reduced  to  the  extent  that  ID 
is  less  than  or  equal  to  IDEV. 
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24 .  Number  of  Load  Types  Card 

The  only  entry  on  the  card  is  NLDTY. 

Comments:  A  load  type  is  one  where  the  relationship  of  the  individual 
loads  within  successive  load  increments  remains  constant. 

25 .  General  Load  Information  Card 

One  card  for  each  load  type. 

Order  of  Entries  Description 

1  NINCRT  (I) 

2  NLOADS  (I) 

Comments:  (a)  For  each  load  type,  this  card  and  the  cards  of  item  26 

are  read  in. 

(b)  The  application  of  prestress  forces  comprises  the 

first  load  increment  in  the  study  of  prestressed  concrete  beam  behaviour. 

(c)  Load  increments  can  be  heavy  in  the  elastic  and  early 
post-cracking  regions,  but  should  become  progressively  lighter  as 
behaviour  becomes  more  inelastic.  Small  load  increments  are  necessary 
close  to  failure  if  the  ultimate  load  conditions  are  to  be  estimated 
within  close  bounds. 

26.  Individual  Load  Description  Cards 

One  card  for  each  nodal  load. 

Order  of  Entries 

1 
2 


Description 
NODER  (I) 
KODER  (I) 


3 


VALUER  (I) 
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27 .  Number  of  Boundary  Conditions  Card 

The  only  entry  on  this  card  is  NDISPL. 

28 .  Individual  Boundary  Conditions  Cards 

One  card  for  each  boundary  condition. 

Order  of  Entries  Description 

1  NODED  (I) 

2  KODED  (I) 

3  VALUED  (I) 

29.  Number  of  Concrete  Elements  Display  Card 

The  only  entry  on  this  card  is  PNELT  (integre) 

30.  Concrete  Elements  Display  Card 

All  the  element  numbers  of  those  elements  whose  centroidal 
stresses  will  be  displayed,  are  listed  on  this  card.  (integres) 

31.  Number  of  Steel  Mesh  Elements  Display  Card 

The  only  entry  on  this  card  is  PNELSM  (integre) . 

32 .  Steel  Mesh  Elements  Display  Card 

All  steel  mesh  elements  to  be  displayed  have  their  element 


numbers  listed  on  this  card  (integres) . 


. 
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33.  Number  of  Reinforcement  Elements  Display  Card 

The  only  entry  on  this  card  is  PNREO  (integre) . 


34 .  Reinforcement  Elements  Display  Card 

Similar  to  card  30,  but  concerning  reinforcement  elements. 


35.  Number  of  Nodal  Deflections  Display  Card 

The  only  entry  on  this  card  is  PDTOT. 


36.  Nodal  Deflections  Display  Card 

Similar  to  card  30,  but  with  respect  to  nodal  deflections. 


37.  Printout  Control  Card 


All  data  entries  are  contained  on  one  card. 

Order  of  Entries  Description 


1 

2 

3 

4 

5 

6 
7 


IDEFLN 

IC0N3 

ICONPR 

IMESH 

IREO 

ILOAD 

ITLOAD 


8 


ISTIF 
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The  data  file  is  now  complete.  If  too  few  data  cards  have 
been  included,  the  end  of  file  card  will  be  read,  and  execution  will 
immediately  cease.  If  too  many  cards  have  been  inserted  in  the  data 
file,  the  error  will  be  apparent  in  the  echo  check  printout  of  the 
data. 
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EXPERIMENTAL  PROGRAM  TEST  RESULTS 

Note:  Units  of  kips,  inches,  and  degrees  are  used  throughout. 

1 .  Beam  Dimensions  and  Reinforcement  Details 

Refer  to  Figures  4.3  and  4.4. 

2 .  Loading  Systems 

Three  different  patterns  of  application  of  the  torque  and 
bending  moment  loads  were  employed  in  the  testing  of  the  seven  beams, 
and  are  shown  in  Figure  1-1. 

3 .  Prestress  Levels  and  Shrinkage  Stresses 

Refer  to  Table  4.3  for  prestress  levels. 

Between  transfer  and  testing,  considerable  concrete  shrinkage 
occurred,  stressing  the  longitudinal  conventional  steel  in  compression 
and  surrounding  concrete  in  tension.  Only  the  concrete  in  the  flanges 
is  assumed  to  develop  shrinkage  stresses. 
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Beams 

Shrinkage  Stresses 

Rl 

R2 

R3 

R4  R5 

T1  T2 

Tension  in  Top 

Flange  Concrete 

.115 

.135 

.125 

.095  .12 

.12  .075 

Tension  in  Bottom 
Flange  Concrete 

.203 

.23 

.216 

.117  .2 

.25  .25 

Compression  in 

Top  Reinforcement 

13.5 

15.75 

14.6 

11.1  14. 

15.  9. 

Compression  in 

Bottom  Reinforcement 

24.0 

27.0 

25.5 

19.9  25. 

15.  15. 

TABLE  1-1  SHRINKAGE  STRESSES 


4 .  Test  Loading  and  Deformation  Results 

In  the  following  tables,  the  rotation  values  are  in  degrees 
and  were  measured  over  a  30  inch  central  beam  length.  The  deflection 
measurements  are  the  pure  bending  vertical  displacements  of  the  central 
cross-section. 

The  torque  and  bending  moments  are  those  values  at  the  central 
cross-section.  The  shear  values  are  those  close  to  the  central 
cross-section . 

5.  Reinforcement  Stresses 


Refer  to  Figures  5.6  to  5.12. 
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Load 

Increment 

No. 

Loading 

Deformations 

System 

T 

P 

Bending 

Moment 

Torque 

Shear 

Central 
Def In. 

Rotation 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

2.0 

116 

48 

.024 

.036 

3 

4.0 

216 

96 

.048 

.068 

4 

6.0 

314 

143 

.073 

.0575 

5 

8.0 

414 

192 

.1045 

.08 

6 

10.0 

510 

238 

.139 

.10 

7 

10.78 

552 

259 

.182 

.182 

8 

12.0 

611 

288 

.239 

.295 

9 

12.89 

656 

310 

.29 

.385 

10 

14.04 

713 

337 

.4 

.4776 

11 

A 

14.96 

0.0 

759 

359 

0.0 

.45 

.585 

12 

16.06 

813 

385.4 

.54 

.667 

13 

17.0 

859 

408 

.62 

.768 

14 

17.97 

908 

431 

.69 

.881 

15 

19.0 

959 

456 

.86 

1.427 

16 

19.657 

991 

472 

.97 

1.306 

17 

21.1 

1062 

506 

18 

21.65 

1090 

520 

19 

21.05 

1060 

505 

Note:  Linear  transducers  were  removed  after  load  increment  No.  15. 


TABLE  1-2  BEAK  Rl  TEST  RESULTS 


« 
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Load 

Loading 

Deformations 

Increment 

No. 

System 

T 

P 

Bending 

Moment 

Torque 

Shear 

Central 
Def In. 

Rotation 

1 

0 

0.0 

0.0 

2 

2 

66.4 

.0167 

3 

4 

136 

.03 

4 

6 

206 

.058 

5 

0.0 

8 

276 

0.0 

.078 

0.0 

6 

10 

346 

.097 

7 

12 

417 

.117 

8 

13 

452 

.128 

9 

14 

488 

.14 

10 

1.04 

555 

25 

.163 

.001 

11 

1.92 

600 

46 

.18 

.0032 

12 

3.1 

657 

74 

.21 

.009 

13 

4.03 

703 

97 

.24 

.017 

14 

5.07 

755 

122 

.29 

.043 

15 

A 

6.0 

797 

142 

0.0 

.36 

.104 

16 

7.0 

850 

168 

.45 

.2 

17 

7.9 

896 

190 

.52 

.26 

18 

9.0 

949 

216 

.6 

.33 

19 

10 

998 

240 

.68 

.39 

20 

10.9 

1044 

262 

.77 

.46 

21 

11.3 

1064 

271 

.816 

.52 

22 

12. 

1097 

287 

.91 

.61 

23 

12.19 

1107 

292 

.97 

.68 

24 

12.75 

1135 

306 

25 

13.35 

1165 

321 

26 

13.9 

1192 

334 

27 

14.5 

1222 

348 

28 

14.58 

1226 

350 

29 

15.06 

1249 

361 

30 

15.49 

14.0 

1271 

372 

Note:  Transducers  were  removed  after  load  increment  No.  23. 


TABLE  1-3  BEAM  R2  TEST  RESULTS 
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Load 

Increment 

No. 

Loading 

Deformations 

System 

T 

P 

Bending 

Moment 

Torque 

Shear 

Central 

Def In. 

Rotation 

1 

0.0 

0.0 

0.0 

0.0 

2 

2 

72 

1 

.015 

3 

2 

72 

1 

.016 

4 

4 

145 

2 

.03 

5 

6 

218 

3 

.043 

6 

8 

291 

4 

.054 

7 

0.0 

10 

364 

0.0 

5 

.067 

0.0 

8 

12 

437 

6 

.078 

9 

14 

510 

7 

.09 

10 

16 

583 

8 

.102 

11 

18 

656 

9 

.115 

12 

20 

729 

10 

.13 

13 

22 

802 

11 

.15 

14 

2.04 

854 

49 

.177 

0.0 

15 

4.1 

906 

99 

.205 

.007 

16 

5.8 

947 

138 

.25 

.02 

17 

C 

8.0 

1003 

192 

.325 

.08 

18 

9.8 

1049 

236 

.425 

.16 

19 

11.86 

1099 

285 

.52 

.246 

20 

12.8 

1123 

307 

.59 

.32 

21 

13.9 

1151 

334 

.64 

.41 

22 

14.84 

1174 

356 

.69 

.4 

23 

15.8 

1199 

380 

.77 

.46 

24 

17 

1229 

409 

.85 

.54 

25 

17.3 

1235 

415 

.89 

.575 

26 

17.81 

1248 

428 

.93 

.6 

27 

18.2 

1258 

437 

1.0 

.72 

28 

19.2 

1283 

461 

29 

19.6 

1293 

471 

30 

20.14 

1306 

483 

31 

20.76 

1322 

498 

32 

21.3 

1336 

510 

33 

21.65 

1344 

519 

34 

22.2 

22 

1357 

532 

11 

Note:  Linear  transducers  were  removed  after  load  increment  No.  27. 


TABLE  1-4  BEAM  R3  TEST  RESULTS 
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Load 

Increment 

No. 

Load ing 

Deformations 

System 

T 

P 

Bend ing 
Moment 

Torque 

Shear 

Central 

Defln. 

Rotation 

1 

0.0 

0.0 

0.0 

0.0 

2 

2 

73 

1 

.01 

3 

4 

146 

2 

.023 

4 

6 

220 

3 

.036 

5 

0.0 

8 

293 

0.0 

4 

.05 

0.0 

6 

10 

366 

5 

.063 

7 

12 

440 

6 

.076 

8 

14 

513 

7 

.089 

9 

16 

587 

8 

.102 

10 

18 

660 

9 

.116 

11 

20 

733 

10 

.131 

12 

4.04 

831 

97 

.167 

13 

6.1 

883 

147 

.197 

Similar 

14 

C 

7.9 

928 

190 

.244 

to 

15 

9.84 

976 

236 

.33 

R3* 

16 

12.1 

1032 

290 

.45 

17 

13.95 

1079 

335 

.55 

18 

22 

1152 

11 

.61 

19 

24 

1222 

12 

.6725 

20 

26 

1297 

13 

.78 

21 

28 

1366 

14 

.90 

22 

30 

1403 

15 

.99 

23 

31 

1443 

15.5 

1.09 

24 

32 

1484 

16 

1.28 

25 

33 

1511 

16.5 

1.51 

26 

34 

1554 

17 

1.665 

27 

35 

1591 

17.5 

2.02 

28 

1 

13.95 

36 

1624 

335 

18 

*  Rotation  results  for  this  short  torsion  loading  interval  are  similar  to  those 
corresponding  results  for  beam  R3. 


Note:  Linear  transducers  were  removed  after  load  increment  No.  27. 


TABLE  1-5  BEAM  R4  TEST  RESULTS 
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Load 

Loading 

Deforti 

nations 

Increment 

No. 

System 

T 

P 

Bending 

Moment 

Torque 

Shear 

Central 

Defln. 

Rotation 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

1.96 

49 

47 

.02 

.006 

3 

3.92 

98 

94 

.043 

.009 

4 

6.15 

154 

148 

.068 

.067 

5 

8 

200 

192 

.085 

.0125 

6 

9 

250 

240 

.103 

.0325 

7 

11 

276 

265 

.113 

.036 

8 

11.9 

298 

282 

.122 

.043 

9 

13 

324 

311 

.131 

.063 

10 

14.2 

355 

341 

.173 

.156 

11 

14.84 

371 

356 

.196 

.21 

12 

15.5 

388 

372 

.217 

.25 

13 

15.9 

398 

382 

.233 

.272 

14 

16.7 

418 

401 

.253 

.31 

15 

17.2 

431 

414 

.273 

.34 

16 

B 

17.4 

0.0 

435 

417 

0.0 

.294 

.39 

17 

18 

451 

433 

.307 

.42 

18 

18.5 

463 

444 

.327 

.437 

19 

19 

476 

457 

.354 

.49 

20 

19.5 

487 

468 

.383 

.53 

21 

19.9 

497 

477 

.4 

.57 

22 

20.5 

513 

492 

.426 

.61 

23 

20.93 

523 

502 

.44 

.64 

24 

21.6 

540 

519 

.47 

.69 

25 

21.82 

546 

524 

.486 

.73 

26 

22.3 

558 

536 

.504 

.76 

27 

22.84 

571 

548 

.53 

.8 

28 

23.3 

583 

559 

.55 

.84 

29 

23.87 

597 

573 

.583 

.89 

30 

24.38 

609 

585 

.617 

.97 

31 

25 

626 

601 

.65 

1.06 

32 

25.6 

640 

615 

.683 

1.167 

TABLE  1-6  BEAM  R5  TEST  RESULTS 
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Beam 

Load 

Loading 

Deformations 

Increment 

No. 

System 

T 

P 

Bending 

Moment 

Torque 

Shear 

Central 

Def In. 

Rotation 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

1.19 

77 

28.6 

.0225 

.018 

3 

2.03 

118 

49 

.04 

.031 

4 

2.95 

164 

71 

.062 

.032 

5 

4.0 

215 

97 

.085 

.008 

6 

5.08 

269 

122 

.109 

.02 

T1 

7 

A 

6.07 

0.0 

319 

146 

0.0 

.14 

.052 

8 

6.75 

352 

162 

.215 

.067 

9 

8.13 

420 

195 

.316 

.13 

10 

10.0 

511 

239 

.565 

.38 

11 

11.0 

563 

264 

.68 

.53 

12 

11.3 

577 

271 

.743 

.63 

13 

9.5 

490 

229 

1 

0.0 

0.0 

0.0 

2 

2 

72 

.029 

3 

4 

145 

.063 

4 

0.0 

6 

216 

0.0 

.1 

0.0 

5 

7 

252 

.12 

6 

8 

289 

.143 

7 

9 

322 

.166 

8 

10 

361 

.197 

9 

11 

396 

.24 

10 

.98 

445 

24 

.32 

- 

T2 

11 

A 

2 

495 

48 

0.0 

.38 

- 

12 

3 

546 

73 

.49 

.1 

13 

4 

594 

96 

.68 

.035 

14 

4.53 

620 

109 

.74 

.156 

15 

5 

647 

121 

.8 

.225 

16 

5.45 

666 

131 

.86 

.26 

17 

6 

690 

142 

.955 

.286 

18 

6.43 

715 

154 

1.06 

.4 

19 

6.9 

738 

166 

1.25 

.528 

20 

7.4 

763 

178 

1.53 

.79 

21 

7.94 

789 

191 

22 

8.2 

801 

196.5 

23 

8.13 

11 

798 

195 

Notes:  Linear  transducers  were  removed  after  load  increment  13  for  Tl,  and  increment  20 
for  beam  T2 


TABLE  1-7  BEAMS  T1  AND  T2  TEST  RESULTS 


(a)  Loading  system  A 
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